AD  A  1 24  702  DECENTRALIZED  CONTROL  OF  A  LARGE  SPACE  STRUCTURE  AS  1/flL 

APPLIED  TO  THE  CSOL  2  MODE  L ( U )  AIR  FORCE  INST  OF  TECH 
WRIGHT-PATTERSON  AF8  OH  SCHOOL  OF  ENGI . .  E  S  ALDRIDGE 
UNCLASSIFIED  OEC  82  AF I T/GA/AA/82D -  1  F/G  77/2  NL 


MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  (UNtAU  Of  STANDARDS  1663  » 


ADA  13  4  7  v* 


AFIT/GA/AA/82D-1 


DECENTRALIZED  CONTROL 
OF  A  LARGE  SPACE  STRUCTURE 
AS  APPLIED  TO  THE  CSDL  2  MODEL 

THESIS 


AFIT/GA/AA/82D-1  Edward  S.  Aldridge 

2d  Lt  USAF 


Approved  for  public  release;  distribution  unlimited 


AFIT/GA/AA/82D-1 


DECENTRALIZED  CONTROL 
OF  A  LARGE  SPACE  STRUCTURE 
AS  APPLIED  TO  THE  CSDL  2  MODEL 


THESIS 


Presented  to  the  Faculty  of  the  School  of  Engineering 
of  the  Air  Force  Institute  of  Technology 
Air  University 

in  Partial  Fulfillment  of  the 
Requirements  for  the  Degree  of  _ 


Graduate  Astronautical  Engineering 
December  1982 


Approved  for  public  release;  distribution  unlimited 


Preface 


To  ay  thesis  advisor.  Dr.  R.  A.  Calico,  1  aa  grateful  for  your 
guidance  and  patience  as  advisor  and  your  knowledge  and  insight  as  in¬ 
structor.  Tour  support  in  the  course  of  this  work  was  invaluable.  I 
hope  the  results  contained  herein  aay  be  of  value  to  you. 

To  ay  typist,  Shirley,  thank  you  for  your  saile  and  your  typing 
services,  the  latter  of  which  aade  this  final  product  possible. 

Finally,  to  ay  GA  and  GAE  classmates,  thanks  for  your  friendship, 
support  and  empathy.  You're  the  ones  who  aade  AF1T  a  worthwhile 
experience. 

Edward  S.  Aldridge 


11 


Contents 


Preface . 11 

List  of  Figures . v 

List  of  Tables . vi 

Abstract . 

I.  Introduction  .  1 

II.  Model  Configurations  .  5 

CSDL  1  Model .  5 

CSDL  2  Model .  9 

III.  System  Model . 21 

Equations  of  Motion  ......  . 21 

Control  Model  . . 23 

Modal  Control . 27 

N  Controllers  .  3° 

Three  Controllers . 32 

Four  Controllers . 37 

Sensor/Actuator  Requirements  .  39 

IV.  Transformation  Technique  .  42 

V.  Computer  Model . . 

VI.  Investigation . . 

Outline . 54 

Modal  Selection  and  Grouping . 55 

CSDL  1 . 56 

CSDL  . . 56 

VII.  Results . 60 

VIII.  Conclusions . 76 

IX.  Recommendations . 78 

Bibliography  .  79 

Appendix  A:  CSDL  1  NASTRAH  Analysis . 80 

ill 


Contents 


Appendix  B; 
Appendix  C 
Appendix  D 
Appendix  E 
Vita  .  .  . 


CSDL  2  NASTRAN  Analysis  . 

Main  Prograa  Listing  . 

Program  Subroutines  . 

Line-of-Sight  and  Defocus  Algorithm 


U 


figure 

1 


CSDL  1  Syatea  Modal .  7 

2  CSDL  2  Conceptual  Diagram . 10 

3  CSDL  2  Syatea  Model . 11 

4  CSDL  2  Optical  Syataa  . . 13 

5  CSDL  2  Support  Truaa  Dimensions . 14 

6  CSDL  2  Support  Structure  Nodal  Placement  .  15 


table 

I 

II 

III 

IV 

V 

VI 

VII 

VIII 

n 

X 

XI 

XII 

XIII 

XIV 

XV 

XVI 

XVII 

XVIII 


List  of  Tables 

*S8i 


CSDL  Hod*  Coordinates  .  6 

Kay  Results  of  NAST8AH  Eigenvalue  Analysis  on  Nominal 
CSDL  1  Modal . 8 

.  Initial  Conditions  Applied  to  CSDL  1 . 9 

CSDL  2  Mode  Coordinates . 16 

Luaped  Mass  Location  and  Distribution . 17 

Key  Results  of  HASTRAN  Eigenvalue  Analysis  on  Nominal 
CSDL  1  Model . <18 

CSDL  2  Sensor /Actuator  Locations  and  Orientations  .  .  *20 


Spillover  Elimination  for  Decoupling  Three  Controllers  .  35 
Spillover  Elimination  for  Decoupling  Four  Controllers  .  39 


Relative  Angles  Between  CSDL  1  Modes  .  57 

Relative  Angles  Between  CSDL  2  Modes  . . 58 

CSDL  1  Overall  Eigenvalue  Analysis  - 

3  Controllers:  8  Modes . 62 

Upper  and  Lower  Triangular  Transformation  Comparison  of 

Controller  Eigenvalues  .  63 

CSDL  1  Overall  Eigenvalue  Analysis  -  3  Controllers: 

12  Modes . 65 

CSDL  1  Overall  Eigenvalue  Analysis  -  4  Controllers: 

12  Modes . . . 66 

CSDL  2  Overall  Eigenvalue  Analysis  -  3  Controllers: 

q±  -  20 . 69 


CSDL  2  Overall  Eigenvalue  Analysis  -  3  Controllers: 

'rigid  "  2*  'flexible  "  1000  . 

CSDL  2  Controller  Eigenvalue  Analysis  -  3  Controllers: 

'rigid  “  2*  'flexible  “  1000  . 

CSDL  2  Overall  Eigenvalue  Analysis  -  4  Controllers: 

'rigid  "  2*  'flexible  “  500  . 


vi 


XIX 


74 


Abstract 


Modern  optimal  control  methods  ere  used  to  develop  a  multiple  input 
multiple  output  controller.  The  controller  is  then  applied  to  two  models. 
The  first  model  is  a  lumped  mass  model  of  a  tetrahedron  consisting  of  four 
unit  masses  interconnected  by  isotropic  massless  rods.  These  rods  are 
assumed  to  be  pin-connected  and  may  undergo  axial  deformation  only.  The 
second  model  is  a  sophisticated  optical  space  structure  more  representa¬ 
tive  of  large  flexible  space  structures  than  the  first  model.  This  model 
consists  of  fifty-nine  nodes  and  twenty-three  lumped  masses.  The  beam 
elements  are  fully  connected  and  may  support  axial,  transverse  and  torsional 
deformations.  NASTKA8  is  employed  to  generate  modal  approximations  of 
both  models,  as  well  as  the  mode  shapes  and  frequencies  of  the  resulting 
modes.  Twelve  modes  are  generated  for  the  first  model.  Of  the  numerous 
modes  available  for  the  second  model,  only  the  first  forty-four  modes  are 
addressed,  and  of  these  twelve  are  implemented  in  the  controller. 

The  control  problem  is  formulated  in  state  vector  form  and  full  state 
feedback  is  implemented.  The  state  is  represented  as  modal  amplitudes  and 
rates  and  the  feedback  gains  are  generated  using  steady  state  optimal  regu¬ 
lator  theory.  State  estimates  are  provided  by  means  of  a  deterministic 
observer.  System  outputs  are  obtained  by  position  sensors  sad  control  Is 
applied  by  point  force  actuators.  The  technique  by  which  "spillover"  is 
eliminated  is  developed  using  the  method  of  singular  value  decomposition. 

Decentralised  control  was  accomplished  using  three  sad  four  controllers 
with  both  models.  Conditions  for  which  tbs  stability  of  each  model  la 
assured  are  developed.  Model  oae  la  run  with  three  controllers 
the  first  eight  modes  to  verify  system  stability.  The  remaining  four 
modes  are  added  as  residuals  to  the  three  controllers.  Model  one  was  also 
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DECENTRALIZED  CONTROL  OF  A  LARGE  SPACE 
STRUCTURE  AS  APPLIED  TO  THE  CSDL  2  MODEL 

I.  Introduction 

With  the  current  success  in  the  Space  Transportation  System,  the  near 
future  holds  a  significant  increase  in  the  size  of  structures  that  may  be 
employed  in  space.  Projected  dimensions  for  the  large  space  structures 
range  from  tens  to  thousands  of  meters  in  size,  taking  advantage  of  the 
low  gravity  environment  to  make  the  systems  cost  effective.  Tubular, 
lightweight  truss  members  make  these  structures  practical,  but  also  make 
them  very  flexible.  The  increases  in  size  and  flexibility  lead  to  an  over 
all  increase  in  the  number  of  low  frequency  modes  that  may  be  contained 
within  the  control  system  bandwidth.  Control  of  such  structures  then  be¬ 
comes  Increasingly  difficult  as  the  dimensions  of  the  controllers  increase 
Since  active  control  is  performed  by  on-line  computers,  larger  controller 

dimensions  result  in  slower  control  response.  To  keep  the  controllers 

.  ** 

dimensionally  realistic,  modeling  of  the  satellite  and  its  structural 
modes  becomes  a  prime  concern.  A  discrete  structural  analysis  of  a  large 
space  structure  may  Include  from  one  hundred  to  several  hundred  of  these 
modes.  Unfortunately,  the  accuracy  of  the  modal  information  obtained  from 
such  an  analysis  decreases  with  increasing  mode  number.  Such  modeling 
Inaccuracies  could  result  in  overall  system  instability  when  not  properly 
accounted  for. 

Of  the  various  control  techniques  available,  modern  state-space  con¬ 
trol  theory  appears  to  be  best  suited  for  application  to  large  flexible 
space  structures,  in  light  of  the  problems  with  off-line  computing 
accuracy  and  on-line  computing  speed.  These  stats  specs  controllers  make 


use  of  reduced  order  finite  element  structural  models  to  minimize  compu¬ 
tational  burdens.  This  may  generate  modeling  and  reduction  errors,  but 
allows  this  method  to  be  easily  applied  to  any  of  a  wide  variety  of  large 
flexible  space  structures.  As  implied  earlier,  the  number  of  structural 

modes  any  single  controller  can  handle  is  limited  by  computational  consi- 

% 

derations,  but  these  limits  may  be  extended  by  using  multiple  controllers 
within  the  system,  each  controller  performing  independently. 

Even  with  an  expanded  number  of  controllers,  the  number  of  modes  that 
may  be  controlled  is  small  compared  to  the  number  of  modes  that  exist  for 
a  given  structure.  Therefore,  the  selection  of  modes  to  be  controlled  must 
be  made  carefully.  Only  those  modes  affecting  performance  need  be  con¬ 
trolled.  The  terms  "controlled"  and  "critical"  will  be  used  interchange¬ 
ably  for  the  modes.  The  remaining  uncontrolled  modes  fall  into  three 
categories:  suppressed,  residual  and  unmodeled.  Obviously,  for  a  large 
space  structure,  the  number  of  structural  modes  approaches  infinity.  Natural 
damping  in  the  structure  will  normally  prevent  instabilities  arising  from 
the  higher  frequency  modes,  so,  for  model  simplicity,  these  are  truncated 
and  left  as  unmodeled  modes.  The  remaining  uncontrolled  modes  are  modeled 
modes.  Of  these,  some  may  have  destabilizing  affects  due  to  spillover 
and  therefore  have  to  be  made  transparent  to  the  controller.  These  are 
called  suppressed  modes,  the  last  mode  group  is  modeled,  uncontrolled  and 
unsuppressed.  These  are  the  residual  modes  and  may  move  freely  when  con¬ 
trol  is  applied.  They  may  become  more  stable,  less  stable  or  unstable  due 
to  control  and  observation  spillover  from  the  critical  modes. 

It  must  be  remembered  that  even  though  specific  modes  are  actively 
controlled,  the  residual  and  unmodeled  modes  still  exist  and  will  contami¬ 
nate  the  observation  (sensor)  data.  Therefore,  the  controlled  and  uncon- 
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trolled  modes  are  coupled.  Balas  (Ref  1)  calls  this  coupling  "observation 
spillover".  Likewise,  any  control  applied  to  the  critical  modes  may  excite 
one  or  several  uncontrolled  modes.  This  form  of  coupling  is  referred  to  as 
"control  spillover",  and  Balas  goes  on  to  state  that  either  or  both  types 
of  spillover  may  result  in  overall  instabilities.  He  proposes  a  state 
variable  feedback  controller  which  relies  on  narrow  bandpass  filters 
to  eliminate  observation  spillover  by  filtering  out  suppressed  mode  frequen¬ 
cies  from  controller  input  data. 

Sesak  (Ref  2)  proposed  the  use  of  a  singular  perturbation  technique  to 
develop  an  appropriate  feedback  controller  and  eliminate  instabilities  due 
to  spillover.  Coradetti  (Ref  3)  later  expanded  Sesak' s  approach  and  con¬ 
cluded  that,  in  the  limiting  sense  with  an  infinite  penalty  against  any 
spillover,  the  singular  perturbation  technique  is  equivalent  to  finding  a 
transformation  matrix  which,  when  applied  to  the  feedback  gains,  will  drive 
the  spillover  terms  to  zero.  This  transformation  matrix  is  determined  by 
performing  a  singular  value  decomposition  of  the  control  and  observation 
matrices  (Ref  4) .  When  the  transformation  technique  is  coupled  with  the 
modem  state-space  control  technique,  an  effective  method  is  obtained  for 
eliminating  spillover,  and  works  equally  well  on  control  and  observation 
spillover.  Moreover,  even  when  such  spillover  is  not  detrimental  to  the 
overall  system  stability,  its  elimination  can  only  enhance  the  system 
performance. 

The  intent  of  this  thesis  is  to  apply  the  aforementioned  control 
techniques  in  developing  a  control  system  consisting  of  three  or  more 
decentralized  controllers.  This  control  system  will  be  applied  to  a  lumped 
mass  tetrahedron  model  generated  by  the  Charles  Stark  Draper  Laboratory, 
Inc.,  (CSDL),  hereafter  called  the  CSDL  1  model.  Calico  and  Janlszewski 
(Ref  5)  applied  the  described  technique  to  the  CSDL  1  model  using  a  single 
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controller  end  eliminating  only  observation  spillover.  Later,  Calico  and 
Miller  (Ref  6)  expanded  the  system  to  a  dual  controller,  and  showed  that 
only  observation  spillover  elimination  is  not  sufficient  for  higher  order 
controllers.  More  detailed  results  are  given  by  Miller  (Ref  7)  for  the 
dual  controller  case.  Therefore,  it  has  been  demonstrated  that  this  tech¬ 
nique  is  appropriate  for  this  model.  This  present  study  will  use  three  con¬ 
trollers  for  this  model.  System  performance  .111  be  evaluated  by  eigen¬ 
value  analysis  of  the  closed  loop  system.  Next,  the  triple  controller 
system  will  be  expanded  to  accomodate  the  second  CSDL  model— a  three- 
mirror,  optical  space  system.  Again,  an  eigenvalue  analysis  will  determine 
the  control  system's  closed  loop  performance.  Line  of  sight  pointing 
accuracy  and  defocus  are  performance  criteria  that  are  mentioned  for 
information,  but  will  not  be  addressed  in  this  investigation.  In  applying 
the  control  method,  for  both  models,  position  sensors  are  used  to  determine 
modal  amplitudes  while  point  force  actuators  provide  the  state  variable 
feedback  control. 

The  following  sections  will  detail  both  of  the  CSDL  models  and  their 
finite  element  representations.  Afterward,  the  modal  control  and  matrix 
transformation  methods  will  be  discussed.  Finally,  the  computer  program 
implementation  and  results  will  be  presented. 
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II.  Model  Configuration 

Illustration  and  demonstration  of  controller  design  methods  for  large 
space  structures  has  always  been  a  difficult  problem.  The  very  nature  of 
large  space  structures  prevents  the  development  of  a  simple  textbook  exam¬ 
ple.  In  response  to  this  problem,  the  Charles  Stark  Draper  Laboratory 
(CSDL),  Inc.,  of  Cambridge,  Massachusetts,  developed  two  (paper)  models  for 
research.  Both  models  are  used  in  this  investigation.  Finite  element 
representations  of  the  models  are  generated  by  the  NASTRAN  computer  pro¬ 
grams.  Presentations  of  the  models  and  their  eigenvalue  analyses  follow. 
CSDL  1  Model 

The  first  model  used  is  a  lumped  mass  tetrahedron  and  is  referred  to 
as  the  CSDL  1  model.  This  model  was  selected  for  its  simplicity  as  well 
as  its  similarity  to  basic  large  space  structures  under  consideration, 
from  both  a  structures  and  a  control  point  of  view.  The  tetiahedral  struc¬ 
ture  is  the  building  block  of  most  large  space  structure  design  concepts. 

It  provides  a  low  order  model  to  which  control  systems  may  be  easily 
applied  due  to  the  small  number  of  modes  present.  Also,  response  charac¬ 
teristics  exhibited  by  the  model  are  very  similar  to  those  observed  in 
large  space  structures.  This  is  probably  the  simplest  model  available 
which  behaves  much  like  a  large  space  structure. 

The  finite  element  model  of  the  structure  is  depicted  in  Fig  1.  The 
structure  has  twelve  members  joined  at  ten  nodes.  The  truss  members  are 
considered  massless  and  are  pin-connected  at  the  nodes,  so  that  only  axial 
forces  are  exerted  (no  bending  moments).  The  masses  are  equal— one  unit 
each— and  are  located  at  the  first  four  nodes  (the  vertices  of  the  tetrs- 
hedron  proper).  Since  each  mass  is  asaiaed  to  have  three  translational 
degrees  of  freedom,  the  system  has  twelve  structural  modes.  The  last  six 
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nodes  fora  the  ground  connections  for  the  three  right-angled  bipods  which 
support  the  tetrahedral  truss.  This  ground  base  provides  a  reference  from 
which  a  line  of  sight  nay  be  established.  The  coordinates  for  the  ten  nodes 
are  given  in  Table  I.  The  reference  frame  origin  for  the  coordinates  is 
placed  directly  below  the  apex  in  the  plane  of  nodes  five  through  ten.  Six 
pair  of  collocated  force  actuators  and  position  sensors  are  used  on  this 
model. 

Table  I 

CSDL  1  Node  Coordinates 


Node 

X 

Z 

£ 

1 

0.0 

0.0 

10.165 

2 

-5.0 

-2.887 

2.0 

3 

5.0 

2.887 

2.0 

4 

0.0 

5.7735 

2.0 

5 

-6.0 

-1.1547 

0.0 

6 

-4.0 

-4.6188 

0.0 

7 

4.0 

-4.6188 

0.0 

8 

6.0 

-1.1547 

0.0 

9 

2.0 

5.7735 

0.0 

10 

-2.0 

5.7735 

0.0 

The  key  results  of  an  eigenvalue  analysis  of  this  model  are  presented 
in  Table  II.  Listed  are  the  generalised  mass  and  stiffness  and  the  natural 
frequencies  of  each  structural  mode.  The  eigenvectors  for  each  mode  ere 
presented  in  Appendix  A.  This  data  was  obtained  from  a  NAS TRAN  eigenvalue 
analysis. 

e  - 
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Table  II 

Key  Result*  of  NAS TRAN  Eigenvalue  Analysis  on  CSDL  1  Model 


Mode 

Generalized 

Mass 

Generalized 

Stiffness 

rad 

u  sec 

IS^2 

£J  sec 

1 

1. *01+00 

1.37E+00 

1.17E400 

1.37E400 

2 

1.0E+00 

2.15E+00 

1.47E400 

2.1SE400 

3 

1.0E+00 

8.79E+00 

2.97E400 

8.79E400 

4 

1.0E+00 

1.26E401 

3.56E+00 

1.26E+01 

5 

1.0E+00 

1.48E401 

3.85E+00 

1.48E+01 

6 

1.0E+00 

2.63E401 

5.15E400 

2.65E+01 

7 

1.0E+00 

3.22E+01 

5.68E+00 

3.22E401 

8 

1.01400 

3.26E+01 

5.71E+00 

3.26E+01 

9 

1.01400 

7.99E+01 

8.94E+00 

7.99E401 

10 

1.01400 

1.06E+02 

1.01E+01 

1.06E402 

11 

1.01400 

1.19E402 

1.09E401 

1.19E+02 

12 

1.01400 

1.95E402 

1 . 40E+01 

1.9SE402 

Eigenvalue  analysis  of  the  nodal  movements  will  give  an  indication  of 
the  performance  of  the  system  control.  Line  of  sight,  based  on  the  x-y 
motion  of  mode  one,  is  an  important  performance  parameter.  Initial  condi¬ 
tion*  may  be  applied  to  the  model  to  develop  a  time  history  of  the  system 
response.  The  initial  conditions  for  this  model  are  listed  in  Table  III. 
As  a  first  cut  on  the  three  controller  design,  the  line  of  sight  pointing 
performance  will  not  be  addressed  in  this  investigation.  Nevertheless, 
the  development  of  the  error  terms  to  which  the  initial  conditions  are 
applied  will  be  explained  along  with  the  system  equations  of  motion. 


8 


Table  III 


Initial  Conditions  Applied ‘to  CSDL  1 


Displaci 


1 

-.001 

-.003 

2 

0.006 

0.010 

3 

0.001 

0.030 

4 

-.009 

-.020 

5 

0.008 

0.020 

6 

-.001 

-.020 

7 

-.002 

-.003 

8 

0.002 

0.004 

9 

0.000 

0.000 

10 

0.000 

0.000 

11 

0.000 

0.000 

12 

0.000 

0.000 

CSDL  2  Modal 

Tha  second  aodel  under  consideration  is  a  "generic”  optical  space 
structure  which  has  a  behavior  ouch  closer  to  that  of  a  large  space 
structure  than  the  first  discussed.  This  aodel  is  referred  to  as  CSDL  2. 
Figure  2  shows  a  conceptual  view  of  the  structure  and  Fig  3  is  a  finite 
eleaent  representation  of  the  aodel. 

CSDL  2  is  a  non-trivial  aodel  representing  a  wide-angle,  three-wirror , 
optical  space  systea  The  two  aajor  coaponents  of  the  systea  are  the  opti¬ 
cal  support  struct  and  the  equlpaant  section.  The  optical  support  struc¬ 
ture  consists  of  the  upper  alrror  support  truss,  the  lower  airror  support 
trus,  and  the  aster ing  truss.  The  upper  airror  support  truss  contains  the 
prlaary  airror  (convex  shaped)  and  the  tertiary  airror  (concave  shaped). 

The  lower  airror  support  truss  contains  the  secondary  airror  (flat)  and  the 
focal  plane  (laage  receiving  device).  The  as  ter ing  truss  aaintains  the 

a 

airror  separation  and  is  the  key  section  w&dn  exaaining  defocus.  The  opti- 
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Figure  3.  CSDL  System  Model 


cal  support  structure  and  mirror  placement  arm  shown  in  Figure  4.  Attached 
to  the  lower  side  of  the  lower  mirror  support  truss  Is  the  'equipment  sec¬ 
tion  which  consists  of  the  control  package,  modeled  as  a  rigid  body,  with 
two  cantilevered  flexible  solar  panels.  The  full  structure  Is  approximately 
twenty-eight  meters  high  and  has  a  mass  of  9300  kilograms.  The  structural 
dimensions  are  shown  in  Figure  5. 

The  finite  element  model  of  the  structure  contains  fifty-nine  node 
points,  but  the  actual  structure  has  only  fifty-one  nodes.  The  extra  nodes 
were  added  to  provide  more  detail  in  the  modeling  of  the  mirrors  and  equip¬ 
ment  section  (Kef  8).  The  coordinates  of  the  nodes  are  given  in  Table  IV, 
and  the  placement  of  the  nodes  In  the  support  structure  are  shown  In  Fig 
6.  Unlike  the  first  model,  the  truss  members  are  fully  joined  so  that  bend¬ 
ing  and  torsion  are  allowed.  The  truas  elements  are  made  of  graphite-epoxy 
and  assumed  to  be  massless.  The  system  mass  Is  limited  at  twenty-three  nodes 
and  distributed  as  shown  In  Table  V.  The  largest  mass  Is  located  In  the 
equipment  package,  as  would  be  expected. 

The  key  results  of  an  eigenvalue  analysis  performed  on  this  model  are 
listed  In  Table  VI.  The  generalised  mass  and  stiffness,  as  well  as  the 
natural  frequency  of  the  first  forty-four  structural  nodes  Is  given.  Of 
the  modes  listed,  those  with  an  asterisk  were  used  In  this  study.  The 
associated  eigenvectors  for  each  mode  are  presented  In  Appendix  B.  Again, 
this  data  Is  obtained  via  the  KASTIAff  computer  program. 

This  model  makes  use  of  twenty-one  pairs  of  collocated  force  actuators 
ami  position  sensors.  A  list  of  sensor /actuator  locat loos  and  orientations 
is  provided  in  Table  VII. 

As  in  the  C3DL  I  model,  eigenvalue  analyses  will  provide  control  per¬ 
formance  information.  Liae-of-slgfat  performance,  as  well  as  defocus  along 
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Table  VI 


r  Results 

of  NASTRAN  Eigenvalue  Analysis 

on  Nominal 

CSDL  2  Mode: 

Generalized 

Generalized 

rad 

.2 

rad^ 

Mode 

Mass 

Stiffness 

u>  sec 

0  sec 

1-6* 

l. OOE+OO 

0.0 

0.0 

0.0 

7* 

1.00E+00 

5. 128E-01 

7.161E-01 

5 . 128E-01 

8 

1.00E+00 

8.521E-01 

9. 231E-01 

8.521E-01 

9 

1. OOE+OO 

8.835E-01 

9.399E-01 

8.835E-01 

10 

1.00E+00 

1. 212E+00 

1. 101E+00 

1. 212E+00 

11 

1. OOE+OO 

8. 189E+00 

2.862E+00 

8.189E+00 

12* 

1.00E+00 

1. 266E+01 

3 . 502E+00 

1.226E+01 

13* 

1. OOE+OO 

1.403E+01 

3. 746E+00 

1 . 403E+01 

14 

1.00E+00 

1.492E+01 

3 . 863E+00 

1.492E+01 

15 

1. OOE+OO 

' .599E+01 

3.998E+O0 

1.599E+01 

16 

1. OOE+OO 

1.625E+01 

4.032E+00 

1 . 625E+01 

17* 

1. 00E+00 

2 . 623E+01 

5. 122E+00 

2.623E+01 

18 

1. 00E+00 

2 . 630E+01 

5.128E+00 

2 . 630E+01 

19 

1. OOE+OO 

2.677E+01 

5.174E+00 

2.677E+01 

20 

1. 00E+00 

3. 310E+01 

5.753E+00 

3. 310E+01 

21* 

1. OOE+OO 

3.730E+01 

6.197E+0 0 

3 . 7  30E+01 

22* 

1. OOE+OO 

5 . 301E+01 

7. 281E+00 

5 . 301E+01 

23 

1. OOE+OO 

9 . 498E+01 

9. 746E+00 

9.498E+01 

24* 

1. OOE+OO 

1. 241E+02 

1 . 114E+01 

1.241E+02 

25 

1. OOE+OO 

1. 999E+02 

1.414E+01 

1 . 999E+02 

26 

1. OOE+OO 

2 . 001E+02 

1.416E+01 

2.001E+02 

27 

1. OOE+OO 

4. 654E+02 

2. 157E+01 

4 . 654E+02 

28* 

1. OOE+OO 

4. 705E+02 

2.169E+01 

4.705E+02 

29 

1 . GOE+OO 

6. 182E+02 

2.468E+01 

6 . 182E+02 

30* 

1. OOE+OO 

6. 275E+02 

2 . 505E+01 

6.275E+02 

31 

1. OOE+OO 

6. 481E+02 

2 . 546E+01 

6 . 481E+02 

32 

1. OOE+OO 

7.428E+02 

2.725E+01 

7.428E+02 

33 

1.  OOE+OO 

1.  700E+03 

4.123E+01 

1.700E+03 

34 

1. OOE+OO 

2. 568E+03 

5.067E+01 

2.568  .+03 

35 

1. OOE+OO 

2.821E+03 

5. 311E+01 

2.821E+03 

36 

1. OOE+OO 

3.095E+03 

5. 563E+01 

3.095E+03 
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Table  VI.  continued 


Key  Results  of  NASTRAN  Eigenvalue  Analysis  on  Nominal  CSDL  2  Model 


Mode 


Generalized 

Mass 


Generalized 

Stiffness 


rad 

u  sec 


n 


sec 


37 

1.00E400 

3. 205E+03 

5.661E+01 

3 . 205E+03 

38 

LOOE+OO 

4. 221E+03 

6.497E+01 

4.221E+03 

39 

1.00E+00 

4. 380E+03 

6. 618E+01 

4. 380E+03 

40 

1.00E400 

5 . 266E+03 

7.257E401 

5.266E+03 

41 

1.00E400 

5. 358E+03 

7 . 320E+01 

5 . 358E+03 

42 

1.00E400 

5. 360E+03 

7.321E+01 

5 . 360E+03 

43 

1.00E+00 

5. 361E+03 

7.322E+01 

5.361E403 

44 

1.00E+00 

5. 368E+03 

7. 327E+01 

5 . 368E+03 

*denotes  modes  actively  controlled  in  this  study.  Of  the  rigid 
body  modes,  only  4,  S,  and  6  were  controlled. 


Table  VII 


CSDL  2  Sensor /Actuator  Locations  and  Orientations 


Pair 

Node 

X 

L 

£ 

1 

9 

0 

1 

0  (In 

2 

9 

0 

0 

1  direction 

3 

10 

0 

0 

1  cosines) 

4 

11 

1 

0 

0 

5 

11 

0 

1 

0 

6 

11 

0 

0 

1 

7 

12 

0 

0 

1 

8 

27 

1 

0 

0 

9 

27 

0 

1 

0 

10 

27 

0 

0 

1 

11 

28 

0 

0 

1 

12 

29 

0 

1 

0 

13 

29 

0 

0 

1 

14 

30 

0 

0 

1 

15 

32 

0 

0 

1 

16 

33 

0 

0 

1 

17 

34 

1 

0 

0 

18 

34 

0 

1 

0 

19 

34 

0 

0 

1 

20 

35 

0 

1 

0 

21 

35 

0 

0 

1 

z-axis , 

are  important  analysis  factors. 

but  as 

a  first  cut,  will  not 

be  directly  addressed  in  this  Investigation.  Instead,  the  eigenvalue 
analysis  will  provide  information  on  controller  maintenance  of  modal 
stability  and  on  controller  Independence  (decoupling) . 

The  controller  development  on  which  this  study  is  based  will  now  be 


presented 


III.  System  Model 


ImtioM  of  Hocloc 

k »  prtMottd  by  Calico  and  Miller  (baft  6  and  7),  the  system  aodel 
aay  b«  developed  from  tl«  vibrational  aquation*  of  notion  for  a  largo 
spac*  structure  given  generally  as 

Mg  ♦  Eg  ♦  Kg  ■  Du  (1) 

where 

M  ■  nan  synnetric  ness  aatrix 
E  “  nxn  synnetric  danping  matrix 
K  *  nxn  symetrlc  stiffness  matrix 

D  ■  nxn  matrix  of  nodal,  attitude- evaluated  actuator  locations 
g  *  nxl  generalised  coordinate  vector 
u  “  nxl  control  input  vector 
Introducing  the  nxn  nodal  matrix  *  for  Eq  1,  such  that 

g  ■  *  n  (2) 

where  n  is  the  n-vector  of  nodal  coordinates,  Eq  1  nay  be  written  as 


\ 

I 

n  + 

\ 

2(u 

n  + 

\ 

2 

w 

\ 

\ 

\ 

^  «■! 

T  _ 
♦  Du 


(3) 


the  and  t« 


being  natural  frequencies  and  danping  coefficients, 
respectively,  of  the  specific  nodes.  The  properties  of  the  nodal  matrix 
♦  are  such  that  the  coefficients  of  Eq  3  are  given  by 

Tv 


\ 


'ii 


2C*» 


-  A* 


*  #te# 


•  ♦**♦ 


(4) 
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F 


[■ 


2  Cw 


x 

* 
i 

m 

"  2 
<*> , 


whara 

R  1 

turn  identity  matrix 
-  nxn  diagonal  daaping  aatrix 
■  nxn  diagonal  aacrix  of  aiganvaluaa  of  Eq  1 
Equation  3  aa y  now  b*  convartad  into  a  acata  spaca  ropraaancacion  of 
Cha  systaa,  givan  by 

a 

x  -  Ax  +  Bu  (5) 

in  which 

A  ■  nxn  plane  aatrlx 
B  -  nxa  input  aatrix 
x  -  nxl  stata  vactor 
u  -  nxl  control  input  vactor 
Thusa  syataa  paraaatars  ara  of  tha  fora: 

0  ♦  I 


A  - 


N  1  *  |"v  ” 

2 

•  -2c<n 

L  NJ  •  L  \ 


B  - 


(6) 


x  ■ 


Tha  coaplata  stata,  howavar ,  is  noraally  not  availabla,  so  Eq  5  aust  ba 
supplaaantad  by  an  output  aquation.  Stata  spaca  fora  gi vas  tha  sanaor 
output  as 
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y  -  c?q  +  Cyq  (7) 

when  both  position  (p  subscript)  end  velocity  (v  subscript)  sensors  are 
used.  Expressing  this  in  ststs  vector  x: 

y  -  Cx  (8) 


where 


(9) 


Equations  S  and  8  fore  the  large  space  structure  model  available  to  the 


control  designer.  These  equations  will  be  further  explained  so  they  will 
hold  more  significance  when  being  applied  to  modal  control  of  flexible 


structures. 


Control  Model 

The  full  structural  model  is  represented  by  the  2n-d imens iona 1  state 
vector  x.  As  noted  esrller,  it  is  impossible  to  model  all  of  the  possible 
modes  for  a  complex  structure,  and  of  those  modeled,  even  fewer  will  be 


actively  controlled.  Assuming  that  multiple  controllers  are  available, 
each  controlling  a  small  subset  n^  of  nodes,  as  in  this  investigation,  the 
state  vector  nay  be  simply  represented  by 


The  x^  terms  represent  2n vectors  of  nodal  amplitudes  and  velocities  as 
defined  by  the  last  of  Eq  6  controlled  by  the  1  th  controller  of  II  con¬ 
trollers  present.  The  5  represents  a  2nr~vector  of  residual  nodes  and 

—  • 
the  the  x  represents  a  2n  -vector  of  usmodeled  modes. 


As  defined  earlier,  the  unmodeled  modes  are  those  which  exist  but  are 
beyond  the  number  of  nodes  in  the  structural  model.  These  will  no  longer 
appear  in  the  derivations.  The  residual  nodes  are  those  which  are  nodeled 
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but  not  controlled.  The  controlled  nodes  ere  those  which  require  ective 
control  in  order  to  obtain  satisfactory  system  response.  The  selection 
of  the  aodes  to  be  controlled  and  their  assignment  to  one  of  the  N  con¬ 
trollers  is  left  to  the  control  designer. 

It  should  be  noted  et  this  point  that  suppressed  nodes  were  not 
directly  referred  to  above  even  though  they  were  defined  earlier  in  the 
text.  They  have  not  been  ignored,  in  fact  they  are  included  within  the 
controlled  nodes,  in  the  following  manner:  In  a  aultlple  controller 
design,  the  individual  controller  actively  controls  those  nodes  assigned 
to  it  and  "ignores"  the  residual  nodes  defined  above.  But  the  nodes 
assigned  to  the  other  controllers  still  interact  with  this  individual 
controller,  causing  control  and/or  observation  spillover.  In  the  process 
of  controlling  the  system,  each  controller  contributes  to  the  elimination 
of  observation  and  control  spillover  in  the  system.  Thus,  each  controller, 
in  effect  "suppresses"  the  modes  contained  in  the  other  controllers.  In 
other  words,  the  controlled  modes  of  one  controller  are  the  suppressed 
modes  of  another  controller.  Therefore,  the  suppressed  modes  are  con¬ 
tained  implicitly  within  the  controlled  modes.  So,  like  the  unnode led 
nodes,  the  suppressed  nodes  exist  in  the  system,  but  will  not  be  men¬ 
tioned  any  longer  In  the  derivations  since  they  are  included  implicitly 
in  the  controlled  nodes. 

Continuing  with  the  derivation,  the  notation  of  Eq  10  nay  be  used 
to  express  the  stats  equations  as  follows 

“i  “  Vi  4  *i“  1-1,2 . ,  W  (11) 

*r  *  Ar*r  ♦  *r“  (l2) 
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J  »  I  C.x.  +  C  x 

i.i  1 1 


liter*  th*  A,  B  and  C  matrices  arc 


'  N 


~ 

-2cuj 


J  *  1|  2,  ...»  N,  r 


c.  ■  C  ,  #  * 

3  p3  *, 


c  ,  * 

vj 


3  m  2*2,  «••*  N,  r 


3  -  1.  2,  ....  H, 


T  T 

Moreover,  Che  lower  partition  of  Eq  IS,  the  •  and  •  matrices,  are  of 
the  fora 

*\  -  (17) 

♦\  -  *r  (18) 


where 


^t’jk  *  *1  ‘  dlk 


<Vjk  ‘  *1  '  drk 


The  are  the  column  vector*  of  the  matrix  ♦  and  the  d^^  and  d^  are  the 
column  vector*  of  the  and  metric**,  respectively.  Using  th*  form* 
given  in  Iqs  17  to  20,  the  and  B*Dr  matrices  may  be  represented  as 
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Ti  - 


f  • 
r 


<♦  i>ll 

^i?12  •  •  •  • 

<*i>l  n 

a 

(*  1>21 
• 

• 

(Pi>22 

<*i>2  n 

a 

e 

e 

**1^2 

(*‘V . 

^rHl 

(Vl2 

<*r>l  n 

a 

<*r>21 

a 

e 

a 

C*r^22 

<*r>2  n 

a 

(*r)n  1 
r 

(*r}nr2 

(*r}n  n 
r  a 

where  n.  is  Che  number  of  modes  in  Che  i  ch  concroller,  n  is  Che  number  of 
i  *  r 

residual  modes,  and  n&  is  che  number  of  acCuaCors  employed.  This  allows 
Eq  15  Co  be  rewriCCen  as 

0 

B, .  j  -  1,  2,  ....  N,  r  (23) 


In  simpler  Cerms,  che  rows  of  Che  pj  macrices  represent  che  amplicude  of 

each  sCrucCural  mode  along  che  line  of  acCion  of  each  acCuaCor  location. 

The  dimension  of  Che  ^  matrices  is  n^  x  making  the  dimension  of 

Che  Bj  macrices  2n^  x  n^  when  Che  upper  null  part it ion  is  Included.  Like- 

wide  lc  can  be  seen  chat  the  C  .*  and  C  .♦  partitions  of  Eq  16  are  of 

PI  vj 

dimension  n^  x  n^  where  n^  is  the  number  of  sensors  employed.  This  makes 
che  dimension  of  the  matrices  n#  x  2n^. 

Examining  the  matrix  more  closely,  the  and  terns  are  the 
position  and  velocity  coefficient  matrices,  respectively,  of  the  sensors 
employed,  aesumlng  that  both  poeitlon  and  velocity  sensors  are  used. 
However,  in  this  study,  only  position  sensors  are  used.  This  makes  che 


into  sero  matrices  so  that  Eq  16  now  beconas 
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* 


CJ  •  [CPJ  ♦ ;  °]  J  ■ 2- 2 . "• r  <24> 

Furthermore,  when  collocated  sensor /actuators  are  employed  with  the  same 

alignment,  this  simplifies  even  more.  In  this  special  case 

T 

(25) 

(26) 

(27) 


so  that 


r.,.1 

m 

i — 

Hd 

♦ 

1 _ 

pj 

J  J 

L  J 

IV. 

m 

L  J 

”  T  “ 

pj 

J 

L  J 

c. 

m 

T 
f , 

J 

j 

’] 


J  -  1,  2, 


,  N,  r 


This  simplicity  is  the  prime  advantage  of  using  position  sensors  only. 

T 

As  was  pointed  out  in  the  Bj  matrix,  the  columns  of  the  matrix  in 
Cj  represent  the  amplitudes  of  each  structural  mode  at  each  sensor  location 
along  the  line  of  the  sensor. 

The  equations  thus  derived  are  very  general  in  form  and  are  indepen¬ 
dent  of  structural  complexity.  Only  the  matrix  dimensions  will  vary 
depending  on  the  number  of  sensors,  actuators  and  modes  studied.  This 
general  development  can  lead  ore  to  understand  the  wide  variety  of 
structures  to  which  the  following  analysis  may  be  applied. 

Modal  Control 

The  controller  design  for  N  controllers  will  be  based  upon  the  model 
given  by  Eqs  11  and  13  which,  restated,  are 


«1  -  \  +  »t  « 


N 

’  *  &  Vi  +  Vr 


The  state  feedback  control  desired  is  of  the  form 
N 

5  -  i  g  ; 

i-i  1  1 


(id 

(13) 


(28) 


where  are  the  control  gain  matrices.  The  development  of  the  matrices 
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will  be  pruwttd  shortly. 

To  form  this  active  control  u,  complete  knowledge  of  the  state  vector 
Is  needed.  Unfortunately,  direct  neasurenent  of  the  entire  state  vector 


is  1 ■possible  and  the  only  neasure  of  x  Is  the  neasurenent  y  given  by 
the  sensors.  As  a  result,  it  Is  necessary  to  develop  a  state  est tutor 
which  can  taka  the  observations  y  and  produce  an  estimate  of  x.  This  can 
be  done  by  employing  an  observer  of  the  fora 


*1  *  Ai*l  +  Blu  +  Ki(y  “  yi> 


7±  “  (30) 

where  x^  are  the  estimated  states,  are  the  estimated  outputs  and  are 
the  observer  gain  matrices.  The  matrices  are  chosen  such  that  the  error 

la  the  state  estimate 

\  (31) 

is  asymptotically  stable,  how,  the  control  vector,  in  terns  of  the  esti¬ 


mated  state  Is  given  by 

*  -  ?  G.  i.  (32) 

1-1  1  1 

Iquations  11,  13,  29,  30  and  32  represent  the  control  problem  for  a  large 


apace  structure. 

Before  proceeding  any  further,  time  will  be  taken  now  to  develop  the 
control  gala  matrices  6^  and  observer  gain  matrices  K^.  For  this,  linear 
optimal  regulator  theory  (Bef  9)  is  used.  The  control  gain  matrix  G  is 

e 

derived  first,  starting  by  defining  a  quadratic  performance  index  J 


h  that 


J  •  *1  f  (**  1  Q4««  ♦  ■  T»J«) 


(33) 


where 


Q  is  an  n  x  n  positive  semideflnite  weighting  matrix 
t  is  n  i  x  i  positive  definite  weighting  matrix 


It  is  desired  to  minimize  this  index  subject  to  Eq  11.  Then,  the  optimal 
solution  to  the  minimization  problem  is 

Gt  -  -R^1  St  (34) 

where  is  the  solution  to  the  steady  state  matrix  Riccati  equation: 

T  **  1 

SiAl  +  Ai  Si  "  SiBiRi  ®iT  Si  +  Qi  *  0  (35) 

Realizing  that  the  eigenvalues  of  the  matrix  (A^-K^C^)  are  the  same  as 

j  XX 

the  eigenvalues  of  its  transpose  -C^K^  ),  a  similar  development  may 
be  used  for  the  observer  gain  matrices  An  equation  for  the  system 

similar  to  Eq  11  can  be  written  using  the  state  w: 


*  V 


'i  *i 


where  the  g  is  the  control  input  given  by 


«i  ■  riT  “i 


Again,  using  linear  optimal  regulator  theory,  a  quadratic  performance  index 


is  defined: 


m 

J  '  H  f  <;1T  Qobl  ;i  +  *1  Robi  *1> 


where  Q  and  R  .  are  weighting  matrices  as  defined  for  Eq  33,  but  are  not 


necessarily  the  sanm  exact  matrices  for  the  6^  matrices.  It  may  be 
desirable  to  weight  the  observation  data  more  or  less  than  the  control 

feedback. 

Continuing,  Eq  38  is  minimised  subject  to  Eq  36  and  the  observation 


gala  matrix  Is  given  by 


K.  -  +R  ”1  C.  P. 
i  obi  i  i 


(39) 


where  p  is  the  solution  to  the  steady  state  matrix  Riccati  equation: 


piait  +*ipi  -piciTRoi,r1  cipi +  w°  <4° 

It  nay  be  noted  that  Eq  40  is,  in  effect,  the  transpose  of  Eq  35,  which 
follows  since  the  observation  gain  matrix  is  developed  from  the  transpose 
of  the  system  matrix  (A^  -  K^C^) . 


Equations  34  and  39  now  form  the  control  gains  and  observation 
estimator  gaii.s  to  be  used  in  their  respective  controllers,  and  are 
determined  such  that  each  controller  is  stable.  However,  due  to  coupling 
these  controllers  and  observers  are  not  independent.  Therefore,  even 
though  the  and  matrices  keep  their  individual  controllers  stable, 
the  overall  system  may  be  unstable. 

This  system  instability,  even  with  stable  controllers,  can  be  seen  in 
the  following  development.  To  begin,  an  N-controller  system  will  be 
illustrated.  Then  a  three-  and  four-controller  system  will  be  shown.  It 
will  become  obvious  that  controller  stability  alone  will  not  guarantee  over¬ 
all  system  stability.  The  one-  and  two-controller  systems  were  demonstrated 
by  Miller,  and  he  proceeded  as  far  as  deriving  the  three  controller  case. 

N  Controller 

For  a  multiple  input  multiple  out  controller,  the  state  equations  are 
given  by  Eq  11  as 

x^  *  A^x^  *  ®iu  i  *  1,  2,  ...,  N  (41) 

but  for  reasons  stated  earlier,  an  observer  of  the  form  given  by  Eqs  29 
and  30  is  used 

Xi  “  Aixi  +  Bi“  +  Ki^  “  yi^  i  -  1,  2,  ...,  N  (42 ) ^ 
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■  C^x^  i  ■  1,  2,  ....  N  (43) 

where  cho  subscript  1  denotes  the  1  th  controller,  the  estimated  state 
vectors  and  y  the  estimated  output  vectors. 

The  observation  gain  matrices  are  chosen  such  that  the  state 
estimate  errors 


ei  -  xi  -  xt 


i  -  1,  2,  N 


(44) 


approach  zero  for  large  time.  The  control  is  then  given  by 
N 

u  -  1  G  x 

i-1  1  1 


(45) 


Equations  42,  43,  and  44  may  then  be  combined  with  the  state  equation 
given  in  Eq  41  to  obtain  the  state  estimate  errors: 

•  a  • 

l±  -  x±  -  x±  =  (A±-  KiCi)ei  +  l  KjC.-X,  i  -  1,  2,  ...,  N  (46) 

3*1  '  3  ’  2 . "•  r 

Now,  using  the  state  equations  given  in  Eq  41,  along  with  the  control 

stated  in  Eq  45  above,  the  controlled  state  equations  are  given  by 

N 

*  (A±  +  BiG1)xi  +  BiGiei  +  I  B^x^  i  -  1,  2,  ....  N  (47) 


The  states  x^  and  errors  in  the  states  e^  may  be  collectively 
evaluated  by  the  controlled  system  state  presented  as  an  augmented  state 
vector  of  Eqs  46  and  47.  This  vector  z  is  given  by 


f-T-T-T-T  -T-T-t") 

(X1  »  el  »  x2  ’  e2  . ’  *N  »  eN  ’  xr  ) 


(48) 
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Writing  out  the  closed  loop  state  equations,  in  terms  of  z,  gives 


A^  B,G 


1"1 


B1G1 


—  BlGi 


A^K.C 


ri 


—  Kici 


BiGi 


B 


KiCl 


i  1 
0 


B  G, 
n  1 


B  G. 
n  1 


B  G. 
n  i 


K  C. 
n  1 


K  G 
n  i 


B  G. 
r  1 


BrGl 


B  G, 
r  i 


BlGi 

—  BlGn 

BlGn  0 

0 

“  Vn 

0  K1Cr 

BiGi 

-  B.G 

1  n 

KiCr  0 

Ai-Kici 

K.C 
i  n 

0  0 

BnGi 

-  A  +B  G  B  G  0 

n  n  n  n  n 

0 

0 

A  -K  C  K  C 
n  n  n  n  i 

B  G. 
r  i 

-  B  G 

r  n 

B  G  A 

r  n  r 

The  system  matrix  may  be  either  upper  or  lover  block  triangularlzed. 
For  upper  block  triangular  form,  the  spillover  elimination  required  is 


BiGj  *  o 


and 


j  -  1,  2,  ....  N-l;  i  -  j  +  1 . N 
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KiCj 


and  for  lower  block  triangular  form 


BiGj  “  ° 


and 


i  -  1,  2,  ....  N-l;  j  -  i  +  1,  ....  N 


(51) 


KiCj  *  ° 


The  block  structure  of  Eq  49  can  be  more  easily  seen  in  specific  exasr- 
pie,  therefore  the  three  and  four  controller  cases  will  now  be  examined. 
Three  Controllers 


Setting  N  ■  3  and  following  the  form  given  in  the  previous  development 


for  N  controllers,  the  state  equations  for  a  three  controller  system 
are  given  as 


x,  ■  +  BjU 


x2  -  A2x2  +  B2u 


X3  “  A3*3  +  B3u 


where  the  subscripts  designate  the  controller  described  by  the  equation. 
The  observer  to  be  used  is 

*i  ■  Vi  +  v  +  Ki(y  -  V  1  -  1,2,  3  (5! 


yt  -  i  -  1,  2,  3  (56) 

Again,  the  observer  gain  matrices  are  chosen  such  that  the  error  in  the 
state  estimates 


ei  *  *i  “  xi 


i  -  1,  2,  3 


are  asymptotically  stable.  And  now  the  control  applied  is  given  by 


U  "  Glxl  +  G2x2  +  G3x3 

Using  Eqs  52  through  57,  the  system's  state  estimate  errors  may  be 
shown  as  described  in  Eq  46  for  N  ■  3: 


*2  "  x2  “  x2  "  +  K2C3X3  +  *2Crxr 

•  •  • 

e3  ■  Xj  -  Xj  ■  (AyKjCjJej  +  K^C^x 
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upper  or  lower  block  trlangularization  nay  be  attempted,  la  either  case 
the  resulting  eigenvalues  will  be  the  same  since  they  are  determined  by  the 
block  diagonal  terms.  The  observation  spillover  and  control  spillover  terms 
to  be  eliminated  for  both  schatws  arc  presented  in  Table  VXII. 

Table  VIII 

Spillover  Elimination  for  Decoupling  Three  Controllers 
Upper  Trlangularization  or  Lover  Tr iangular last ion 


B2G1  *  0 

B1G2  “  0 

B3G1  ‘  ° 

B1G3  ’  ° 

B3G2  “  ° 

B2G3  "  ° 

«aci-° 

*1C2  “  0 

K3C1  *  ° 

J* 

uT 

■ 

o 

*3C2-0 

¥3"° 

There  are  actually  two  methods  of  approaching  the  elimination  of  spill¬ 
over  in  the  three  controller  system.  The  first  of  these  deals  with  selective 
placement  of  tha  sensors.  If  the  sensors  are  positioned  in  a  proper  manner, 
the  model  amplitude  matrix  •,  and  as  a  result,  the  system's  B  and  C  matrices, 
will  be  of  the  form  that  a  selective  arrangement  of  controlled  modes  will 
make  one  controller orthogonal  to  the  other  two  controllers.  In  effect,  this 
reduces  the  problem  to  a  two  controller  aystem  since  the  two  controllers  will 
have  no  effect  on  a  controller  orthogonal  to  them.  Then,  only  the  spillover 
terms  between  the  two  eon-orthogonal  controllers  need  be  considered.  This 
was  specifically  demonstrated  by  Miller  in  his  examination  of  a  three  con¬ 
troller  system. 

To  present  an  example,  aastaM  the  controller  three  modes  are  ortho¬ 
gonal  to  the  controller  one  and  two  modes.  Then  all  cross  terms  between 
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on*  and  thr**  and  two  and  Chrae  vill  be  aero,  which  effectively  reduces 
Table  VIII  to  the  spillover  elimination  requirements: 
for  upper  triangulat  form 

B2Gx  -  0  and  K2C1  -  0  (68) 

and  for  lower  triangular  form 

BjG2  -  0  and  KjC2  -  0  (69) 

It  can  be  shown  that  these  terms  are  the  seme  as  those  required  to  block 
triangularlse  the  two  controller  system,  thus  demonstrating  the  three  controller 
reduction  to  a  two  controller  system.  However,  since  the  intent  of  this  study 
la  to  examine  a  three  and  four  controller  system,  this  approach  will  not  be 
addressed  further. 

The  second  approach  proposes  that  in  the  transformation  a  matrix  (or 
Tj)  be  found  which  will  drive  both  and  (K^  and  K^C^)  to  zero,  and 

a  matrix  T^  (or  T^)  be  found  which  will  drive  both  B^  and  B^G^  (B^G,  and 
B2G3)  to  zero.  This  process  will  be  explained  in  the  next  chap^r,  however, 
for  now,  it  is  not  always  possible  to  find  transformation  matrices  which 
will  reduce  the  spillover  terms  to  zero.  On  the  other  hand,  it  is  also  not 
always  possible  to  position  the  sensors  such  that  there  are  two  orthogonal 
sets  of  modes  to  reduce  the  system  to  two  controllers.  Therefore,  a 
compromise  is  in  order:  If  two  clearly  orthogonal  sets  of  modes  are  not 
available,  the  modes  should  be  assigned  so  they  are  as  nearly  orthogonal 
as  possible.  Then  the  transformation  is  performed  on  the  full  set  of  spill¬ 
over  terse  given  in  Table  VIII.  This  investigation  will  assume  the  modes 
are  not  fully  orthogonal,  but  still  assign  the  modes  to  make  the  controllers 
as  orthogonal  as  possible,  and  attempt  to  find  transformation  matrices 
which  will  make  each  spillover  term  in  Table  VIII  approximately  zero. 
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Four  Controllers 


<L 


Th«  four  controller  development  is  identical  to  that  of  three  con¬ 
troller!,  so  the  basic  results  only  will  be  presented.  The  state  equations 
arc  given  by  setting  M  *  4  in  Eq  41,  and  the  observers  are  given  by  Eqs  42 
and  43,  again  setting  N  ■  4.  The  state  estlaete  error  is  given  by  Eq  44 
and  the  resulting  control  applied  is 

“  "  Glxl  +  G2x2  *  C3X3  +  G4x4 

To  avoid  repeating  the  equations  already  presented,  suffice  it  to  say 
that  the  state  estimate  errors  are  of  the  fora,  given  for  the  first  con¬ 
troller 

• 

”  *1  "  (A^-K^C^)e^  +  k^C2x2  +  ^1G3X3  *S.G4X4  +  ^lGrXr 

and  the  full  state  equations  are  of  the  fora,  again  given  for  the  first 
controller 

*1  -  (Aj+BjG^x^  +  BjGji^  +  BjG2*2  ♦  BjGjij  +  BjG^  (71) 

The  reaaining  three  controllers'  errors  and  states  aay  be  derived  froa 
Eqs  46  and  47,  respectively. 

These  aay  be  combined  into  the  controlled  system  state  vector  z, 
defined  as 


aad  the  closed  loop  system  equations  aay  be  given  as 
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•IN 


V»1G1 

B1G1 

B1G2 

B1G2 

B1G3 

B1G3 

B1G4 

B1G4 

0 

0 

Wi 

KlC2 

0 

Kl°3 

0 

K1C4 

0 

KlCr 

Vi 

B2G1 

A2+b2C2 

B2G2 

B2C3 

B2G3 

B2G4 

B2G4 

0 

*2C1 

0 

0 

V  K2C2 

*2C3 

0 

*2C4 

0 

Vr 

B3G1 

B3°l 

B3G2 

B3G2 

A3+B3g3 

B3G3 

B3G4 

B3G4 

0 

K3C1 

0 

K3C2 

0 

0 

A3~K3C3 

K3C4 

0 

K3Cr 

Vl 

B4G1 

B4G2 

B4G2 

V’3 

B4G3 

A4+B4G4 

B4G4 

0 

Vl 

0 

K4C2 

0 

V'3 

0 

0 

A4“K4C4 

K.C 

4  r 

BrGl 

BrGl 

BrG2 

BrG2 

BrG3 

BrG3 

B  G. 
r  4 

BrG4 

A 

r 

The  spillover 

terms  to 

be  eliminated  are 

given  in 

Table  IX, 

and 

correspond  Co  Che  condicions  given  in  Eqs  SO  and  SI.  IC  is  apparenc  Chls 


(73) 


syscem  is  more  difficult  Co  block  criangularize  than  the  three  controller 
system.  However,  a  prime  advantage  to  using  additional  controllers  is  the 
number  of  modes  to  be  controlled  may  be  divided  among  more  controllers. 

This  reduces  the  order  of  each  controller  thus  reducing  the  burden  on  the 
computer.  This  is  especially  important  in  solving  the  matrix  Riccati  equa¬ 
tions,  Eqs  35  and  40,  since  the  computational  burden  is  approximately  the 
cube  of  the  order  of  the  equation  (Ref  3). 
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Table  IX 

Spillover  Elimination  for  Decoupling  Four  Controllers 
Upper  Triangulariaation  or  Lover  Trlangularization 


Vi  - 0 

K2C1 

-  0 

B1G2 

-  0 

*S.°2  ~  0 

B3G1  “  0 

Vl 

-  0 

B1G3 

-  0 

hc3m0 

b3g2  -  0 

K3C2 

-  0 

B1G4 

-  0 

KiC4-° 

B4G1  *  0 

*4C1 

-  0 

B2G3 

-  0 

^2C3  *  0 

Va  - 0 

K4C2 

-  0 

B2G4 

-  0 

¥4‘° 

B4G3  "  ° 

*4C3 

-  0 

B3G4 

-  0 

^4-0 

Sensor /Actuator  Requirements 

As  mentioned  for  three  controllers,  in  order  to  perform  spillover 
suppression,  one  or  more  gain  matrices  must  be  made  orthogonal  to  N-l 
B  or  C  matrices.  For  example,  from  Eq  50,  to  satisfy  the  expression 
for  BjGj  the  columns  of  must  be  simultaneously  orthogonal  to  the  rows 
of  through  B^.  In  order  words,  the  columns  of  must  be  in  the 
null  space  of  the  matrix  B2N  where  is  defined  as 


B2N  " 


The  null  space  of  B,M  has  dimension  P-„  given  as 


f2*  ’  -  r2K> 


where 


n  »  number  of  actuators 
a 

r2N  -  rank  of  B2N  -  min  (n2  +  n3  +  . . .  +  nN,  nft) 

Therefore,  G^  has  P2N  columns. 

The  number  of  actuators  must  exceed  the  rank  of  B2n  in  order  for 
*2N  Gi  *  °*  Otherwise  the  system  is  overspecified  and  no  transformation 
matrix  exists  which  will  drive  the  B^  G^  to  zero.  If  the  rows  of  B2N 
are  linearly  independent ,  then  the  number  of  actuators  needed  is  given  as 


n  > 
a 


N 

£ 

i-2 


(76) 


and  if  the  rows  are  not  linearly  independent,  aa>T2n‘  can  be  seen 
that  the  other  control  gain  matrices  will  have  a  sufficient  number  of 
actuators  is  the  inequality  in  Eq  76  is  met*.  A  similar  study  shows  that 
the  number  of  sensors  needed  is,  for  C 2JJ  with  linearly  Independent  columns 


ns  * 


N-l 

£ 

1-1 


(77) 


a.'**  for  roluens  that  are  not  linearly  independent,  n#>  r2N  where  r 


2N 


here  is  the  rank  of  C2R. 


Likewise,  for  the  lower  block  triangular  conditions  given  in  Eq  51, 
the  actuator  and  sensor  requirements  can  be  shown  to  be 


n.  * 


n.  > 


N-l 

£ 

i-1 

N 

£ 

i-2 


(78) 


(79) 


for  linearly  independent  rows  and  columns  of  the  B2g  and  C2N  matrices, 
respectively. 

It  should  be  pointed  out  that  satisfying  the  Inequalities  given  In 
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Eqs  76  through  79  may  actually  require  more  actuators  and  sensors  than 
indicated.  As  an  example,  consider  a  thirty  node  model.  Using  three 
controllers,  each  with  ten  modes,  the  above  inequalities  require  at  least 
twenty  sensors  and  twenty  actuators  to  assure  decoupled  system  stability. 
However,  to  control  and  observe  the  system,  at  least  one  more  sensor  and 
actuator  would  be  required.  These  conditions  must  be  set  in  order  to 
generate  the  transformation  matrices  and  to  implement  the  controllers. 
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IV.  Transformation  Technique 

It  has  bean  mentioned  several  times  that  the  close  loop  state  equa¬ 
tions,  Eqs  49,  67  and  73,  will  be  put  into  block  triangular  form  by 
the  selective  elimination  of  control  and  observation  spillover  tens. 
However,  the  exact  details  of  this  spillover  elimination  have  been 
neglected  until  now.  The  following  will  describe  the  generation  of  the 
transformation  matrices,  which  were  referred  to  specifically  as  T  and  T 
in  the  previous  section.  The  T  matrix  is  a  transformation  matrix  for  the 
control  spillover  and  the  T  matrix  is  a  transformation  matrix  for  the 
observation  spillover. 

In  a  single  controller  case,  it  can  be  seen  that  the  spillover  terms 

which,  when  eliminated,  will  assure  system  stability  are  B  G  or 

8 

KCg,  where  the  s  subscript  designates  modes  to  be  suppressed.  An  issMd- 

lately  obvious  solution  to  this  is  G  *»  0  or  K  *  0.  But  this  solution  will 

make  the  respective  controllability  term  BcG,  or  observability  term 

KC  ,  equal  to  zero  also.  Therefore,  this  solution  is  unsatisfactory.  'x>\e 
c 

transformation  method  generates  a  solution  which  for  a  single  controller, 
is  subject  to  the  conditions: 


BgG  -  0  (80) 

or 

KC  -  0  (81) 

s 

while  maintaining 

B  G  *  0  (82) 

c 

or 

KC  d  0  (83) 

c 


It  would  also  be  desirable  to  apply  the  additional  constraint  to  the 
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residual  sodas 


B  G  -  0 
r 

(84) 

Kr  -  0 

(85) 

however ,  due  to  the  large  number  of  structural  nodes  present  In  the  model, 
this  constraint  is  not  realistic  and  will  be  ignored  in  this  development. 
The  effects  of  the  residual  spillovers  may  be  minimized  by  the  careful 
selection  of  modes  designated  as  residual  or  suppressed,  so  as  to  create 
a  frequency  separation  between  the  residuals  and  the  bandwidth  of  the 
controller. 

For  a  multiple  controller,  the  conditions  given  in  Eqs  80  and  81 
apply,  but  now  the  B  and  C  matrices  may  take  on  the  form  illustrated 
by  Eq  74  and  will  be  referred  to  as  the  and  matrices.  Instead  of 

discussing  all  of  the  possible  combinations  of  B^  and  C^N  for  N  con¬ 
trollers,  take  as  an  example  the  first  condition  given  in  Eq  50,  that  is 

%±Gj  -  0  j  -  1,  2,  ...,  N-l;  i  -  j  +  1,  ....  N  (86) 

Given  H  controllers,  will  have  to  be  made  orthogonal  to  N-l  matrices. 
The  N-l  B?  matrices  may  be  combined  into  a  single  matrix  such  that 


[±] 


Therefore,  one  of  the  conditions  to  be  mat  is  B^  •  0.  In  other  words, 
tbs  matrix  myst  be  transformed  such  that  its  colons  are  orthogonal  to 
tbs  rows  of  Bjji*  or^  as  stated  in  the  previous  section,  G^  must  be  in  the 
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null  space  of  B2JJ -  This  is  the  most  difficult  case  for  N  controllers  and 
is  presented  only  to  describe  how  the  multiple  matrix  is  set  up.  The 
remainder  of  the  derivation  will  be  in  terms  of  a  generic  Bg  matrix  which 
represents  and  B2  alike. 

The  transformation  matrix  sought  will  be  referred  to  as  T  and  will  be 
such  that 

B  T  -  0  (88) 

s 

B  has  the  row  dimension  of  n  (the  number  of  modes  to  be  suppressed)  and 
3  IQ 

the  column  dimension  of  n  (the  number  of  actuators).  T,  therefore,  has 

a 

dimensions  of  n  by  n  -  n  .  If  there  are  fewer  actuators  than  linearly 
a  a  m 

independent  modes,  then  no  solution  matrix  T  exists.  The  system  is  over- 
specified  in  twis  case,  meaning  there  are  more  equations  than  unknowns. 

If  the  number  of  modes  and  actuators  are  equivalent,  then  the  system  is 
stable  but  uncontrollable  (or  in  the  case  of  the  KC  matrices,  the  system  is 
stable  but  unobservable  assuming  an  equal  number  of  sensors  and  actuators). 
The  actuators  (sensors)  are  saturated  with  maintaining  stability  alone. 

Simply  stated,  the  conditions  given  in  Eq  76  and  77  must  be  met  in 
order  to  generate  a  transformation  matrix. 

Tj  Illustrate  the  result  of  applying  the  transformation  matrix  described 
above,  consider  the  following  system  of  controlled  and  suppressed  modes: 

x  -Ax  +  B  u  (89) 

c  c  c  c 

x  «  A  x  +  B  u  (90) 

s  s  s  s 

where 

u  -  G  x  (9D 

c  c 

The  Bsu  term  of  Eq  90  is  a  control  spillover  term  which  may  be  adversely 
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affected  by  the  control  applied  to  Eq  89.  The  elimination  of  this  term 

requires  the  use  of  a  transformation  matrix  T  such  that 

B  T  *  0  (92) 

s 

while  maintaining 

BcT  4  0  (93) 

This  transformation  matrix  may  be  used  to  define  a  new  control  v  as 

u  -  Tv  (94) 

Inserting  this  expression  into  the  state  equations  given  in  Eqs  89 
and  90  yields 


x  -  A  x  +  B  Tv  (95) 

x  *  A  x  +  B  Tv  (96) 

s  s  s  s 

Letting  B  T  «  B*  and  knowing  that  B  T  a  0  the  new  system  is  described  by 
c  s 

« 

x  -Ax  +  B*v  (97) 

c  c  c  c  v  ' 


x  “Ax  (98) 

s  s  s 

in  which  no  controller  spillover  exists.  The  new  control  vector  will  be 
shown  to  be 

v  -  G*x  (99) 

c  c 

With  this  general  overiew  of  the  tranformation  process  and  its  results, 

the  development  of  the  matrix  T  will  now  be  discussed. 

The  major  tool  used  obtaining  this  result  is  called  the  Singular 

Value  Decomposition  (Ref  4).  The  matrix  to  be  decomposed  is  B  which  has 

s 

dimensions  n  x  n  and  can  be  described  by 

IQ  A  J 


B#  -  WIVT 


(100) 
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and 


W  is  an  (n  x  n  )  orthogonal  matrix  of  left  singular  vectors 
mm 

V  is  an  (n  x  n  )  orthogonal  matrix  of  right  singular  vectors 

&  3l 


E  is  an  (n  x  n  )  matrix  with  the  s  singular  values  of  B 
m  a  3 

in  the  first  s  entries  along  the  main  diagonal  and  zeroes 
in  all  other  positions : 


(101) 


n  x  n 
m  a 


such  that 


S  - 


(102) 


The  total  number  of  singular  values  present  is  equal  to  the  rank  of  the 

B  matrix,  and  they  are  all  non-negative.  Assuming  B  is  of  full  rank,  then 
s  s 

s  ■  min  (n  ,  n  ). 

a  m 

By  partitioning,  the  W  matrix  can  be  defined  by 


where 


(103) 


W_  is  a.i  n  x  s  matrix  of  left  singular  vectors  associated 

8  ID 


with  the  non-zero  singular  values 


W  is  an  n  x  r  matrix  of  left  singular  vectors  associated 
r  m 

with  the  zero  singular  values 

and 

s  +  r  *  n  (104) 

m 

Similarly,  by  partitioning,  the  V  matrix  can  be  defined  by 

V  -  V  V  (105) 

s  p 

where 

V  is  an  s  x  n  matrix  of  right  singular  vectors  associated 
S 

with  the  non-zero  singular  values 

V  is  a  p  x  n  matrix  of  right  singular  vectors  associated 

P  a 

with  the  zero  singular  values. 

and 

s  +  p  *  n  (106) 

d 

Remembering  the  V  matrix  is  orthogonal  and  noting  then  that 

V  T  V  -  0  (107) 

s  p 

the  decomposed  matrix  may  be  written  as 

BT-WSVTV-BV  -0  (108) 

S  8  S  p  S  p 

which  leads  to  the  conclusion  that  the  transformation  matrix  desired  is 
given  in  the  matrix  of  right  singular  values  associated  with  the  zero 
singular  values: 

T  -  V  (109) 

P 

where  T  i  0. 

Once  the  transformation  matrix  is  found,  implementation  is  relatively 
simple.  Equation  94  defined  a  as  the  new  control  input  vector.  This 
is  now  given  as 
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The  G*  matrix  Is  found  in  the  same  manner  as  in  the  Modal  Control  section 
of  Chapter  III  In  which 

G*  -  -R*"1  B*T  Si  (111) 

where 

R*  T 

i  -  Ti  (112) 

R^  is  the  positive  definite  weighting  matrix  in  Eq  33 

B*  -  (113) 

and  is  the  solution  to  the  matrix  Riccati  equation: 

Si  Ai  +  AiT  Si  "  Si  Bi  R*-1  B*T  Si  +  Q±  "  0  (114) 

Simple  manipulation  of  Eq  111  will  show  that  the  transformed  gain  matrix  is 
finally  given  by 

G*  -  Tt  Gt  (115) 

Substitution  of  this  back  into  the  state  equations  yields  a  closed  loop 

system  which  is  block  triangular  with  no  control  spillover. 

This  technique  may  be  paralleled  to  obtain  a  r  transformation  matrix 

for  the  observation  gain  to  eliminate  observation  spillover.  Substituting 
T  T 

C  for  B  ,  K  for  G  and  T  for  T  will  give  the  same  results  with  T  equal 
8  8 

to  Vp.  Ohtervise,  a  new  derivation  using  Cg  and  K  will  give  the  result 
that  the  observation  gain  transformation  matrix  T  is  equal  to  the 
transpose  of  the  matrix  of  left  singular  vectors  associated  with  the 
zero  singular  values  of  C&.  Here,  as  alluded  to  earlier,  the  number  of 
sensors  must  be  greater  than  the  number  of  modes  to  be  suppressed. 


V. 


Simplicity  and  flexibility  were  prime  consider at lone  in  the  develop¬ 
ment  of  the  computer  programs.  The  general  format  given  by  Miller  was 
followed  aa  it  provided  a  very  simple  progression  of  logic. 

Two  programs,  using  the  same  logic,  were  created.  The  same  control 
techniques  are  applied,  however ,  the  spillover  elimination  schemes  used 
are  different.  One  provides  decoupled  control  by  reducing  the  closed 
loop  state  equation  matrix,  given  in  form  by  Eq  49,  to  an  upper  block 
triangular  form.  The  other  reduces  the  matrix  to  a  lower  triangular 
form.  It  was  deemed  unnecessary  to  combine  the  two  into  one  program  with 
an  input  selecting  one  or  the  other  since  the  eigenvalue  results  provided 
by  both  programs  are  Identical.  In  fact,  by  renumbering  the  controllers 
in  the  spillover  elimination  portion  of  one  program,  the  opposite  trans¬ 
formation  is  achieved.  For  example,  in  the  upper  trinagular  three  controller 
system,  the  spillover  terms  to  be  suppressed  are  given  in  the  first  column 
of  Table  VIII.  By  inverting  the  first  column  subscripts  such  that  1,  2, 

3  becomes  3,  2,  1,  it  can  be  seen  that  the  resulting  spillover  terms  are 
identical  to  those  in  the  second  column  for  a  lower  triangular  transforma¬ 
tion.  The  subscript  inversion  is  equivalent  to  renumbering  the  controllers 
in  the  transformation  portions  of  the  program.  For  this  reason,  only  the 
program  for  the  upper  block  triangular  controlled  system  is  presented. 

This  is  listed  in  Appendix  C.  The  subroutines  which  support  this  program 
are  given  in  Appendix  D.  Several  other  subroutines  are  called  but  not 
listed.  These  are  provided  by  the  International  Mathematical  and 
Statistical  Library  (IMSL). 

Since  program  flexibility  is  desired,  race  the  nodal  data  is  loaded, 
the  program  is  designed  to  make  any  number  of  runs  with  different  para- 
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asters  for  each  run.  The  parameters  that  aay  be  varied  by  the  operator 
include:  using  a  three  or  a  four  controller  ay a tea,  which  nodes  are 
assigned  to  each  controller,  what  control  and  observer  freighting  values 
to  assign  to  each  node,  and  what  initial  systea  damping  ratio  is  applied. 

The  prograa  aay  read  the  input  data  froa  initialization  asaignawnts 
within  the  prograa  or  froa  a  permanent  data  file.  In  either  case,  the 
prograa  operates  as  if  it  were  interactive  by  prompting  for  input  and 
then  echoing  the  data  read  in.  TVs  makes  the  output  very  easy  to  interpret 
by  allowing  the  user  to  trace  the  computer's  progress  through  the  execu¬ 
tion  of  the  prograa. 

The  prograa  is  initialized  by  inputting  the  nuaber  of  controllers 
desired  and  then  the  nuaber  of  modes  in  each  controller.  If  three  con¬ 
trollers  are  used,  the  nuaber  of  residual  modes  will  be  requested,  other¬ 
wise  the  fourth  controller  systea  is  run  without  residuals.  Next  the 
nuaber  of  actuators  and  sensors  are  input,  along  with  the  modal  damping 
ratio  c.  For  the  models  studied,  the  sensors,  actuators  and  open  loop 
damping  applied  were  obtained  froa  previous  studies.  CSDL  1  was  tested 
with  six  actuators,  six  sensors,  and  a  daaplng  ratio  of  0.005  (Ref  7), 
and  CSDL  2  was  tested  with  twenty-one  actuators,  twenty-one  sensors, 

and  a  damping  ratio  of  0.01  (Ref  10).  The  program  will  then  read  froa 

T 

a  permanent  file  the  aatrix  (♦  D)  of  modal  amplitudes  at  each  actuator 
location,  followed  by  the  transpose  of  the  aatrix  of  modal  amplitudes  at 
each  sensor  location.  In  this  study,  these  two  matrices  are  identical 
since  colocated  pairs  of  actuators  and  sensors  are  eaployed.  The  sensor 
modal  amplitude  aatrix  is  input  in  transposed  fora  so  that  the  aatrix 
for  a  colocated  systea  aay  be  copied  directly  froa  the  actuator  modal 
aaplitude  r-.irix.  However,  these  are  left  as  separate  entries  in  the 


event  the  actuators  and  sensors  are  not  colocated.  Finally,  the  nodal 
frequencies  are  read  in  from  the  permanent  file.  After  this  preload  of 
data,  the  desired  run  is  made  by  specifying  which  modes  are  to  be  con¬ 
trolled  by  each  controller  and  which  modes  are  to  be  left  as  residuals, 
along  with  the  desired  control  and  observer  weighting  values  for  each 
mode.  Additional  runs  may  be  made  in  the  same  job  by  specifying 
different  modal  arrangements,  weighting  assignments  or  controller 
conf igura  t ions . 

Program  execution  actually  begins  with  the  formation  of  the  A,  B, 

C  and  weighting  matrices  for  each  controller.  This  is  conveniently  done 
by  subroutines  which  read  the  required  data  for  the  modes  specified.  These 
subroutines  allow  the  operator  to  change  the  size  of  the  controllers  as 
needed  slsn>ly  by  specifying  the  number  of  modes  to  be  placed  in  each 
controller. 

Once  the  initial  matrices  are  formed,  the  control  and  observation  feed¬ 
back  gain  matrices,  G^  and  K^  respectively,  are  determined  using  a  series 
of  subroutines  which  generate  a  numerical  solution  to  the  matrix  Riccatl 
equation.  These  sophisticated  routines  were  created  by  Kleinman  (Ref  11) 
and  so  are  known  as  the  Kleinman  routines.  The  and  K^  matrices,  along 
with  the  parameter  matrices.  A,  B,  and  C,  arc  then  combined  to  form  the 
closed  loop  system  matrix,  as  given  in  form  by  Eq  49.  This  particular 
program  develops  the  three  controller  system  with  residuals  as  in  Eq  67, 
but  unlike  Eq  73,  does  not  Include  residual  terms  in  the  four  controller 
systems. 

The  eigenvalue  analysis  of  the  system  is  performed  next,  making  use  of 


the  XSML  routine  EIGRF  which  determines  the  eigenvalues  of  real  non-sym 
trie  matrices.  First,  the  eigenvalues  of  the  overall  system  matrix  are 


c 


generated  to  show  the  stability  of  the  full  system.  Then  the  eigenvalues 
of  the  A  +  BC  and  A  -  KC  matrices  for  each  controller  are  found.  These 
values  demonstrate  the  effectiveness  of  the  Individual  controller  and 
show  which  modes  were  affected  most.  Last,  for  the  residual  modes, 
if  any,  the  eigenvalues  of  the  A  residual  matrix  are  given. 

Spillover  elimination  is  the  next  step  In  the  control  algorithm.  This, 
too,  varies  with  controller  configuration,  but  also  varies  with  the  form  of 
triangularlzation  selected.  The  3  or  4  controller  variation  is  accomoda¬ 
ted  by  the  program,  bit  the  type  of  triangularlzation  selected  determines 
which  of  the  two  programs  is  used.  This  selection  is  more  a  matter  of 
operator  preference  than  system  requirement,  though,  since  the  upper  and 
lower  block  triangular  reductions  yield  results  identical  to  several  deci¬ 
mal  places.  Regardless  of  the  program  selected,  the  spillover  elimination 

is  a  rather  lengthy  portion  of  Che  program.  The  modes  to  be  suppressed  are 

T 

formed  into  non-zero  B  and  C  matrices  of  the  form  given  in  Eq  87.  The 

S  8 

IMSL  routine  LSVDF  is  used  to  perform  a  singular  value  decomposition  on  these 
matrices.  By  using  the  left  singular  vectors  associated  with  the  zero 

T 

singular  values  of  B  and  C  ,  the  transformation  matrices  T.  and  r.  are 

3  9  XX 

formed.  The  program  then  loops  back  with  the  transformation  matrices  and 
these  are  applied  as  discussed  in  Chapter  IV  to  create  new  gain  matrices, 

G*  and  K*.  A  new  closed  loop  system,  which  has  the  selected  block  trian¬ 
gular  form,  is  generated  and  the  eigenvalue  analysis  is  repeated.  In  this 
analysis,  all  of  the  individual  controller  eigenvalues  should  be  matched 
by  identical  elgenvaluea  in  the  overall  closed  loop  system,  whereas,  in  the 
first  analysis  there  may  be  minimal  correlation,  depending  on  system 
coupling. 
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The  difference  in  sizes  of  the  two  nodels  examined  using  these  pro¬ 
grams  demonstrate,  to  some  extent,  the  high  degree  of  flexibility  of  the 
control  method  applied.  The  only  change  that  needs  to  be  made  to  adapt 
the  method  for  another  structure  is  the  basic  matrix  dimensioning  in  the 

program.  Nothing  else  has  to  be  altered,  as  long  as  the  system  model 
* 

can  be  defined  by 


and 


x  ■  Ax  +  Bu 


(116) 


y  -  Cl  (117) 

as  discussed  in  Chapter  III.  Thus,  the  ease  of  application  of  the  control 
method  heretofore  described  can  be  seen.  Now,  the  performance  of  control 
method  will  be  examined  using  the  programs'  eigenvalue  analysis  results. 
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VI.  Investigation 

Control  of  a  large  flexible  space  structure  Is  a  complicated  task  that 
Is  best  f one  step  at  a  time.  Therefore,  a  systematic,  building-block 
approach  was  used  to  conduct  the  study.  Miller's  work  concentrated  on  one 
and  two  controller  systems.  This  investigation  is  the  next  logical  step  in 
the  process,  concentrating  on  the  three  and  four  controller  systems. 

Outline 

The  initial  phase  of  this  study  concentrated  on  the  control  of  the 
CSDL  1  model  using  three  controllers  with  no  residual  terms.  The  first 
eight  structural  modes  were  used  in  this  analysis,  and  the  control  weight¬ 
ing  matrix  was  fixed  at  -  20  £  I  J. 

Next,  the  residuals  were  added  to  the  system,  utilizing  the  last  four 
structural  modes  of  th?  CSDL  1  model.  With  this,  the  full  twelve  mode 
model  was  implemented. 

The  residuals  were  then  incorporated  into  a  controller  to  form  a 
four  controller  system  without  residuals.  Again,  the  full  twelve  mode 
model  was  used  to  test  the  new  system. 

Once  satisfactory  results  were  obtained  using  the  CSDL  1  model, 
the  program  was  expanded  in  dimensions  to  accomodate  the  CSDL  2  model. 

A  twelve  mode  subset  of  the  modeled  modes  was  used  for  this  segment  of 
the  study.  The  selection  of  this  subset  of  modes  will  be  explained  in 
the  next  section.  Since  twelve  modes  were  used  in  both  models,  the  in¬ 
crease  in  dimensions  was  dictated  by  the  number  of  sensors  snd  actuators 
used  in  the  second  model.  Therefore,  all  matrices  were  expanded  to  a 
minimum  dimension  of  twenty-one,  as  can  be  seen  in  the  dimensioning  por¬ 
tion  of  the  program  in  Appendix  C. 

Zt  was  desired  to  control  all  twelve  nodes  of  the  CSDL  2  model. 


therefore,  no  residuel  modes  were  included  in  this  phase  of  the  study 


Three  controllers  were  used  with  a  control  weighting  matrix  of  * 


To  obtain  desired  control  characteristics ,  the  values  of  the 
control  weighting  matrix  were  varied  and  the  resulting  eigenvalue  move¬ 
ments  were  observed.  Three  of  the  twelve  modes  are  rigid  body  modes, 
therefore  the  weightings  on  these  modes  were  varied  first,  keeping 
the  value  of  20  for  all  of  the  flexible  modes.  Then  the  values  of  the 
weightings  for  the  flexible  modes  were  varied  individually  to  see  what 
effect  this  had  on  the  system's  controllability  and  observability. 

Finally,  the  CSDL  2  model  was  run  using  the  four  controller  system. 
The  control  weightings  were  varied  just  as  in  the  three  controller  system 
to  determine  if  there  were  any  significant  changes  in  the  system's  con¬ 
trollability  or  observability.  The  results  of  the  preceeding  outline 
will  be  presented  shortly. 

Modal  Selection  and  Grouping 

One  of  the  more  difficult  tasks  in  controlling  a  large  space  struc¬ 
ture  is  the  determination  of  which  structural  modes  are  to  be  actively 
controlled  and  which  are  left  as  residuals.  Factors  which  affect  this 
selection  process  include  sensor/actuator  placement  and  alignment, 
controller  bandwldths,  and  control  constraints  such  as  llne-of -sight 
tolerances. 

For  multiple  controllers,  an  additional  step  has  to  be  taken.  This 
is  the  assignment  of  the  modes  to  be  controlled  to  minimise  the  control 
effort.  Organising  modes  into  compatible  groupe  can  be  done  by  a 
simple  examination  of  the  angles  between  the.  vectors  of  nodal  amplitudes 
(the  angles  between  the  rows  of  the  •  D  matrix).  Defining  these  vectors 


5S 


as  ,  the  angles  aay  be  found  from  the  equation  for  the  dot  product  of 
two  vectors,  given  by 


*1  •  ■ 

♦i 

cos  0 


ij 


(118) 


The  nodes  are  then  grouped  so  no  two  orthogonal  nodes  are  in  the  sane 
controller.  In  sone  cases  the  grouping  of  orthogonal  nodes  in  the  sane 
controller  is  unavoidable.  Angles  snaller  than  +45  degrees  were  used 
in  this  study  when  possible,  however,  in  sone  cases  the  Unit  was 
extended  to  +70  degrees. 

CSDL  1.  Selection  of  the  nodes  to  be  controlled  in  the  first  nodel 
was  not  a  difficult  task  as  there  are  only  twelve  nodeled  nodes  for  the 
structure.  As  noted  earlier,  the  first  eight  nodes  were  originally 
selected  for  active  control.  The  last  four  were  initially  left  as 
residuals,  but  were  later  used  as  controlled  nodes.  Hence,  the  full 
systen  was  eventually  controlled  with  no  residual  nodes.  The  assign— nt 
of  nodes  to  controllers  was  nade  based  on  the  relative  angle  node  angles 
given  in  Table  X.  The  actual  groupings  are  presented  case  by  case  in  the 
next  section. 


CSDL  2.  The  selection  of  nodes  for  the  second  nodel,  on  the  other 
hand,  was  not  a  simple  task.  There  art  well  over  one  hundred  catalogued 
■odes  at  present,  however,  only  a  —all  subset  of  the  nodeled  nodes 
were  selected  for  active  control.  The  first  forty -four  nodes  are 
generally  used  to  denonstrate  the  structural  behavior  and  a  subset  of 
twelve  nodes  was  selected  fron  these.  The  specific  nodes  were  chosen 
as  a  result  of  a  study  conducted  by  Lockheed  and  sponsored  by  tons  Air 
Develop— nt  Center  (Xef  10).  This  study  used  a  High  Authority  Control 
— thod  and  based  its  nodal  selection  upon  the  line-of -sight  and  defocus 
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Te»  le  I 

Relative  Angles*  Between  CSDL  1  Modes 


Mod* 

1 

2 

3 

4 

5 

6 

1 

0.00 

90.00 

90.00 

66.40 

90.00 

90.00 

2 

90.00 

0.00 

64.33 

90.00 

96.61 

90.00 

3 

90.00 

64.22 

0.00 

90.00 

85.25 

90.  ($0 

4 

66.40 

90.00 

90.00 

0.00 

90.00 

90.00 

5 

90.00 

96.61 

85.25 

90.00 

0.00 

90.00 

6 

90.00 

90.00 

90.00 

90.00 

90.00 

0.00 

7 

33.23 

90.00 

90.00 

99.63 

90.00 

90.00 

8 

90.00 

50.12 

111.44 

90.00 

84.38 

90.00 

9 

90.00 

147.39 

145.07 

90.00 

77.97 

90.00 

10 

51.52 

90.00 

90.00 

14.88 

90.00 

90.00 

11 

90.00 

84.91 

61.82 

90.00 

23.47 

90.00 

12 

90.00 

108.20 

89.58 

90.00 

11.63 

90.00 

Mode 

7 

8 

9 

10 

11 

12 

1 

33.23 

90.00 

90.00 

51.52 

90.00 

90.00 

2 

90.00 

50.12 

149.39 

90.00 

84.91 

108.20 

3 

90.00 

111.44 

145.07 

90.00 

61.82 

89.58 

4 

99.63 

90.00 

90.00 

14.88 

90.00 

90.00 

5 

90.00 

84.38 

77.97 

90.00 

23.47 

11.63 

6 

90.00 

90.00 

90.00 

90.00 

90.  GO 

90.00 

7 

0.00 

90.00 

90.00 

84.75 

90.00 

90.00 

8 

90.00 

0.00 

97.28 

90.00 

92.36 

92.83 

9 

90.00 

97.28 

0.00 

90.00 

98.79 

68.88 

10 

84.75 

90.00 

90.00 

0.00 

90.00 

22.09 

11 

90.00 

92.36 

98.79 

90.00 

0.00 

30.09 

12 

90.00 

92.83 

68.88 

22.09 

30.09 

0.00 

**11  eagles  in  degrees 
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Table  IX 


Relative  Angles*  Between  CSDL  2  Modes 


Mode 

4 

5 

6 

7 

12 

13 

4 

0.00 

89.09 

90.00 

90.01 

90.09 

123.20 

5 

89.99 

0.00 

67.47 

175.39 

116.85 

90.05 

6 

90.00 

67.47 

0.00 

115.31 

0.00 

89.99 

7 

90.00 

175.39 

115.31 

0.00 

60.67 

89.69 

12 

90.00 

116.85 

104.01 

60.67 

0.00 

89.69 

13 

123.20 

90.05 

89.99 

39.95 

89.69 

0.00 

17 

107.48 

89.97 

90.04 

90.02 

90.02 

61.22 

21 

90.05 

106.57 

116.99 

73.19 

102.77 

90.07 

22 

90.01 

116.23 

91.57 

66.49 

68.57 

89.87 

28 

39.16 

87.89 

89.41 

92.45 

96.61 

104.90 

30 

81.12 

87.94 

89.27 

92.35 

95.63 

73.29 

Mode 

|  17 

21 

22 

24 

28 

30 

4 

104.48 

90.05 

90.01 

90.01 

39.16 

81.12 

5 

89.97 

106.57 

116.23 

90.16 

87.89 

87.84 

6 

90.04 

116.99 

91.57 

80.13 

89.41 

89.27 

7 

90.02 

73.19 

66.49 

90.68 

92.45 

92.35 

12 

90.01 

102.77 

68.57 

111.87 

96.61 

95.63 

13 

61.22 

90.07 

89.87 

90.15 

104.90 

73.29 

17 

0.00 

89.91 

90.15 

89.96 

80.56 

127.80 

21 

89.91 

0.00 

93.88 

38.05 

93.50 

93.10 

22 

90.15 

93.88 

0.00 

103.22 

90.51 

90.74 

24 

89.96 

38.05 

103.22 

0.00 

94.15 

93.50 

28 

80.56 

93.50 

90.51 

94.15 

0.00 

89.82 

30 

127.80 

93.10 

90.74 

93.50 

89.92 

0.00 

*all  angles  in  degrees 
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equations  given  by  a  TRW  study  (Ref  12).  The  line-of-sight  and  de¬ 
focus  algorithm  is  presented  in  Appendix  E. 

The  modes  used  in  this  investigation  contain  three  rigid  modes — 
labeled  4,  5,  and  6 — and  9  flexible  body  modes — labeled  7,  12,  13, 

17,  21,  22,  24,  28,  and  30.  These  modes  were  selected  by  Lockheed  as 
having  the  greatest  impact  on  line-of-sight  and  defocus,  based  on  che 
equations  given  in  Appendix  E,  and  so  were  adopted  for  this  study. 

Applying  the  relationship  given  by  Eq  118,  the  angles  between  the 
twelve  modes  are  given  in  Table  XI.  The  actual  groupings  selected  will 
be  presented  case  by  case  in  the  following  section. 


f 
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VII.  Results 


For  the  CSDL  1  model,  the  open  loop  damping  applied  was  0.005. 
Earlier  studies  of  the  model  showed  that  a  closed  loop  damping  of  0.10 
on  each  mode  was  needed  to  meet  pointing  requirements.  Therefore,  a 
minimum  of  ten  percent  damping  became  the  desired  parameter  for  accept¬ 
able  system  performance.  Due  to  unresolved  problems  with  the  computer 
subroutine  for  a  time  response,  line-of-sight  performance  could  not  be 
evaluated,  hence,  the  closed  loop  damping  is  the  only  numerical  perfor¬ 
mance  parameter  available. 

Calico  and  Miller  determined  that  a  control  weighting  matrix  of 

f\  1 


Q  *  20 


I 


was  sufficient  to  meet  the  ten  percent  damping  require- 


L  'J 

ment  for  one  and  two  controller  systems.  This  Q  matrix  was  adopted  as  an 
initial  value  in  this  study. 

From  Table  X,  the  most  favorable  modal  groupings  for  a  three  con¬ 
troller  are  given  by 

Group  1:  2,  3,  8,  9 

Group  2:  5,  6,  11,  12 

Group  3:  1,  A,  7,  10 

Mode  6  was  found  to  be  mutually  orthogonal  with  every  other  mode.  Its 
placement  is  based  on  the  recommendation  that  the  last  four  modes  be 
treated  as  residuals.  Without  mode  6,  the  second  group  would  consist  of 
mode  5  alone.  Therefore  mode  6  is  placed  in  the  second  group  and  the 
resulting  groupings  are 
Group  1:  2,  3,  8 

Group  2:  5,  6 

Group  1*  4.  7 
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Residuals:  9,  10,  11,  12 


The  first  test  was  performed  without  residuals  to  insure  that  the 
transformation  portion  of  the  program  would  successfully  block  triangular- 
ize  the  system  through  spillover  elimination.  The  resulting  eigenvalue 
analysis  is  presented  in  Table  XII.  The  overall  system  eigenvalues 
are  given  and  are  arranged  by  their  respective  controller  assignments. 

The  individual  modes  may  be  identified  by  the  imaginary  parts  of  the  eigen¬ 
values,  as  these  are  approximately  equal  to  the  modal  frequencies.  This 
particular  system  was  stable  and  very  well  damped  (average  damping  of 
thirty  percent)  before  transformation.  After  the  transformation,  stabi¬ 
lity  was  maintained,  as  expected,  and  there  was  a  notable  loss  in  damping 
on  three  modes.  There  was  an  overall  movement  of  eigenvalues  to  the 
right,  meaning  after  the  transformation  there  was  less  system  stability. 
Nevertheless,  the  transformation  succeeded  in  reducing  the  system  to  a 
block  triangular  form,  which  was  the  intent  of  this  step. 

It  should  be  pointed  out  that  the  results  given  were  obtained  using 
an  upper  block  triangular  transformation.  A  lower  block  triangular  trans¬ 
formation  yields  the  same  overall  system  results,  but  the  controller 
eigenvalues  are  presented  differently.  When  examining  the  A  +  BG  and 
A  -  KC  eigenvalues  for  each  controller,  those  eigenvalues  for  A  +BG 
in  an  upper  triangular  system  are  equal  to  those  for  A  -  KC  in  a  lower 
triangular  system.  Likewise,  those  eigenvalues  for  A  -  KC  in  an  upper 
triangular  system  are  equal  to  those  for  A  +  BG  in  a  lower  triangular 
system.  Table  XIII  presents  the  separate  controller  eigenvalues  for 
the  system  given  in  Table  XII.  The  upper  and  lower  block  triangular 
systems  are  presented  side-by-side  to  show  controller  relationships. 
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Table  XII 

CSDL  1  Overall  Eigenvalue  Analysis  -  3  Controllers 
Modal  Assignment 

Controller  1:  2,  3,  8  Controller  3:  1,  4,  7 

Controller  2:  5,  6  Residuals:  None 


Overall  System  Eigenvalues 

Before  Transformation  After  Transformation 

Controller  1 


-1.5814 

+  5.591861 

5 

M 

0.283 

-1.4503  + 

5.526401 

Z 

a 

0.262 

-1.4394 

+  5.369981 

Z 

a 

0.268 

-1.5096  + 

5.502781 

Z 

a 

0.274 

-.82780 

+  2.995711 

5 

a 

0.276 

-0.0895  + 

2.963611 

t 

a 

0.030 

-1.0125 

+  2.720861 

Z 

a 

0.372 

-0.9389  + 

2.840611 

c 

a 

0.319 

-.44731 

+  1.506291 

Z 

a 

0.297 

-0.4183  + 

1.447051 

4 

a 

0.289 

-.54119 

+  1.395841 

4 

m 

0.387 

-0.4969  + 

1.442171 

C 

a 

0.345 

Controller  2 

-.16127  +  4.909651  z  -  0.328  -1.6127  +  4.909651  Z  -  0.328 
-1.6126  +  4.909621  z  -  0.328  -1.3253  +  4.876321  z  -  0.272 
-1.2041+3.837021  c  -  0.314  -0.7383+3.872761  ;«  0.191 
-1.3672  +  3.488991  z  -  0.392  -0.0192  +  3.848341  z  -  0.005 


Controller  3 


-.15673  +  5.460781 

C  -  0.287 

-1.5672  +  5.460751 

C  -  0.287 

-1.567‘>  +  5.460751 

C  -  0.287 

-1.3751  +  5.45811 

Z  -  0.252 

-1.0420  +  3.426331 

C  -  0.304 

-0.5452  +  3.558661 

C  -  0.153 

-1.0420  +  3.426321 

C  -  0.304 

-1.0419  +  3.426321 

Z  -  0.304 

-.34014  +  1.163961 

?  -  0.292 

-0.1766  +  1.169331 

Z  -  0.151 

-.34014  +  1.163961 

Z  m  0.292 

-0.3401  +  1.163961 

Z  m  0.292 
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Table  XXII 


Upper  and  Lower  Triangular  Tranaformatlon 
Comparison  of  Controller  Eigenvalues 
Modal  Assignment 

Controller  1:  2,  3,8  Controller:  1,  4,  7 

Controller  2:  5,  6  Residual:  None 


Transformed  Controller  Eigenvalues 


Upper  Block  Triangular 


Lower  Block  Triangular 


A  +  BG  1 


A  +  BG  1 


-1.4503  +  5.526401 
-0.0895  +  2.963611 
-0.4183  +  1.447051 

A  -  KC  1 

-1.5096  +  5.502781 
-0.9389  +  2.840611 
-0.4970  +  1.442181 


-1.5096  +  5.502781 
-0.9389  +  2.840611 
-0.4970  ±  1.442181 

A  -  KC  1 

-1.4503  +  5.526401 
-0.0895  +  2.963611 
-0.4182  +  1.447051 


A  +  BG  2 


A  +  BG  2 


-1.3253  +  4.876321 
-0.7383  +  3.872761 

A  -  KC  2 

-1.6127  +  4.909641 
-0.0192  +  3.848341 


-1.6127  +  4.909641 
-0.0192  +  3.848341 

A  -  KC  2 

-1.3253  +  4.876321 
-0.7383  +  3.872761 


A  +  BG  3 


A  +  BG  3 


-1.5672  +  5.460751 
-1.0419  +  3.426321 
-0.3401  +  1.163961 

A  -  KC  3 

-1.3751  +  5.458111 
-0.5452  +  3.558661 
-0.1766  +  1.169331 


-1.3751  +  5.458111 
-0.5452  +  3.558661 
-0.1766  +  1.169331 

A  -  KC  3 

-1.5672  +  5.460731 
-1.0419  +  3.426321 
-0.3401  1.163961 
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Next,  Che  four  residual  nodes  are  added  to  complete  Che  aodel.  Ic 
was  desired  Co  see  Che  novenenC  of  Che  residuals,  if  any,  caused  by  Che 
spillover  elininacion.  The  eigenvalue  analysis  for  chis  seep  is  pre- 
senced  in  Table  XIV.  As  before,  chere  is  a  slight  sacrifice  in  che  closed 
loop  danping  during  Che  CransfomaCion,  buC  only  evo  nodes  exhiblced 
losses  below  Che  Cen  percent  nark.  The  first  of  chese  is  due  to  a 
decrease  in  controllability  in  Che  first  controller,  which  is  expected 
to  occur.  The  second  is  believed  Co  be  a  result  of  che  pairing  of  modes 
S  and  6  in  che  sane  controller.  This  is  supported  b>  several  different 
groupings  with  and  without  nodes  5  and  6  together.  The  drastic  reduc- 
Cion  in  the  controllability  and/or  observability  in  mode  S  is  not  seen 
when  modes  5  and  6  are  in  separate  controllers.  Therefore,  this  indi¬ 
cates  a  bad  grouping  of  nodes,  but  system  stability  was  not  totally 
sacrificed.  Overall,  the  requiresnnt  of  ten  percent  damping  in  all 
nodes  was  net.  The  residuals,  although  less  stable,  did  not  become 
unstable,  and  maintained  the  0.005  danping  originally  applied  to  the 
system. 

To  maintain  the  parallel  exanination,  the  same  modal  grouping  was 
then  controlled  using  four  controllers.  As  before,  damping  applied 
was  0.005  and  the  control  weighting  natrix  was  Q  ■  20  £  I  ^  .  The 

resulting  eigenvalue  analysis  is  given  in  Table  XV.  Again,  a  loss  in 
danping  occurs  during  transformation  for  the  first  three  controllers, 
but  the  fourth  controller  actually  show  sons  Improvement  in  danping 
on  the  nodes  which  were  previously  residuals.  Also,  for  the  first  tine 
during  this  study,  the  danping  on  one  node  in  the  first  controller 
dropped  below  the  original  open  loop  damping,  this  being  node  3.  The 
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Table  XIV 


CSDL  1  Overall  Eigenvalue  Analysis  -  3  Controllers 
Modal  Assignment 

Controller  1:  2,  3,  8  Controller  3:  1,  4,  7 

Controller  2:  5,  6  Residual:  9,  10,  11,  12 

Overall  System  Eigenvalues 

Before  Transformation  After  Transformation 

Controller  1 


-1.5082  + 

5.608771 

C 

as 

0.269 

-1.3875 

+ 

5.694181 

r 

•* 

a 

0.244 

-1.1542  + 

5.364101 

c 

a 

0.282 

-.15842  + 

5.362871 

c 

SB 

0.295 

-1.2334  ±  3.187031 

c 

a 

0.387 

-0.0903 

2.960691 

c 

a 

0.030 

-1.1119  + 

2.601291 

SB 

0.427 

-0.9482 

+ 

2.81358* 

c 

m 

0.337 

-0. 3482  + 

1.586071 

c 

mi 

0.219 

-0.3555 

+ 

1.52676* 

c 

m 

0.233 

-0.6579  + 

1.350471 

c 

a 

0.487 

-0.5482 

+ 

1.383561 

c 

m 

0.396 

Controller  2 

-1.6126  + 

4.909651 

c 

- 

0.328 

-1.6127 

+ 

4.909651 

c 

- 

0.328 

-1.6126  + 

4.909621 

c 

- 

0.328 

-1.3253 

+ 

4.876321 

c 

- 

0.339 

-1.0697  + 

4.013421 

c 

m 

0.267 

-0.0192 

+ 

3.843841 

c 

- 

0.005 

-0.8339  + 

3.772431 

c 

a 

0.221 

-1.0692 

+ 

3.35319* 

c 

V 

0.319 

Controller  3 

-1.5124  + 

5.512871 

c 

m 

0.274 

-1.6097 

+ 

5.522971 

c 

a 

0.291 

-1.6183  ± 

5.412921 

c 

a 

0.299 

-1.3595 

+ 

5.388131 

c 

SB 

0.252 

-1.4846  + 

3.404981 

c 

a 

0.436 

-0.7383 

+ 

3.872761 

c 

m 

0.191 

-0.6912  ± 

3.202521 

c 

m 

0.216 

-0.5112 

+ 

3.664041 

t 

m 

0.140 

-0.2689  + 

1.220061 

f. 

m 

0.220 

-0.3381 

+ 

1.187361 

c 

m 

0.285 

-0.4152  ± 

1.123041 

c 

m 

0.370 

-0.1719 

+ 

1.137051 

c 

SB 

0.151 

Residual 

-0.0701  +  13.96481 

5 

a 

0.005 

-0.0699 

+ 

13.96641 

c 

a 

0.005 

-0.0573  + 

10.8978* 

5 

a 

0.005 

-0.0548 

10.9231* 

t 

a 

0.005 

-0.0679  +  10.25741 

c 

a 

0.005 

-0.0517 

+ 

10.30051 

c 

a 

0.005 

-0.0768  +  8.883121 

c 

- 

0.009 

-0.0451 

+  8.92917* 

c 

* 

0.005 
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Table  XV 

CSOL  1  Overall  Eigenvalue  Analysis  -  4  Controller 
Modal  Assignments 


Controller  1:  2,  3,  8 

Controller  2:  5,  6 


Controller  3:  1,  4,  7 

Controller  4:  9,  10,  11 , 


Overall  System  Eigenvalues 

Before  Transformation  After  Transformation 

Controller  1 


-1.5834 

+ 

5.793201 

4 

as 

0.273 

-1.4838 

+ 

5. 62"27i 

4 

-1.4901 

+ 

5.557741 

4 

■ 

0.268 

-1.4925 

+ 

5.345871 

4 

-1.1962 

+ 

3.375961 

4 

m 

0.354 

-0.0120 

+ 

3.851571 

4 

-1.1160 

+ 

2.783591 

4 

m 

0.401 

-0.9282 

+ 

2.790361 

4 

-0.3421 

+ 

1.666201 

4 

m 

0.205 

-0.3831 

+ 

1.595111 

4 

-0.6071 

+ 

1.379051 

4 

m 

0.440 

-0.4653 

+ 

1.424341 

4 

Controller  2 

-1.6127 

+ 

4.909651 

4 

- 

0.328 

-1.6126 

+ 

4.909651 

4 

m 

-1.6126 

+ 

4.909641 

4 

« 

0.328 

-1.6126 

+ 

4.909601 

4 

m 

-1.0361 

+ 

5.150581 

4 

- 

0.201 

-0.0192 

+ 

3.848341 

4 

m 

-1.0151 

+ 

4.244181 

4 

- 

0.239 

-0.3599 

+ 

3.607571 

4 

m 

Controller  3 


-1.6365  +  5.393571 

4  -  0.303 

-1.3751  + 

5.458111 

4  . 

-1.7358  +  5.164751 

4  -  0.336 

-1.5566  + 

5.455551 

4  - 

-1.3556  +  3.644651 

4  -  0.372 

-0.2324  + 

3.559321 

4  • 

-0.9042  +  3.618431 

4  -  0.250 

-0.5452  + 

3.558661 

4  > 

-0.2572  +  1.244191 

4  -  0.207 

-0.1766  + 

1.169331 

4  - 

-0.4006  T  1.126511 

4  -  0.356 

-0.2668  5 

1.1667U 

4  - 

Controller  4 


-0.5085  +  13.99421 

4  -  0,036 

-0.3065  ±  13.97161 

4  - 

-0.2441  +  13.74051 

4  -  0.018 

-0.il37  ±  13.94101 

4  ■ 

-1.2226  +  11.054H 

4  -  0.111 

-0.9592  ±  10.9018* 

4  - 

-1.5415  +  10.36221 

4  -  0.149 

-0.2197  ±  10.8734* 

4  - 

-0.5054  +  10.10031 

4  -  0.050 

-0.0515  t  10.30351 

4  - 

-0.8318  +  9.520681 

4  -  0.087 

-1.2352  ±  10.23071 

4  - 

-1.6044  +  8.979221 

4  -  0.179 

-1.5900  ±  9.158701 

4  - 

-0.8934  +  8.084211 

4  -  0.110 

-1.0199  ±  8.328721 

4  “ 

12 


0.264 

0.279 

0.003 

0.333 

0.240 

0.327 


0.328 

0.328 

0.005 

0.100 


0.252 

0.281 

0.065 

0.153 

0.151 

0.229 


0.026 

0.008 

0.088 

0.020 

0.005 

0.121 

0.174 

0.122 
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efface*  of  Chi*  cannot  be  determined  until  a  time  history  of  system 
response  can  be  generated.  The  other  mode  seriously  affected  vas  mode 
5  again.  This  is  an  expected  result,  as  mentioned  earlier.  However, 
the  overall  results  from  this  four  controller  run  were  surprising,  since 
the  requirement  was  specified  in  an  earlier  derivation  that  the  number 
of  modes  being  suppressed  could  not  exceed  the  number  of  sensors  avail¬ 
able.  Decoupling  controller  1  in  this  case  Involves  suppressing  nine 
modes  with  a  six  sensor  system.  This  is  the  reason  for  the  total  loss 
in  damping  and  probably  accounts  in  large  part  for  the  overall  decrease 
in  closed  loop  damping.  In  all  other  modal  groupings,  the  four  con¬ 
troller  system  failed  to  stabilize  the  system  and  extremely  large,  posi¬ 
tive  eigenvalues  were  given  in  the  overall  system  analysis.  This  parti¬ 
cular  run  was  successful  only  on  a  chance  compatible  grouping. 

For  the  CSDL  2  model,  the  open  loop  damping  ratio  applied  was 
0.010,  twice  that  for  the  CSDL  1  model.  Since  a  multiple  controller  had 
never  been  applied  to  this  model  prior  to  this  study,  an  initial  control 


weighting  matrix  of  Q 


20 


\  ‘ 
I 

N 


was  used  for  a  first  "feel"  at 


controlling  the  system.  The  modal  angles  from  Table  XI  give  the  follow¬ 
ing  grouping  for  a  three  controller  system,  controlling  all  twelve  modes 
Group  1:  5,  6,  7,  21 

Group  2:  4,  13,  17,  30 

Group  3:  12,  22,  24,  28 

(recalling  that  the  moot*  selected  are  not  sequential).  Mode  28  is 
similar  to  mode  6  in  the  first  model  in  the  sense  that  it  is  mutually 
orthogonal  to  all  except  mode  4.  Therefore,  it  is  randomly  placed  with 
the  third  grouping  to  balance  the  modal  distribution.  Control  of  all 
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twelve  modes  la  desired,  therefore,  there  will  be  no  residual  nodes  in 
the  study  of  the  second  aodel.  As  in  the  first  model,  we  will  try  to 
achieve  a  closed  loop  damping  of  ten  percent  on  all  modes.  This  is 
a  simple  starting  point  for  testing  the  application  of  a  multiple  con¬ 
troller  to  the  given  model  and  is  in  no  way  definitive  as  to  controller 
success  or  failure. 

The  eigenvalue  analysis  for  the  first  run  using  the  modal  group¬ 
ings  above  and  Q  »  20  M  is  given  in  Table  XVI.  Although  no 
inherent  instabilities  exist,  it  is  evident  there  is  an  excessive 
amount  of  damping  on  the  rigid  body  modes  (over  one  hundred  percent) 
and  no  increase  at  all  in  the  flexible  modes.  In  some  instances,  there 
is  even  a  decrease  in  the  damping  during  the  transformation,  but  as 
can  be  seen,  no  mode  went  below  the  initial  open  loop  damping  of  one 
percent . 

The  next  step  was  to  equalise  the  damping  between  the  rigid  and 
flexible  modes.  This  was  atempted  by  raising  the  control  weighting 
on  the  flexible  modes  and  decreasing  it  on  the  rigid  body  modes.  In 
the  untransformed  closed  loop  system,  there  is  a  symmetry  which  exists 
only  as  a  result  of  the  colocation  of  sensors  and  actuators.  This 
symmetry  is  apparent  in  the  matrix  Riccatl  solutions  to  4- 
and  Ai  -  K^,  in  the  state  feedback  matrices  and  and  in  the 
controller  eigenvalues  of  A^  +  B^G^  and  A^  -  However,  when  the 

control  weightings  are  varied  as  described,  this  symmetry  is  completely 
lost  in  those  controllers  containing  both  rigid  and  flexible  modes. 

This  loss  of  symmetry  is  not  understood,  since  the  weighting  matrices 


Table  XVI 


CSOL  2  Overall  Eigenvalue  Analysis  -  3  Controllers 

q£  -  20 

Modal  Assignments 

Controller  Is  5,  6,  7,  21  Controller  3:  12,  22,  24,  28 

Controller  2:  4,  13,  17,  30  Residual:  None 

Overall  System  Eigenvalues 

Before  Transformation  After  Transformation 


Controller  1 


-.06463  +  6.107291 

4  - 

0.011 

-.06305  +  6.107131 

4 

ai 

0.010 

-.06495  +  6.106961 

4  - 

0.011 

-.06479  4-  6.107111 

c 

m 

0.011 

-.04041  +  0.717351 

4  - 

0.056 

-.03367  +  0.716161 

4 

m 

0.047 

-.04187  +  0.715031 

c  - 

0.058 

-.04115  +  0.716151 

4 

m 

0.057 

-.24749  +  0.242111 

c  - 

1.022 

-.25463  +  0.239891 

4 

m 

1.061 

-.26163  +  0.238291 

c  - 

1.098 

-.24170  +  0.228841 

c 

m 

1.056 

-.20131  +  0.199561 

4  - 

1.009 

-.20632  4-  0.198141 

4 

m 

1.041 

-.21119  +  0.196891 

c  - 

1.072 

-.18949  4-  0.183101 

c 

m 

1.035 

Controller  2 


-.25049 

4- 

25.04831 

C 

m 

0.010 

-.25049 

+ 

25.04831 

4 

m 

0.010 

-.25049 

+ 

25.04831 

C 

m 

0.010 

-.25049 

+ 

25.04831 

4 

m 

0.010 

-.05144 

4- 

5.121381 

c 

m 

0.010 

-.01543 

+ 

5.121381 

4 

m 

0.010 

-.05144 

+ 

5.121381 

c 

m 

0.010 

-.05144 

+ 

5.121381 

4 

m 

0.010 

-.08045 

4- 

3.745181 

c 

m 

0.021 

-.07811 

± 

3.745021 

4 

m 

0.021 

-.08042 

+ 

3.745181 

c 

m 

0.021 

-.08044 

+ 

3.745161 

4 

m 

0.021 

-.29226 

4- 

0.271471 

c 

m 

1.077 

-.28999 

+ 

0.268191 

4 

m 

1.081 

-.29512 

+ 

0.270791 

4 

m 

1.090 

-.22818 

+ 

0.217131 

4 

m 

1.051 

Controller  3 

-.21691 

+ 

21.69031 

4 

m 

0.010 

-.21691 

+ 

21.69031 

4 

m 

0.010 

-.21691 

+ 

21.69031 

4 

m 

0.010 

-.21691 

+ 

21.69031 

4 

m 

0.010 

-.14844 

4 

11.14031 

4 

m 

0.013 

-.12045 

+ 

11.14031 

4 

m 

0.011 

-.14763 

4- 

11.13951 

4 

m 

0.013 

-.14801 

■f 

11.13991 

4 

m 

0.013 

-.07731 

+ 

7.280261 

4 

m 

0.011 

-.07444 

+ 

7.280231 

4 

m 

0.010 

-.07710 

+ 

7.280071 

4 

m 

0.011 

-.07717 

+ 

7.280201 

4 

m 

0.011 

-.04368 

+ 

3.501921 

4 

m 

0.012 

-.03802 

+ 

3.501731 

4 

m 

0.011 

-.04277 

t 

3.501111 

4 

m 

0.012 

-.04313 

+ 

3.501681 

4 

m 

0.012 
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believed  to  be  a  numerical  incongruity  in  the  matrix  Riccati  solution 
subroutine,  MRIC  (Appendix  D),  but  this  is  not  confirmed. 

To  bypass  this  problem  until  its  effects  on  the  system  performance 
can  be  determined,  the  modal  groupings  were  rearranged  to  combine  the 
rigid  body  modes  into  one  controller.  As  a  result,  the  groupings  became: 
Group  1:  4,  5,  6 

Group  2:  7,  13,  17,  21,  30 

Group  3:  12,  22,  24,  28 

The  control  weightings  could  then  be  adjusted  freely  for  the  flexible 
and  rigid  body  modes  because  they  are  totally  decoupled.  The  rigid 
body  mode  control  weighting  was  adjusted  from  0.02  to  20  by  orders  of 
magnitude  and  the  flexible  body  mode  control  weighting  was  adjusted 
from  20  to  15000  by  approximate  doubles  of  the  previous  value.  A  repre¬ 
sentative  eigenvalue  analysis  is  presented  in  Table  XVII.  The  values 
of  are  presented  such  that 

i 

(119) 

a 

qj_ 

There  is  an  unstable  mode  present  with  this  modal  grouping  and  it 
exists  only  within  the  overall  system  since  the  individual  controllers 
are  stable.  But  the  transformation  succeeded  in  stabilizing  the  mode  and 
even  increased  the  damping  on  the  mode.  Damping  on  two  modes  achieved 
the  ten  percent  desired  (modes  7  and  13),  not  to  mention  the  rigid  body 
modes,  all  of  which  are  critically  damped.  Several  other  modes  had  signi¬ 
ficant  gains  in  their  damping  (modes  12,  21,  22  and  24).  However,  other 
modes  sacrlf lad  soma  of  their  damping  during  the  transformation,  therefore 


Table  XVII 


CSDL  2  Overall  Eigenvalue  Analysis  -  3  Controllers 


■  2-°  ’fl.xibL.  '  1000 
Modal  Assignments 


Controller  1:  4,  5,  6  Controller  3:  12,  22,  24,  28 

Controller  2:  7,  13,  17,  21,  30  Residuals:  None 

Overall  System  Eigenvalues 

Before  Transformation  After  Transformation 

Controller  1 


-.16355 

+ 

0.64915i 

4 

m 

0.252 

-.16093  + 

0.156911 

4 

* 

1.026 

-.58333 

+ 

0.579451 

4 

m 

1.007 

-.14062  + 

0.137921 

4 

m 

1.020 

-.11817 

+ 

0.244371 

4 

m 

0.484 

-.11454  + 

0.113071 

4 

m 

1.013 

-.16155 

+ 

0.156631 

4 

m 

1.031 

-.10480  + 

0.103671 

4 

m 

1.011 

-.14044 

+ 

0.137901 

4 

m 

1.018 

-.08908  + 

0.088381 

4 

m 

1.008 

-.11469 

+ 

0.133241 

4 

m 

1.013 

-.00248  + 

0.002481 

4 

m 

1.000 

Controller  2 

-.25049 

+ 

25.04831 

4 

m 

0.010 

-.25050  + 

25.04831 

4 

m 

0.010 

-.25049 

+ 

25.04831 

4 

m 

0.010 

-.25050  + 

25.04831 

4 

m 

0.010 

-.14827 

+ 

6.128691 

4 

m 

0.024 

-.12806  + 

6.106161 

4 

m 

0.021 

-.18024 

+ 

6.092721 

4 

m 

0.030 

-.14468  + 

6.105811 

4 

m 

0.024 

-.06326 

+ 

5.122681 

4 

m 

0.012 

-.06036  + 

5.121281 

4 

m 

0.012 

-.05924 

+ 

5.118961 

4 

m 

0.012 

-.06072  + 

5.121261 

4 

m 

0.012 

-.60535 

+ 

3.799601 

4 

m 

0.139 

-.42295  + 

3.725171 

4 

m 

0.114 

-.40053 

+ 

3.570231 

4 

m 

0.112 

-.48612  + 

3.718481 

4 

m 

0.131 

-.60203 

+ 

1.056931 

4 

m 

0.570 

-.23565  + 

0.727251 

4 

m 

0.324 

+.47643 

+ 

0.676371 

4 

m 

— 

-.02015  + 

0.716121 

4 

m 

0.028 

Controller  3 

-.21691 

+ 

21.69031 

4 

m 

0.010 

-.21691  + 

21.69031 

4 

m 

0.010 

-.21691 

+ 

21.69031 

4 

m 

0.010 

-.21691  + 

21.69031 

4 

m 

0.010 

-.72335 

+ 

11.13831 

4 

m 

0.065 

-.34237  + 

11.13581 

4 

m 

0.031 

-.67321 

+ 

11.08391 

4 

m 

0.061 

-.69781  + 

11.11941 

4 

m 

0.063 

-.21047 

+ 

7.290881 

4 

m. 

0.029 

-.13175  + 

7.279451 

4 

m 

0.018 

-.18001 

+ 

7.258771 

4 

m 

0.025 

-.19501  + 

7.278091 

4 

m 

0.027 

-.21479 

+ 

3.524611 

4 

m 

0.061 

-.11033  + 

3.500431 

4 

m 

0.032 

-.15069 

+ 

3.445141 

4 

m 

0.044 

-.18141  + 

3.497921 

4 

m 

0.052 
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there  is  a  trade-off  for  the  added  gain.  Regardless  of  the  weighting 
factor  applied,  modes  28  and  30  never  achieved  any  gain  in  damping. 

Table  XVIII  lists  the  eigenvalues  of  each  controller  for  the  run. 

The  instability  shown  in  the  overall  system  eigenvalues.  Table  XVII,  is 
absent  in  the  individual  controller.  This  is  a  perfect  example  of 
the  point  mentioned  repeatedly  in  Chapter  III  that  stable  controllers 
do  not  insure  a  stable  system.  Spillover  terms  do  have  a  noticeable 
effect.  This  situation  might  never  have  been  seen  in  this  study  if  the 
modal  groupings  had  not  been  changed. 

Other  items  of  note  from  Table  XVIII  include  the  display  of  con¬ 
troller  symmetry  mentioned  earlier.  In  the  untransformed  system  the 
control  and  observation  portions  of  each  controller  yield  identical 
eigenvalues.  The  loss  of  this  similarity  led  to  postponing  the  use  of 
different  control  weightings  for  each  mode.  Also  visible,  in  the  trans¬ 
formed  system,  is  the  loss  of  controllability  in  Controller  1  and  the 
loss  of  observability  in  Controller  3.  This  is  characteristic  of  all 
multiple  controllers.  If  a  lower  block  triangular  transformation  were 
used,  the  loss  of  controllability  would  be  in  Controller  3  and  the  loss 
of  observability  would  be  in  Controller  1.  This  may  be  seen  by  the  rela¬ 
tionship  of  the  two  systems  demonstrated  earlier  in  Table  XIII.  In  any 
case,  the  first  and  last  controller  in  any  multiple  controller  will  exper¬ 
ience  a  loss  of  controllability  or  observability. 

The  final  step  was  to  apply  the  four  controller  system  and,  again, 
attempt  to  achieve  ten  percent  closed  loop  damping  on  all  modes.  It 
was  deemed  best  to  keep  the  rigid  body  modes  in  one  controller  and 
distribute  the  remaining  modes  to  the  last  three  controllers.  Since 
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Table  XVIII 


CSDL  2  Controller  Eigenvalue  Analysis  -  3  Controllers 


^rigid  ^  ^flexible 

Modal  Assignments 


1000 


Controller  1:  44  5,  6 

Controller  2:  7,  13,  17,  21,  30 

Controller 

Before  Transformation 
A  +  BG  1 

-.16093  +  0.156911 
-.14062  +  0.137921 
-.11454  +  0.113071 

A  -  KC  1 

-.16903  +  0.1569U 
-.14062  +  0.137921 
-.11454  +  0.113071 

A  +  BG  2 

-.25050  +  25.04831 
-.16457  +  6.105331 
-.06103  +  5.121261 
-.50479  +  3.716321 
-.28750  +  0.731871 

A  -  KC  2 

-.25050  +  25.04831 
-.16457  +  6.105331 
-.06103  +  5.121261 
-.50479  +  3.716321 
-.29750  +  0.731871 

A  +  BG  3 

-.21691  +  21.69031 
-.69781  +  11.11941 
-.19501  +  7.278091 
-.18141  +  3.497911 

A  -  KC  3 

-.21691  +  21.69031 
-.69781  +  11.11941 
-.19501  +  7.278091 
-.18141  +  3.497911 


Controller  3:  12,  22,  24,  28 

Residual:  None 

Eigenvalues 

After  Transformation 

A  +  BG  1 

-.10480  +  0.103671 
-.08908  +  0.088381 
-.00248  +  0.002481 

A  -  KC  1 

-.16093  +  0.156911 
-.14062  +  0.137921 
-.11454  +  0.113071 

A  +  BG  2 

-.25050  +  25.04831 
-.12806  +  6.106161 
-.06072  +  5 . 121261 
-.48612+3.718481 
-.23565  +  0.727251 

A  -  KC  2 

-.25050  +  25.04831 
-.14468  +  6. 105811 
-.06036  +  5.121281 
-.42295  +  3.725171 
-.02015  +  0.716121 

A  +  BG  3 

-.21691  +  21.69031 
-.69781  +  11.11941 
-.19501  T  7.278091 
-.18141  +  3.497921 

A  -  KC  3 

-.21691  +  21.69031 
-.34237  +  11.13571 
-.13175  +  7.279441 
-.11033  +  3.500431 
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Table  XIX 


CSDL  2  Overall  Eigenvalue  Analysis  -  4  Controllers 


q  .  . ,  “2.0 

rigid 


q ...  -  500 

flexible 


Modal  Assignments 


Controller  1:  4,  5,  6 

Controller  2:  13,  17,  30 


Controller  3:  12,  22,  24,  28 

Controller  4:  7,  21 


Overall  System  Eigenvalues 


Before  Transformation 


Controller  1 


After  Transformation 


-.4836 

+  0.617861 

4 

» 

0.783 

-.1609  + 

0.156911 

4 

3 

1.025 

-.1619 

+  0.452021 

4 

3 

0.358 

-.1406  + 

0.137921 

4 

3 

1.019 

-.1183 

+  0.177301 

4 

s 

0.667 

-.1145  + 

0.113071 

4 

» 

1.013 

-.1616 

+  0.156671 

4 

3 

1.031 

-.1048  + 

0.103671 

4 

s 

1.011 

-.1399 

+  0.138151 

4 

= 

1.013 

-.0891  + 

0.08838 i 

4 

m 

1.008 

-.1150 

+  0.113531 

4 

= 

1.012 

-.0025  + 

0.002471 

4 

“ 

1.000 

Controller  2 


-.2505 

+ 

25.04831 

4 

= 

0.010 

-.2505  + 

25.04831 

4 

B 

0.010 

-.2505 

+ 

25.04831 

4 

s 

0.010 

-.2505  + 

25.04821 

4 

- 

0.010 

-.0572 

+ 

5.121951 

4 

s 

0.011 

-.0560  + 

5.121331 

4 

3 

0.011 

-.0556 

+ 

5.120421 

4 

0.011 

-.0560  + 

5.121321 

4 

= 

0.011 

-.4162 

+ 

3.777561 

4 

* 

0.110 

-.3002  + 

3.735471 

4 

3 

0.080 

-.2979 

+ 

3.654451 

4 

3 

0.082 

-.3443  + 

3.732171 

4 

3 

0.092 

Controller  3 


-.2169 

+ 

21.69031 

4 

- 

0.010 

-.2169  + 

21.69031 

4 

3 

0.010 

-.2169 

+ 

21.69031 

4 

3 

0.010 

-.2169  + 

21.69031 

4 

= 

0.010 

-.5136 

+ 

11. 1409i 

4 

3 

0.046 

-.3791  + 

11.13461 

4 

- 

0.034 

-.4861 

+ 

11.11221 

4 

3 

0.044 

-.4795  + 

11. 1308i 

4 

« 

0.043 

-.1546 

+ 

7.285341 

4 

3 

0.021 

-.1345  + 

7.279401 

4 

3 

0.018 

-.1399 

+ 

7.270241 

4 

3 

0.019 

-.1366  + 

7.279361 

4 

3 

0.019 

-.1488 

+ 

3.514131 

4 

3 

0.042 

-.1080  + 

3.500491 

4 

3 

0.031 

-.1134 

+ 

3.473691 

4 

- 

0.033 

-.1107  + 

3.500421 

4 

3 

0.032 

Controller  4 


-.1164 

+  6.116561 

4 

“  0.019 

-.0745  + 

6.107011 

4 

“  0.012 

-.1317 

+  6.099931 

4 

-  0.022 

-.1241  + 

6.106241 

4 

“  0.020 

-.4323 

+  0.850261 

4 

“  0.508 

-.2064  + 

0.724851 

4 

-  0.285 

-.3671 

+  0.632771 

4 

-  0.5,80 

-.0072  + 

0.716041 

4 

-  0.010 

the  three  controller  application  had  no  residuals,  a  direct  conversion 
of  the  residuals  to  a  fourth  controller  was  not  possible,  therefore  a 
little  modal  rearrangement  was  in  order.  To  maintain  a  similarity  to 
the  previous  tests,  two  modes  were  shifted  from  the  second  group  to 
form  a  fourth  group.  The  final  groupings  were  then  given  by 
Group  1:  4,  5,  6 

Group  2:  13,  17,  30 

Group  3:  12,  22,  24,  28 

Group  4:  7,  21 

Using  a  control  weighting  of  2.0  for  the  rigid  body  modes  and  500  for 
the  flexible  body  modes,  the  overall  eigenvalue  analyses  for  a  four  con¬ 
troller  system  are  given  in  Table  XIX.  As  in  the  previous  case,  the 
rigid  body  damping  was  increased  for  the  rigid  body  modes.  In  the 
flexible  modes,  no  appreciable  loss  in  damping  occured,  except  in  the 
fourth  controller  where  the  expected  loss  in  observability  was  found. 

A  blanket  ten  percent  in  damping  was  unachieveable  by  uniform  increases 
in  the  control  weighting,  again  suggesting  a  reexamination  of  the  modal 
groupings.  Also,  modes  28  and  30  were  unaffected  by  any  value  of 


from  10  to  15000. 


VIII.  Conclusions 


This  investigation  demonstrated  the  feasibility  of  using  multiple 
controllers  in  maintaining  system  stability.  It  was  shown  that  some 
closed  loop  damping  was  sacrificed  in  transforming  the  system  to  a 
block  triangular  form,  however,  the  effects  of  that  loss  have  not  been 
examined.  Modal  grouping  played  a  very  important  part  in  the  system 
stability  achieved  and,  in  the  four  controller  system,  allowed  an  other¬ 
wise  uncontrollable  configuration  (CSDL  1)  to  be  stabilized. 

The  use  of  angles  between  modal  amplitude  vectors  is  a  convenient 
method  for  initial  grouping  of  modes,  but  the  rank  of  the  B  and  C 
matrices  should  be  examined  closely.  If  they  are  not  of  full  rank,  the 
modal  groups  may  not  be  fully  compatible  as  indicated  by  zero  entries 
in  the  non-zero  singular  values  of  the  singular  value  decomposition. 

Loss  of  controllability  and  observability  in  the  first  and  last 
controllers  become  more  noticeable  with  the  addition  of  more  controllers. 
This  may  be  the  result  of  the  modal  assignments  used. 

The  inability  to  affect  the  last  two  modes  of  the  CSDL  2  model 
suggest  that  the  sensor  and  actuator  placement  may  not  be  suitable  for 
controlling  those  modes.  This  may  be  resolved  by  repositioning  the 
sensors  available  or  adding  sensors. 

For  those  cases  involving  residuals,  the  residuals  were  not  ser¬ 
iously  destabilized,  although  their  general  movement  was  to  become 
less  stable.  This  movement  to  the  right  is  contrary  to  what  was 
desired  and  needs  more  study. 

Overall,  the  goal  of  effecting  ten  percent  closed  loop  damping  was 
successful  on  the  CSDL  1  model  and  unsuccessful  on  the  CSDL  2  model. 


although  significant  increases  were  obtained  on  moat  modes.  Total 
controller  decoupling  was  achieved  while  maintaining  controller  and 
system  stability  for  both  models.  The  performance  of  the  controller, 
as  well  as  could  be  determined  without  running  a  forty-four  mode  simu¬ 
lation,  was  satisfactory,  but  did  not  meet  all  expectations. 


IX.  Recommendations 

There  are  several  directions  that  may  be  taken  at  this  point.  The 
intent  of  this  study  was  to  examine  the  application  and  performance  of 
three  and  four  decentralized  controllers  on  the  CSDL  2  model,  using  the 
CSDL  1  model  as  a  check  on  the  controller  algorithm.  A  re-examination 
of  the  modal  assignments  is  in  order  to  find  grouping  which  are  more 
compatible.  This  compatibility  may  be  determined  by  receiving  non¬ 
zero  singular  values  for  the  singular  value  decomposition  or  simply 
from  the  fullness  of  rank  of  the  matrices.  Another  direction  that 
may  be  taken  is  to  consider  adding  a  sensor  specifically  for  modes  28 
and  30  or  redistributing  the  existing  sensors  to  observe  these  two 
modes  more  directly.  These  changes  are  suggested  to  improve  the  obser¬ 
vability  and  controllability  of  the  existing  system. 

A  time  history  of  the  controller  response  would  be  invaluable  at 
this  time,  as  this  is  the  major  performance  criteria.  This  investigator 
was  unable  to  complete  such  a  response.  Additionally,  the  program  used 
may  be  examined  for  means  to  minimize  core  memory  requirements.  Fin¬ 
ally,  the  controllers  may  be  expanded  to  run  a  higher  number  of  modes. 
The  first  forty-four  modes  are  usually  used  as  a  fair  system  representa¬ 
tion.  These  suggestions  expand  upon  the  work  done  to  date. 

All  of  the  above  are  either  necessary  or  desirable.  At  a  minimum, 
these  should  be  accomplished  before  the  feasibility  of  implementing 
this  system  can  be  fully  evaluated. 
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♦TP  Matrix  (12  x  6) 
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♦TP  Matrix  (12  x  21) 


Actuator  (Sensor) 


Mode 

1 

2 

3 

4 

5 

6 

4 

.004775 

.004876 

-.04429 

0.0 

.004775 
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AD  A  1 24  702 


DECENTRALIZED  CONTROL  OF  A  LARGE  SPACE ‘ STRUCTURE  AS 
APPLIEO  TO  THE  CSDL*  2  MODE L  ( U )  AIR  FORCE  INST  OF  TECH 
WRIGHT-PATTERSON  AFB  OH  SCHOOL  OF  ENGI . .  E  S  ALDRIDGE 
UNCLASSIFIED  DEC  82  API T/GA/AA/82D- 1  F/G  22/2  NL 


MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BUREAU  OF  STANDARDS -1963- A 


oooon  non 


PROGRAM  A C0SS2 (INPUT • OUTPUT .TAPES »  TAPE6  > 


this  program  generates  an  upper  triangular  transformation 

REAL  Al(21«2l).A2(21*21)*A3(21*21), A 4(21*21) 

REAL  81(21*21 ) *62(21 *21) «B 3(21 *21) «B4 (21 *21) 

REAL  C1(21*21)*C2(21*21)»C3(21*21)*C4(21*21) 

REAL  CTCCK21  ,21)  *CTCC2(21  *2 1 )  *CTCC3  (  21 ,2 1 )  ,CTCC4  ( 21  * 21 ) 
REAL  SAT( 21,21) ,SAT2 (21 *21 ) , SAT3 ( 21 *2 1 ) • S AT 4 ( 21  * 21 ) 

REAL  AKC( 21*21) *ACT (2 1  *  21 ) *BCG (2 1 , 21 ) »KCC  (21 ,21 ) 

REAL  P(21«21)*S( 21*21) 

REAL  OA1 (21*21 )*0A2(21 *21 ) *QA3<21,21 ) ,QA4(21*21> 

REAL  0061 (21*21) *Q OB 2(21 *21)* 0063(21 *21 >*9084(21,21) 

REAL  ACG1(21,21)*ACG2<21*21)*ACG3(21«21)*AC64(21*21) 

REAL  ABG1  (21 ,21 ) « ABG2(21 .21 ) *A8G3(21 ,2 1 ) , ABG4(21 ,2 1 ) 

REAL  GAIN  1(21 *21) *GAtN2(21 , 21 > *GAI N3( 2 1  ,21 > .SAT N4 ( 21 ♦ 21 > 
REAL  KT1(21*21)*KT2 (21 *21 ) *KT3(21«21 ) «KT4(21*21) 

REAL  K0B1 (21 *21 ) « K0B2 (21 ,2 1 ) *K0B3 (21 *2 1 ) * K084 ( 21 *2 1 ) 

REAL  GAMH  A4 (21  *  21 ) *GAMMA2 ( 2 1 ,21 ) .GAMMA 3 (21  ,21 > 

REAL  T2(2 1 ,21 ) ,T 3 ( 21 , 2 1 ) ,T1 ( 21 *21) 

REAL  TRT(21 ,21 ) , TEN( 21 ,21 ) *CT (21 *21 ) * V(2 1 ,21 ) 

REAL  RK(2 1,21 >,RK4(21 ,21 ) *RK2 (21 *21 >  *RK3 (21 *21  ) 

REAL  RG2(  21*21)  ,RG3(21  *21  )  *RG1  (21  *21  ) 

REAL  MAJM(58,50)*O(17),XO(58),«(17)  .TOL.DT 

REAL  2ETA , AA( 17 ) , B8 (17 ) ,S ING(21 ) , XTR( 2 1 .21 > ,X1 (58 > 

REAL  EAT(58,58)*EAT2(58,58)« WORK (56,58),ST0R(21,21 > 

REAL  PHIA (21,21 )* PHI S(21, 21 ), MODE (2,21 ),INIT( 4,21) 
INTEGER  N*N2,NC1,NC2,NC3*NC12*NC22*NC32*NR,NR2 
INTEGER  r Cl(17) ,IC2(17) ,IC3(17) *I,J,K,L,M,KK,LL,MM 
INTEGER  0EC,Q,NACT,NS£N.rR(17) ,IER,SKIP,NC0L,NC0L1 
INTEGER  N  OA  »NGI M.NDAl »NDI Ml  * ZZ , E2 ,E3 , E4 ,Pl ,P2, P3 
COMPLEX  Z  (58)  ,U1( 12) 

COMMON /HA INA/NOA.NOAl .WORK 

COMMON/MAINB/NCOL.NCOLl 

COMMON/MA INl/NDIM* NO I Ml ,TEN,X(3364) 

COMMON/MA I N2/ST  OR 
COMMON/HA IN 3/XT R 
COMMON/SA  VE/TdOO)  *TS(  100  ) 

COMMON/ IN  OUT/KOUT , TAPE 

COMMON/NU  M/ICI ,IC2,IC3,IR , NC1,NC2,NC3  ,NR 


INITIALIZATIONS 


NOIM  =  12 
N0IM1  s  13 
NCOL  =  21 
NC0L1  =22 
NOA  =  58 
NDA1  =  59 
KOUT  =  6 
TAPE  =  9 
8  =  0 


87 


onoon  non  n n n o o n n  oo 


IER  =  0 
ZZ  =  0 


PRINTM////)* 

PRINT** •  ********  ********************************  ******* 
PRINT***  *****  ***** 
PRINT***  *****  34UPPER-SESIDUAL  ***** 
PRINT***  *****  ***** 
PRINT***  *****  B  L  C  C  K  ***** 
PRINT*. •  *****  ***** 
PRINT***  *****  C  $  D  L  II  ***** 
PRINT*.*  *****  ***** 
PRINT***  ********  *********  ********  ************  ********** 
PRINT  •(//>* 


PRINT*  *•  THIS  PROGRAM  GENERATES  A  SOLUTION** 

•  •  USING  AN  UPPER  TRIANGULAR  TRANSFORMATION  » 

PRINTM// //>• 


INITIAL  SELECTION  FOR  THREE  OR  FOUR  CONTROLLERS 


PRINT*,*  FOR  A  THREE  CONTROLLER  RUN,  ENTER  3.  OR.  • 
PRINT*,*  FOR  A  FOUR  CONTROLLER  RUN,  ENTER  4  >• 

REA0C8.O  OEC 

DEC  DEFAULT  SWITCH 

IF  (OEC *N E*4 )  OEC  =  3 
PRINT*,*  • 

PRINT*. •  THIS  IS  A  *,OEC**  CONTROLLER  RUN  • 


PHI  MATRICES  AND  CONTROLLER  ENTRIES 


PRINTM///)* 

IF  (DEC.EQ.3)  THEN 

PRINT*,*  ENTER  NC1*NC2*NC3*NR*NACT .NSEN.ZETA  >• 

ELSE 

PRINT*,*  ENTER  NCI, NC2,NC3,NC4,NACT, NSEN.ZETA  >* 
ENDIF 

REAOC8 , *  >  NCI »NC2*NC3,NR*NACT, NSEN.ZETA 
PRINT* ,NC 1 , NC2,NC 3  *NR ,NACT .NSEN.ZETA 
PRINT***  • 

PRINT*,*  ENTER  THE  *,NACT**  ELEMENTS  FOR  EACH  PHIA  • 
PRINT*,*  • 

N  s  NCI  *  NC2  ♦  NC3  *  NR 
DO  1  1=1, N 

PRINT*, 'ENTER  PHIA  •,!,»  >• 

READC6  *• >  (PHI A (I , J) , J=1 ,NACT ) 
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on  non  oooooo  nno  ooooo 


PRINT •  *,<PHIA<I«J)»J=1»NACT> 

1  CONTINUE 
PRINT • C/Z )  • 

>RINT*,*  ENTER  THE  • «NSEN« ■  ELEMENTS  FOR  EACH  PHIS 
PRINT***  • 

DO  2  1=1* N 

PRINT*, 'ENTER  PHIS  *,I,»  >* 

REA0<8**>  IPHISII *J)*J=1«NSEM> 

PRINT*,*  • ,<PHIS(I«J> ,J=1,NSEN> 

2  CONTINUE 
PRINT •(//)• 


OMEGAS 


PRINT**  *  ENTER  THE  VALUE  FOR  EACH  OMEGA  • 
PRINT*,*  • 

DO  3  I =1  *  N 

PRINT*, 'ENTER  OMEGA  *,I,»  >* 

READ( 8, *>  U<I) 

PRINT*,*  *,W«I> 

DCI)  =  -2.  •  ZETA  *  UCI) 

3  CONTINUE 


20  CONTINUE 


SECONDARY  SELECTION  FOR  THREE  OR  FOUR  CONTROLLERS,  TO 
BE  USED  FOR  RUNS  AFTER  THE  FIRST  JOB 


IF  (Q.EQ.2)  THEN 

PRINT*, •  FOR  A  THREE  CONTROLLER  RUN,  ENTER  3,  OR, 
PRINT*,*  FOR  A  FOUR  CONTROLLER  RUN,  ENTER  4  > 

REA0<8» •>  OEC 

OEC  DEFAULT  SWITCH 

IF  COEC.NE**)  OEC  =  3 
PRINT*,*  • 

PRINT*,*  THIS  IS  A  •  , DEC, *  CONTROLLER  RUN  • 
PRINT*  C //  )  ' 


IF  C0CC.EQ.3)  THEN 

PRINT*, •  ENTER  THE  VALUES  OF  NCI ,NC2 ,NC3 ,NR  >• 
ELSE 

PRINT*, •  ENTER  THE  VALUES  OF  NC1,NC2,NC3«NC4  >• 
ENOIF 

REA0C8, *)  NCI ,NC2«  MC3«NR 
PRINT*, NCl,NC2,NC3, NR 


nnnnnn  no 


PRI NT • (//  )  • 

ENOIF 

C 

C 

PRINT** •  THE  FOLLOWING  NODES  ARE  ENTERED  ACCORDING  TO  THE  • 
PRINT*.*  ORDER  IN  WHICH  THEY  ARE  ENTERED  IN  THE  DATA  FILE  • 
PRINT*. •  AND  NOT  ACCORDING  TO  THEIR  ACTUAL  NODE  NUMBER.  • 
PRINT •  (  // >  • 

PRINT*.*  ENTER  THE  *.NC1.*  CONTROLLER  1  NODES  >• 

READ (8  »  * )  (IC1 (I >  » 1  =  1 «NC1  ) 

PRINT*.*  *,<IC1CI >,I=1.NC1) 

PRINT*,*  * 

PRINT*,*  ENTER  THE  *,NC2,*  CONTROLLER  2  NODES  >• 

REA DC 5,*)  (IC2(I).I=1,NC2) 

PRINT*,*  *,(IC2(I),I  =  1 ,NC2  > 

PRINT*.*  * 

PRINT*,*  ENTER  THE  *,NC3,«  CONTROLLER  3  NODES  >• 

READ  CO,*)  CIC3CI) ,1=1, NC3) 

PRINT*,*  *,CIC3CI  )*I=1  ,NC3) 

PRINT*,*  * 

IF  CDEC.E  Q.3)  THEN 

PRINT*,*  ENTER  THE  *,NR,*  RESIDUAL  NODES  >• 

ELSE 

PRINT*,*  ENTER  THE  *,NR,*  CONTROLLER  4  NODES  >• 

ENOIF 

READCO, *)  ( IR (I ) » 1=1  ,NR  ) 

PRINT*,*  *»(IR(I),I=1,NR) 

PRINT*,*  * 


NCI  2  -  2  *  NCI 
NC22  -  2  *  NC2 
NC32  =  2  *  NC3 
N2  s  2  *  N 
NR2  =  2  *  NR 
IF  CDEC.E  Q.3)  THEN 

M  =  2  *  NCI 2  ♦  2  •  NC22  ♦  2  •  NC32  ♦  NR 2 
ELSE 

N  s  2  *  NC12  ♦  2  *  NC22  ♦  2  *  NC32  ♦  2  *  NR 2 
ENOIF 

100  CONTINUE 

PRINT*,*  TO  PRINT  ALL  OF  THE  MATRICES  ENTER  1,  ELSE  ENTER  0 

READC8.*)  Q 

PRINT*!///)* 


READ  IN  THE  WEIGHTING  MATRIX 
01  AGONAL  VALUE  FOR  EACH  MODE 


PRINT*,*  ENTER  THE  DIAGONAL  VALUES,  IN  NODE  INPUT  * 
PRINT*, •  ORDER,  FOR  THE  CONTROL  WEIGHTING  MATRIX  >• 
KEAOCG,*)  C AA  C I ) , I  si ,N) 

PRINT*, *  * 


r>  o  o  o  o  onooo 


PRINT*, CAA(I), 1=1 ,N> 

PRINT •<// » • 

PRINT*, •  ENTER  THE  DIAGONAL  VALUES,  IN  NODE  INPUT  • 

•  PRINT*,*  ORDER,  FOR  THE  OBSERVER  WEIGHTING  HATRIX  >» 
READ  Cfl,*>  (BB  (I  )  ,  1=1  ,N) 

PRINT*,*  • 

PRINT*, <BB<1),I=1,N> 

PRINT  *<///>• 


FORMING  THE  A,B,C  AND  WEIGHTING  MATRICES 


CALL  FORHAIAl ,0*W,NC1 , NC12 , IC1 ) 

CALL  FORM  B(B1 ,PHI A  »NC1 ,NC12,NACT,IC1 ) 
CALL  FORM CC Cl tPHISfNCl ,NC12,NS£N,IC1) 
CALL  FORM A(A2*D«W« NC2 «  NC2  2  « I C2  ) 

CALL  FORM  BCB2 ,PHI A,NC2,NC22,NACT,IC2) 
CALL  FORM Cf C2 , PHIS, NC2,NC22,NSEN,IC2> 
CALL  FORMA(A3,D,U,NC3,NC32,IC3) 

CALL  FORM BCB3, PHI A,NC3,NC32,NACT,IC3) 
CALL  F0RMC(C3,PHIS,NC3,NC32,NSEN,IC3) 
CALL  F0RMA(A4,0,W,NR,NR2,IR) 

CALL  F0RMB(B4 ,PHI A , NR • NR2  * NACT » IR ) 
CALL  FORM C(C4, PHIS, NR, NR2,NSEN,IR) 
CALL  FORMQtQAl ,AA,HCI ,IC1 ) 

CALL  F0RMQ(Q0B1,B6,NC1,IC1) 

CALL  FORM  Q( QA2  *  AA • NC2 , IC2  > 

CALL  FORM  Q(Q082«B8«NC2 , IC2 ) 

CALL  FORMQ(QA3,AA,NC3,IC3) 

CALL  FORMQ(QOB3,BB,NC3,IC3) 

IF  (OEC«E  Q.4 )  THEN 
CALL  FORMOC  QA4, AA , NR, IR) 

CALL  F0RMQ(Q0B4,BB,NR,IR) 

ENOIF 


PRINTING  THE  A»B,C  AND  WEIGHTING  MATRICES 


IF  (Q.EQ.l)  THEN 

PRINT*,*  THE  CONTROLLER  l  A  HATRIX  IS  • 

CALL  PRNMA1,NC12,NC12) 

PRINT*, •  THE  CONTROLLER  l  B  MATRIX  IS  • 

CALL  PRNT  <B1,NC12«NACT) 

PRINT***  THE-  CONTROLLER  1  C  MATRIX  IS  • 

CALL  PRNT  (C1,NSEN,NC12) 

PRINT*,*  THE  Cl  CONTROL  WEIGHTING  MATRIX  IS  • 
CALL  PRNT (QA1, NCI 2, NCI 2 > 

PRINT*, •  THE  Cl  OBSERVER  WEIGHTING  MATRIX  IS  * 
CALL  PRNT (QOBl «NC12«NC12) 

PRINT*,*  THE  CONTROLLER  2  A  MATRIX  IS  * 

CALL  PRNT ( A2«NC22«NC22) 
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PRINT* THE  CONTROLLER  2  B  MATRIX  IS  • 

CALL  PRNT(B2tNC22iNACT) 

PRINT*.*  THE  CONTROLLER  2  C  MATRIX  IS  • 

CALL  PR NT  (C2.NSEN.NC22) 

PRINT*.*  THE  C2  CONTROL  HEI6HTING  MATRIX  IS  • 
CALL  PRNT CQA2.NC22.NC22I 

PRINT*.*  THE  C2  OBSERVER  WEIGHTING  MATRIX  IS  • 
CALL  PRNT (QOB2.NC22.NC22) 

PRINT*. •  THE  CONTROLLER  3  A  MATRIX  IS  • 

CALL  PRNT (A3.NC32.NC32) 

PRINT*.*  THE  CONTROLLER  3  B  MATRIX  IS  • 

CALL  PRNT (B3.NC32.NACT) 

PRINT*.*  THE  CONTROLLER  3  C  MATRIX  IS  * 

CALL  PRNT  (C3.NSEN.NC32) 

PRINT*,*  THE  C3  CONTROL  WEIGHTING  MATRIX  IS  • 
CALL  PRNT (0A3.MC32.NC32) 

PRINT*, •  THE  C3  OBSERVER  UEIGHTIN6  MATRIX  IS  • 
CALL  PRNT  (Q0B3.NC32.NC32) 

C 

IF  CDEC.EQ.3I  THEN 
C 

IF  (NR.EQ.O)  THEN 

PRINT*,*  NO  RESIDUAL  TERMS  • 

GOTO  115 
ENDIF 
C 

PRINT*,*  THE  A  RESIDUAL  MATRIX  IS  • 

CALL  PRNT (A4.NR2.NR2) 

PRINT*,*  THE  B  RESIDUAL  MATRIX  IS  • 

CALL  PRNT (B4.NR2.NACT) 

PRINT*,*  THE  C  RESIDUAL  MATRIX  IS  * 

CALL  PRNT  (C4.NSEN.NR2) 

ELSE 

PRINT*,*  THE  CONTROLLER  4  A  MATRIX  IS  * 

CALL  PRNT (A4.NR2.NR2) 

PRINT*,*  THE  CONTROLLER  4  B  MATRIX  IS  * 

CALL  PRNT (B4.NR 2.NACT ) 

PRINT*,*  THE  CONTROLLER  4  C  MATRIX  IS  * 

CALL  PRNT  <C4,NSEN,NR2> 

PRINT*,*  THE  C4  CONTROL  WEIGHTING  MATRIX  IS  • 
CALL  PRNT ( QA4.NR2 ,NR2 ) 

PRINT*, •  THE  C4  OBSERVER  WEIGHTING  MATRIX  IS  • 
CALL  PRNT (Q0B4.NR2  »NR2 ) 

ENDIF 

C 

ENOIF 

C 

119  CONTINUE 


THIS  SECTION  GENERATES  THE  RICCATI  SOLUTIONS 
AMO  THE  GAIN  MATRICES  OF  EACH  CONTROLLER 
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IF  (ZZ.EQ.0)  THEN 

CALL  VMULFPIB1 *B1 * NCI 2 » NACT » NC12 *NCOl * NCO L * SA T ,NCOL * TEA ) 
END  IF 

CALL  VHULFMtCl »C1 *NSEN ,NC1 2* NC12, NCOL * NCOL *CTCCl »NCOL • IER > 
120  CONTINUE 
I  Eft  =  0 
TOL  =  0.001 

PRINT** •  THE  FOLLOWING  ARE  THE  MRIC  A*BG  l  INPUTS  • 

PRINT  •  (  //  )  * 

PRINT  *  * •  MATRIX  A1  IS  • 

CALL  PRNT ( A1 * NC 12 »NC12 ) 

PRINT***  MATRIX  SAT  <B1*B1T>  IS  * 

CALL  PRNT CSAT.NC12  .NC12) 

PRINT***  MATRIX  OA1  IS  « 

CALL  PRNT  (GA1 *NC12 «NC12 ) 

PRINT***  NCI 2  s  *  *NC12 
PRINT*!// >• 

CALL  MRIC (NC12*A1 *SAT»QA1 *S* ABG1*T0L*IER) 

IF  (22.EQ.0 )  THEN 

PRINT***  THE  RICCATI  SOLUTION  OF  AC  ♦  BCG  #1  IS  • 

PRINT*.*  IER  =  *  *  I ER 
CALL  PRNT (S*NC12»NC12) 

ENOIF 

CALL  VMULFH!B1,S,NC12,NACT,NC12,NCOL»NCOL,GAIN1 , NCOL, IER) 
IF  CZZ-EQ.1J  THEN 

CALL  VMUL  FM ! T 1 , GA INI »NACT , E4.NC12  > NCOL* NCOL* ST OR, NCOL  * IER ) 
CALL  MMUL  CRG1 *STOR *EA , E 4,NC 12,TEN I 
CALL  MMULCTl.lEN.WACT ,E4,NC12,GAIWl » 

PRINT***  THE  Gl*  6AIN  MATRIX  IS  * 

ELSE 

PRINT*,*  THE  Gl  GAIN  MATRIX  IS  • 

ENDIF 

CALL  PRNT  (GAIN1*NACT  *NC12) 

IER  =  0 
TOL  =  0.001 

CALL  TFR( ACT ,A1*NC12*NC 12*1*2) 

PRINT*,*  THE  FOLLOWING  ARE  THE  MRIC  A-KC  1  INPUTS  • 

PRINT •<//>• 

PRINT*,*  THE  MATRIX  A1  TRANSPOSE  IS  • 

CALL  PRNT (ACT *NC12*NC12) 

PRINT*,*  THE  MATRIX  CTCC1  !C1T*C1>  IS  • 

CALL  PRNT CCTCC1 »NC12»NC12) 

PRINT*,*  THE  MATRIX  OOB1  IS  • 

CALL  PRNT (QOB1 ,NC12,NC12) 

PRINT*,*  NCI 2  =  * , NCI 2 
PRINT* !// >* 

CALL  MRIC INC12*ACT,CTCC1, Q0B1,P,ACG1* TOL* IER) 

IF  (Z2.E0.0)  THEN 

PRINT*,*  THE  RICCATI  SOLUTION  OF  AC  -  KCC  Rl  IS  • 

CALL  PRNT tP,NC12,NC12) 

ENOIF 

CALL  MMUL(C1,P,NSEN,NC12,NC12*KT1) 

PRINT*, •  THE  K1  GAIN  MATRIX  IS  • 

CALL  TFR( KOB1 ,KT1 *NSEN ,NC1 2* 1 *2  ) 


CALL  PRNT ( K0S1 *NC 12 *NSEN> 

125  CONTINUE 

IF  (ZZ.EO.O)  THEN 

CALL  VMUL  FP (B2  *  B2«NC22  *NACT  tNC22*NCOL«NCOL*SAT2*NCOL*IER) 
CALL  VMULFM(C2  *C2  *  NSEN*NC22  *  NC22  * NCOL  ,NC0L,CTCC2*NC0L*IER) 
ENOIF 

140  CONTINUE 
IER  -  0 
TOL  =  0.0  01 

PRINT*  * •  THE  F0LL0UIN6  ARE  THE  URIC  A*B6  2  INPUTS  • 

PRINT •  C //  )  • 

PRINT***  THE  MATRIX  A2  IS  * 

CALL  PRNT (A2*NC22*NC22> 

PRINT***  THE  MATRIX  SAT2  (B2*B2T)  IS  * 

CALL  PRNT (SAT2*NC22*NC22) 

PRINT** •  THE  MATRIX  0A2  IS  * 

CALL  PRNT (QA2 «NC22 «NC22 ) 

PRINT***  NC22  =  *  * NC22 
PRI NT  * (//  )  * 

CALL  MRIC (NC22*A2*SAT2*QA2*S*ABG2«TOL*IER) 

IF  (ZZ.EQ.O)  THEN 

PRINT*,*  THE  RICCATI  SOLUTION  OF  AC  ♦  BCG  «2  IS  * 

PRINT*,*  IER  =  *  * IER 
CALL  PRNT (S,NC22*NC22) 

ENDIF 

CALL  VMUL FM(B2*S*NC22,NACT«NC22*NC0L» NCOL .GAIN2.NCOL, IER) 
IF  (ZZ.EQ.l)  THEN 

CALL  VHULFMCT  2*6A 1N2  *NACT  *  £3  *NC22 ,NCOL  *  NCOL*  S TOR* NCOL  *  IER  > 
CALL  MHUL  (RG2«STOR*£3*E3  «NC22*T EN ) 

CALL  MHUL  C T2 , TEN *N ACT , E3 , NC22 » GAI N2 ) 

PRINT*,*  THE  G2*  GAIN  MATRIX  IS  * 

ELSE 

PRINT***  THE  G2  GAIN  MATRIX  IS  * 

ENOIF 

CALL  PRNT  (GAIN2«NACT,NC22) 

IER  =  0 
TOL  =  0.0  01 

CALL  TFR(ACT*A2*NC22«NC22*1*2) 

PRINT***  THE  FOLLOWING  ARE  THE  MRIC  A-KC  2  INPUTS  • 

PRINT  »(//)* 

PRINT*,*  THE  MATRIX  A2  TRANSPOSE  IS  • 

CALL  PRNT (ACT «NC22*NC22) 

PRINT*.*  THE  MATRIX  CTCC2  CC2T*C2>  IS  • 

CALL  PRNT (CTCC2 ,NC22*NC22  > 

PRINT*,*  THE  MATRIX  9062  IS  * 

CALL  PRNF (0GB2 *NC22*NC22) 

PRINT*,*  NC22  =  • *NC22 
PRINT •(//>• 

CALL  MRIC  <NC22*ACT,CTCC2,QOB2«P*AC62,rOL*I&R> 

IF  (ZZ.EQ.O)  THEN 

PRINT*,*  THE  RICCATI  SOLUTION  OF  AC  -  KCC  »2  IS  • 

CALL  PRNT (P*NC22*NC22> 

ENOIF 

CALL  MHUL (C2 «P*NSEN*NC22 *NC22 *KT 2 ) 
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IF  !Z7.EQ.l)  THEN 

CALL  VMUL  FP!RK2  »GAMMA2»P3  »P3  »NS£N» NCOL »  NCOL ,STOR »NCGL  » IER ) 
CALL  MMUL ! ST0R»KT2»P3»NSEN»NC22 *KCC) 

CALL  MMUL (6AMMA2,  KCC ,  NSEN , P3 , NCZ2 tKT 2 ) 

PRINT** •  THE  K2*  GAIN  MATRIX  IS  •  ' 

ELSE 

PRINT*,*  THE  K2  GAIN  MATRIX  IS  • 

END  IF 

CALL  TFR( K082«KT2* NSEN* NC 22*1*2) 

CALL  PRNT (KQB2,NC22,NSEN) 

145  CONTINUE 

IF  1ZZ.EQ.O.OR.DEC.EQ.3)  THEN 

CALL  VMULFP !B3,B3,NC32*NACT«NC32  *NC0L*NC0L*SAT3«NC0L * IER  ) 
ENDIF 

IF  (ZZ.EQ.O)  THEN 

CALL  VMULFM!C3*C3, NSEN, NC32,NC32»NC0L,NC0L»CTCC3, NCOL ,IER) 
ENDIF 

150  CONTINUE 
IER  =  0 
TOL  =  0.0  01 

PRINT*,*  THE  FOLLOWING  ARE  THE  MRIC  A*BG  3  INPUTS  • 
PRINT*!//  )* 

PRINT*,*  THE  MATRIX  A3  IS  • 

CALL  PRNT CA3,WC32*NC32> 

PRINT*,*  THE  MATRIX  SAT3  CB3*B3T>  IS  * 

CALL  FRNT (SAT3 ,NC32*NC32) 

PRINT*,*  THE  MATRIX  QA3  IS  * 

CALL  PRNT CQA3*NC32,NC32) 

PRINT*,*  NC32  =  • ,NC32 
PRINT*!//  )  * 

CALL  MR IC !NC32 , A3 , S  A T3 ,QA  3 , S ,ABG 3 , TOL , I ER  > 

IF  CZZ.EQ.O)  THEN 

PRINT*,*  THE  RICCATI  SOLUTION  OF  AC  ♦  BCG  #3  IS  * 

PRINT*,*  IER  =  *, IER 
CALL  PRNT <S,NC32,NC32) 

ENOIF 

CALL  VMULFN!B3 , S, NC32 , NAC T , NC32 ,NC0L ,NC0L,GAIN3,NC0L , IER) 
IF  !ZZ.E0 .1 • AN0.DEC.EQ.4)  THEN 

CALL  VMULFM!T3,GAXN3,NACT,E2,NC32,NCOL,NCOL,STOR,NCOL,IER) 
CALL  MMUL!RG3,STOR,E2,E2,NC32,TEM) 

CALL  MMUL  !T3,TEN,NACT ,E2 ,NC3  2,GAI*  3) 

PRINT*,*  THE  63*  GAIN  MATRIX  IS  • 

ELSE 

PRINT*, •  THE  G3  6AIN  MATRIX  IS  * 

ENDIF 

CALL  PRNT  (GAIN3«NACT ,NC32) 

IER  =  0 
TOL  =  0.3  01 

CALL  TFR! ACT , A3 ,NC32 ,NC32  ,1,2) 

PRINT*,*  THE  FOLLOWING  ARE  THE  MRIC  A-KC  3  INPUTS  • 
PRINT*!//  )• 

PRINT*,*  THE  MATRIX  A3  TRANSPOSE  IS  • 

CALL  PRNT !ACT,NC32,NC32) 

PRINT***  THE  MATRIX  CTCC3  !C3T*C3)  IS  • 
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CALL  PRNT <CTCC3*NC32*NC32) 

PRINT***  THE  MATRIX  G083  IS  * 

CALL  PRNT < Q0B3 *NC 32* NC32) 

PRINT***  NC32  =  *  *NC32 
PRINT*!// )• 

CALL  HR IC ( NC32 »ACT *CTCC3 *  Q0B3*P  * ACG3 • TOL* IER 1 
IF  <2Z.EQ.O)  THEN 

PRINT***  THE  RICCATI  SOLUTION  OF  AC  -  KCC  R3  IS  • 

CALL  PRNT  <P<NC32<NC32 1 
ENOIF 

CALL  HHUL  <C3*P *NSEN* NC32* NC32 »KT3  ) 

IF  <ZZ.E«*1)  THEN 

CALL  VHULFP<fiK3«GAMHA3.P2*P2*NSEN»MC3L»'iC3L*STCR»r;COL.IEP  1 
CALL  HHUL  <ST Oft  *K7  2  ,P2  » 'iSEFi  ,NC32  ,KCC> 

CALL  HHUL  CGAHHA3 *KCC« '.SEN  * P2 »NC32  »KT *  ) 

PRINT***  THE  K3*  GAIN  MATRIX  IS  * 

ELSE 

PRINT***  THE  K3  GAIN  MATRIX  IS  « 

ENOIF 

CALL  TFR( K0B3  *KT3  *K3Eft  «NC32  *1 *2  ) 

CALL  PRNT  <K0B3»NC32»NSEN) 

155  CONTINUE 

IF  <DEC.EQ.4>  THEN 

CALL  VMULFP<B4,B4,r;R2*NACT*NR2*NCC*l  »NC3L*SAT4,NCCL»IE~ 1 
IF  <ZZ.E3.01  THEN 

CALL  VHULFM<C4tC4,N3EN*NR2*NS2*NC0L*NC0L*CTCC4*NC0L*IER> 
E^DIF 

160  CONTINUE 
IER  =  0 
T  CL  =  C.0  01 

CALL  HR IC  <NR2*A4*SAT4*QA4*S*  ABG4  *TOL  *1 E3  I 
IF  <ZZ.EQ.O>  ’HEN 

PRINT***  THE  RICCATI  SOLUTION  OF  AC  ♦  BCG  »♦  IS  • 

PRINT***  IER  =  **IEP 
CALL  PRNT  <S«NR2 ,NP2 1 
Ef.DIF 

CALL  VHULFM<B4«S*NR2*NACT  *NR2*NC0L*NC0L*GAIN4*NC3L  »IER ) 
PRINT***  THE  64  GAIN  MATRIX  IS  * 

CALL  PRNT  <6AIN4*NACT  *NR2 1 
IER  =  0 
’ OL  =  0.0  01 

CALL  TFR<ACT*A4*NR2*N’2»1*2> 

CALL  HRIC  <NR2  *  ACT  «CTCC4 *Q0B4 ,P ,ACG4 *T0L *IER > 

CALL  HHUL  <C4,P*NSE»!*NR2*NR2  *KT4  I 
IF  <Z7.EQ .1 )  THEN 

CALL  VHULFP<RK4*6AHMA4*P1«P1,NSEN<NCOL*NCOL*STOR*NCOL<IER1 
CALL  HHUL  <STCR*KT4*P1  *NSE»J«ftR?  *KCC> 

CALL  MMUL  <GAMMA4*KCC*N$EN  *P1 *WR2  *KT4  > 

PRINT***  ’HE  K4*  GAIN  MATRIX  IS  * 

ELSE 

PRINT***  THE  K4  GAIN  MATRIX  IS  • 

Ef  OIF 

CALL  TFR  <  K084  «KT4  *  NSEi.  2«!  .1 
CALL  PRNT  (KC84.NR2  »N  :  Ef' , 
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165  CONTINUE 
Er  OIF 

THIS  SECTION  GENERATES  THE  BLOCK  SEGMEN'S 
OF  MAJM  AND  PUTS  THEM  INTO  THE  MAJM  MATRIX 

THE  THRE£  CONTROLLER  MATRIX  WILL  CONTAIN 
RESIDUAL  TERMS  C SEE  DIAGRAM  BELOW). 

THE  FOUR  CONTROLLER  MATRIX  DOES  NOT 
CONTAIN  RESIDUALS  CYET  )• 

THE  THREE  CONTROLLER  MAJM  MATRIX  WITH 
RESIDUALS  WILL  LOCK  LIKE: 

****  *** 


A1*BG1 

B1G1 

B1G2 

B1G2 

B1G3 

31G3 

0 

0 

A1-KC1 

K1C2 

0 

K1C3 

0 

K1C° 

B2G1 

B2G1 

A2+BG2 

82  G2 

B2G3 

B2G3 

0 

K2C1 

0 

0 

A2-KC2 

K2C3 

0 

K2CB 

B3G1 

B3G1 

B362 

B3G2 

A3*BG3 

B3G3 

3 

K3C1 

0 

K3C2 

0 

0 

A3-KC3 

K3CR 

BRG1 

BR61 

URG2 

BRG2 

BRG3 

3R63 

AR 

•  *  *  » 


K  *  2  *  N  Cl  2 

KK  =  K  ♦  \C22 

L  =  2  *  NC22  ♦  K 

LL  =  L  ♦  NC32 

PI  s  2  *  NC32  ♦  L 

IF  CDEC.EQ.3)  THEN 

MM  =  2 * NC 1 2  ♦  2*NC22  ♦  2*NC32  ♦  NR2 

ELSE 

MM  =  2*NC 12  ♦  2 *NC22  ♦  2*NC32  ♦  2*NR2 
P2  =  Pi  ♦  NR2 
ELDIF 
C 

00  200  1=1* MM 
DO  200  J=  1  *  MM 

200  MAJMC  I  *  J)  =  0.0 
DO  201  I=1«NC12 
DO  2C1  J= 1 *NC12 

201  MAJMC I  * J)  =  ABGl(IfJ) 
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DO  202  1=1 #\C22 
DO  202  J=lf:;C22 

202  MAdMCI*Kt  J+K)  =  ABG2CItd> 

DC  20  3  1=1  *»:C32 

DC  203  J= 1 *NC32 

203  HAdMC I*L* d+L  >  =  A8G3(!fJ> 

CALL  TFRC  AKCtACGl  t  MCI  2  »f.Cl  2  tl  1 2  ) 

DC  204  I  =  1 » *.’ C 1 2 
DC  2C4  J=1,NC12 

204  HAdMC I *NC 12  td+NCl 2 )  =  AKCCItd) 

CALL  TFR<  AKCtACG2t'»C22tMC22tlt2  ) 

DO  205  1=1 «NC22 

DO  205  J=l»NC22 

205  HAdMC I+KK « d+KK)  =  AKCCItd) 

CALL  TFR<  AKC»ACG3»\C32*r-C32tl*2) 

DC  206  I=1»NC32 
DC  206  d=lt*-C32 

206  MAdMCI+LL td*LL>  =  AKCCItd) 

CALL  MMLJL  <B1  tGAI  Nil  *NC12  *NACT  »NC1 2  tBCS  ) 
DO  207  1=1 t  NC12 
DC  207  0= 1 »NC12 

207  MAJMtl* J>.‘ C12)  =  BCGCItd) 

CALL  MMUL (B1,GAI\2 »NC12  »NACT»NC22»BCG1 
DO  238  1=  1  *f.'C12 
DO  2C6  J=  1  tNC22 
M  AdMC  I  »  J+-  K  )  =  8CG(  ltd) 

208  MAdMCI,d«-KK)  =  BCG(IiJ) 

CALL  MMUL  <BltGAI\3*NC12»NACTfNC32*BCG> 
DO  23?  1= 1 «NC12 
DO  20?  d=  1 1 NC32 
HAJMU  td*c>  =  8CGCItd> 

209  MAdMCItd«-LL)  =  BCGCItd) 

CALL  P«ULCB2*GAIM»NC22»NACTfNC12t8CG) 
DO  210  1= 1 *NC22 
DO  21C  J=  1  *’.C12 
HAdMC I*K»  J)  =  BCG ( I  * J ) 

210  MAdMCI*Ktd*fJCl2)  =  BCGCItd) 

CALL  HMUL ( B2 1 GAIN2 »NC22 tN ACT t NC22 t BCG  ) 
DC  211  1= 1 1 NC22 
DO  211  •;=lt.NC22 

211  HAdHCI+Kt  d>KK )  =  BCGCItd) 

CALL  MMl/L  <B2tGA:f:3tNC22tNACT*NC32tBC6  ) 
DO  212  1= 1 «NC22 
DO  212  d= 1 t NC32 
HAdMC I*Kt d*L )  =  BCGCItd) 

212  HAdMC I ♦Kt  d>LL )  =  BCGCItd) 

CALL  MHUL  CBJtGAI'l  1  tfJC  32  tNACT  tNCl  2  tBCS  ) 
DO  213  I=1*MC32 
00  213  d=  1  t'lC12 
MAJHC I *L«  d)  =  BCGCItd) 

213  HAdMC  Z+Ltd+NC12)  =  BCGCItd) 

CALL  VMUL  CB3  tGAI*j2  tNC32  t’JACT  tMC22t3CG) 
DC  214  1= 1 tNC32  ” • ’ 

00  214  J= 1«NC22 
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PAJMC  I^Lf  J-»K)  =  BCG  (I  «  J) 

214  PAJHC*Lt  J*KK>  =  RCG<  I*  J) 

CALL  f*HJL  CB3f GAIN2  «K:C32  tNACT  *NC32  »3CS  I 
OC  215  I=1»nC32 
DO  215  J=  1  »‘JC 32 

215  MAJMC^L*  J*LL>  =  BCG(!#J) 

CALL  WMUL  < KOB1  «C2  * NCI  2  *N3EN  ,f.C22  * KCC) 
DO  216  1= 1 » MC12 
00  216  J=1»NC22 

216  HAJKCI+NC12 tJ>K)  =  KCCCIiJ) 

CALL  PMUL  <K0Bl,C3*\C12».NSEN,NC32tKCC> 
DO  217  1=  1  *«\C1 2 
OC  217  J=  lfNC32 

217  MAJMC Z+'\C  12  »U*L  )  =  KCCCfJ) 

CALL  HMUL  (KCB2tCl  fNC2 2  f NSEN f NC12 f KCC > 
DO  215  I=lf\C22 
DC  21a  J=l.f,C12 

218  HAJKC  +  KKfJ)  =  KC  C  ( I  »  J  ) 

CALL  MMLIL  (K082*C3«NC22tNSEN»^C32*KCC) 
DO  215  I~ 1 «NC22 
00  21c  J=lt‘iC32 
215  PAJMC^KKfJ+L)  =  KCCCfJ) 

CALL  HPCL  (K083,C1  »uC32  fNSEKffcC12  ,KCC  ) 
DO  22C  Isl»NC32 
DO  22 C  J=lfNC12 

220  PAJMC*LL  »J)  =  KCC(IfJ) 

CALL  HMUL  (KGB3, C2  *hC32  fNSENff<C22  »KCC  ) 
OC  221  I=1.NC32 
DO  221  J=1«NC22 

221  MAJHC  +  LL  f  J*K>  =  KCCCfJ) 

CALL  MULL  <B4,GAI\1,.\S2,NACT,NC12*BCG> 
DC  222  I=ltVc2 
DO  222  J=lfNC12 
M  A  JM  ( I ♦ FI »  J )  =  BCG (I  *  J) 

222  HAUM<I*-P1  »JV!C12)  =  BCGCfd) 

CALL  VMUL  CB4  tGAIN  2  f Ka  2  f  NACT  »NC22»BCG> 
DC  223  I=ltN?2 
DO  223  J=ltNC22 
HAJPC+Pl  tJ»K>  =  8CGCfd> 

223  MAJK<I«-P1  f  J-*KK)  =  BCGCfJ) 

CALL  VMUL  (B4,GAI")3  ♦N°2tNACTfNC32»BCG) 
DC-  224  I=1*NP2 
DC  224  d=l*f;C32 
MAJMCI+P1  tJ+L >  =  BCG (  I  « «J ) 

224  MAJMC+P1 fJ+LL)  =  BCG (I »J) 

CALL  MMUL(K0Bl«C4*r.C12*NSENtNR2fKCC) 

OC  225  l-  1  **.C12 
OC  225  d=l*NR2 

225  HAJMC  I*).C  12t J+Pl )  =  KCCCfJ) 

CALL  MMUL  CK0B2*C4,NC22tfiSEN»N«2*KCC) 

DO  226  1=1 t  NC22 
00  226  J= If N9 2  • 

226  MAJP<  i*KK  f  J*cl  )  =  KCCCfJ) 

CALL  MHUL  < K0B3 1 C4  »  NC32  f  .NS EN  t  NR 2  « KCC  > 


on  no 


DO  2 27  I=1*NC32 
DO  227  J=ltNP2 
227  MAJM<I*Li- *J-*P1  > 


KCC  < I » J  > 


IF  (DEC »E  Q*3  )  THEN 


CALL  F0P*1ACA4fD»W*NR#\R2*IR) 

DO  229  I=lfNR2 

DC  229  J=1»NP2 

229  KAJMCI  +  P1  td+Pll  =  AM  I  » J) 

C 

ELSE 

C 

DO  23C  I=l*f;'2 

DC  230  J=1*NR2 

230  KAJM<I*P1  »J-*P1)  =  A8G4<I*JJ 
CALL  7F9C  AKC« ACG4 ,NR2 *NR2 »1 f2 > 

OC  231  1= 1 *NS2 

DC  231  J= 1  *  NR  2 

231  MAJMC  I*P2  t  J+P2)  =  AKC(I«J) 

CALL  HMUL (B1 «GAI^4 ,NC12 ,NACT »NR2 *BCG) 

OC  232  I=1»NC12 

DO  232  JsltNR2 

MAJMC I »  J* PI )  =  8CGCI *  JJ 

232  MAUMI I  »«!♦  32 1  =  BCSCr.JJ 

CALL  !H8UL  <B2»GAI\4«[\C22  fM  ACT  *f,R  2  t  BCD  ) 
DO  233  1= 1 *mC22 
DO  233  d=l»NP2 
HAJHd+K*  J+Pl >  =  BCG< I » J> 

233  *AJM<I*Kt J+P2)  =  BCG<1»J) 

CALL  KMJL  CB3tGAI\4»NC32*NACT*NR2,RCG> 
DO  234  I=1#NC32 
OC  234  J=l,NP2 
*AJM<I«-L*  J+Pl  >  s  BCGC*J) 

234  PAJM(I>L«  J*P2)'  =  BCGfltJ) 

CALL  PHUL  <B4»GAI.\4  »NR2t NACT *NR2  »BC6) 
OC  23b  I=1»NP2 
DO  239  J= 1 » NP  2 

238  MAJPCI^PI »J*C21  =  BCG(I,J> 

CALL  MHUL <K0B4, Cl i NR2 «NSEN tNC12 * KCC ) 
DC  242  I=1,NR2 
DC  242  J= 1 «NC12 

242  MAJM<I»P2 »J>  *  KCCCItd) 

CALL  MMULCK0B4fC2»NR2*NSEN*NC22*KCC> 
OC  243  I* 1 » NR 2 
OC  243  d=l*NC22 

243  MAJM«I»P2 #J*K)  s  KCCCI»J) 

CALL  KMUL  C  K0B4  »C3  » ;.R2  ♦  NSEN  iNC32  *KCC  ) 
00  244  1= 1 *NR2 
OC  244  dslf*JC32 

244  MAJM<I»P2  td+L)  =  KCC(ItJ) 

Ef.DIF 
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ooooo  onooo 


TA»  OAH  C 


N  CU  UC  HAVE  THE  MAdM  MATRIX 


rp-- 


IF  <DEC.£Q*4>  ’HEN 

PRIM*,*  THE  FOUR  COMFOLLER  MAJM  IS  * 

ELSE 

PRINT* ,  •  THE  THREE  CONTROLLER  MAJM  U/RESIDUALS  13  * 
ENDIF 

CALL  PRNT XL<MAJM,wm,h«) 


EIGENVALUE  ANALYSIS  SEC'ION 


PRIM •<///)• 

PRINT*, *  OVERALL  SYSTEM  EIGENVALUES  • 

CALL  EIGTFCMAJM,MM,\OA,0,Z,TEN,NCCL*VOR<,IER> 

PRINT*,*  IER  =  *,IER 
DG  ACC  1=1, MM 

♦00  PRINT*,*  *»ZCI> 

PRINT • <// ) * 

C 

PRINT*,*  EIGENVALUES  OF  AC  ♦  SCG  SYSTEM  1  • 

CALL  EIGRF<AB6l»\C12»NC0L,Q,Ml ,TEN»NCCL»STOR,IE; > 
PRINT*,*  IER  =  •• IER 
00  4C1  1= 1  ,KC12 

♦  01  PRINT*, •  **U1<II 

PRINT*  (//  )  * 

C 

PRINT*,*  EIGENVALUES  OF  AC  -  KCC  SYSTEM  1  * 

CALL  TFRC  AKC, ACG1 , NC12 »*C12 ,1 ,2 ) 

CALL  EIGRF(AKC,NC12*?4COL»0,Ml*TENtNC0Lt3ToR*IER  » 
PRINT*,*  IER  =  *  * IER 
OG  AC 2  1= 1 ,NC12 

♦02  PRINT*,*  *,UI(I> 

PRI NT*  f // > • 

C 

PRINT*,*  EIGENVALUES  CF  AC  ♦  BCG  SYSTE*  2  » 

CALL  EZ63  F(  ABG2  ,NC22,i\COL,0,Wl,TEP,N'COL,STOR,IEP) 
PFINT*,*  IER  =  *,I£R 
DO  AO 3  1  =  1  ,NC22 

♦  03  PPIMT*,*  *,W1U> 

PRINTM//  >* 

C 

PRINT*,*  EIGENVALUES  OF  AC  -  KCC  SYSTEM  2  • 

CALL  TFR(  AKC,  AC62  *NC22  ,'.C22,1,2) 

CALL  EIGR  F  <  AKC  ,NC22,.^C3L,  <3  ,  U1  ,  TEN,  N  COL,  5TCR*IER> 
PPINT*,*  =  »,IER 

00  AO ♦  1=1  ,NC22 

AOA  PRINT*,*  *,W1CI» 

PRINT* (//  )  • 

C 
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PRINT* * •  EIGENVALUES  DF  AC  ♦  3CG  SYSTEM  3  •  " 

CALL  E 163  F  ( ABG3  *\C  32»\CCL,0»Ul»TEN*NCCL»STOR»IER> 

PRINT*,*  IER  =  •  »  IER 
OC  405  I=1,NC32 

405  PRINT*,*  *,W1II> 

PRINT**// >• 

C 

PRINT*,*  EIGENVALUES  OF  AC  -  KCC  SYSTEM  3  • 

CALL  TFRC  AKC  ,  ACG3  , *!C3  2  ,f.C32 , 1 ,2  > 

CALL  EIG°.  F*AKC»NC32»NCOL»C,yl,TEN»NCOL,STQR,IER  I 
PRINT*,*  IER  =  *  *  I ER 
OC  406  I=1,NC32 

406  PRINT*, •  *,UltI) 

PRINT  •<//)• 

C 

IF  <NR*EO*0>  THE\ 

F’PINT* ,  *  NO  RESIDUAL  TERM  EIGENVALUES  * 

GCTO  410 
ENDIF 
C 

IF  (DEC.EQ.4)  THEN 
C 

PRINT*,*  EIGENVALUES  OF  AC  *  8CG  SYSTEM  4  * 

CALL  EIG?FCABG4,N'2*\COL*O»M1,TEN*NC0L*3TOR,IER> 

PRINT*,*  IER  =  *, I ER 
DO  407  1=1, NP2 

407  PRINT*, *  *»U1<II 

PRINT  *  * // >• 

C 

PRINT*,*  EIGENVALUES  OF  AC  -  KCC  SYSTEM  4  • 

CALL  TFR<AKC,ACG4,NF2,NR2,1,2) 

CALL  EIG^F*AKC,NR2,NCCL,0,U1,TEN,NCCL,STCR  ,IER  > 

PRINT*,*  IER  =  * , I ER 
00  406  I=1,N*2 

403  PRINT*,*  *,W1 (I) 

C 

ELSE 

C 

PRINT*,*  EIGENVALUES  OF  THE  A  RESIDUAL  MATRIX  * 

CALL  EIGRFCA4,NR2,NCGL,0,W1,TEN,NC3L,ST3R,IER> 

PRINT*,*  IER  =  • , IER 
00  40S  1  =  1, N«2 

409  PRINT*,*  *,W1*I) 

Ef  OIF 

PRINT**///)* 

C 

410  CONTINUE 

IF  (ZZ*EQ  *1 )  GOTO  600 
C 

c 

c  THIS  SEC” ION  FORMS  THE  TRANSFORMATION  MATRICES* 

C  TO  GET  MAJM  IN  UPPER  TRIANGULAR  FORM,  IT  IS 
C  NECESSARY  TO  DRIVE  THE  B2G1,  83G1,  B3G2,  K2C1, 

C  K3C1  AND  K3C 2  TERMS  TO  ZERO  (THREE  CTLRS). 
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nnnnnonnnnnnnnnnnnoononnnnnnnonnonnnonoon 


UHEN  FOUR  CONTROLLERS  ARE  USED.  *  HEN  MAJM 
MILL  INCLUDE  G4  AND  K4  TERMS.  IT  WILL  THEN 
BE  f'ECESS APT  TO  DRIVE  THE  B4G(I  )♦  A>  D  K4C(I> 
TERMS  TC  ZERO*  AS  WELL. 


AFTER  THE  TRANSFORMATION  IS  COMPLETE* 

THE  THREE  CONTROLLER  MAJM  (WITH  RESIDUALS) 
WILL  LOCK  LIKE: 


***•  *** 


A1*BG1 

B1 61 

B1G2 

B1G2 

B1G3 

B1G3 

0 

0 

A1-KC1 

K1C2 

0 

K1C3 

0 

K1CR 

0 

0 

A2*BG2 

B2G2 

B2G3 

32G3 

0 

0 

0 

0 

A2-KC2 

K2C3 

0 

K2CR 

0 

0 

0 

0 

A3+8G3 

B3G3 

C 

0 

0 

0 

0 

3 

A3-KC3 

K3CR 

BRG1 

BRG1 

BRG2 

BR62 

BR  63 

BRG3 

A° 

*  .  * 


WHERE  THE  MON-ZERO  TERMS  INCLUDE  THE  TRANSFORMATION 
MATRICES. 


*•* 


ON  WITH  THE  TRANSFORMATION  MATRICESC 

FIRST  THE  OBSERVER  GAIN  MATRIX*  K4 
UHEN  USING  FOUR  CONTROLLERS 


IF  (DEC .E 0.4)  THEN 

CALL  TFR( CT*C1*NSEN*NC12  *1 *2  > 

OC  500  1= i *MC1 

DO  500  J=1«NSEN 

500  V(IfJ)  =  CT(I,J) 

CALL  TFR( CT«C2  *NS  EN*NC22  *1 *2  ) 
00  501  1=  1  *NC2 
00  501  J: 1 *NSEN 

501  V(I+NC1« J )  =  CT(ltJ) 

CALL  TFR(CT,C3,NSEN*NC32*1*2  ) 
00  502  I* 1 *NC3 
OC  502  J* 1 *NS£N 

502  VII ♦MC2  +N Cl  *  J )  =  CT( I «  J ) 
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on  o  onooooo 


NPV  =  NCI  ♦  NC2  ♦  NC3 
PRINT  *  »■  V  (C1/C2/C3)  IS  » 

CALL  PRNT  CV,NRV,NSES> 

CALL  LSV3FCV,NC0L»NrV,N3EN,TEN,NCCL,-l,SIN6»ST0R,IER) 
PRINT  *»  •  • 

PRINT*.*  LSVDF  K4  IER  =  *,IER 
PRINT*  (// »• 

PRINT*,*  V  OUT  OF  LSVOF  IS  * 

CALL  PRNT  ( V* NSEN * NSEN ) 

PI  s  NSEN  -  NR V 
IF  (Pl.LT.l)  THEN 
DC  503  1=1,  NSEN 

503  GAMMA  4(1. 1 )  =  VCI.NSEN) 

PI  =  1 

ELSE 

DC  504  1=1, NSEN 
00  504  J=1,P1 

504  GAMMA4CI«J)  =  VCI.U+NRV) 

Ef.OIF 

PRINT***  TRANSFORMATION  MATRIX  6AMMA4  • 

CALL  PRNT  (GAMMA4 . NSEN, PI ) 

CHECK  TC  SEE  THAT  GAMMA4  IS  ORTHOGONAL  TO  01. C2.  AND  C3 

NOTE!  AKC  I-l  THIS  SECTION  IS  JUST  A  WORK  AREA  TO  TEST 

THE  ORTHOGONALITY  CF  CT  *  GAM«A.  IN  ALL  CASES  IT 
SHOULD  BE  A  BLOCK  ZERO  MATRIX . 

CALL  TF?(CT»C1  , N ' EN . *<C1 2 . 1  * 2  ) 

CALL  MMUL  (CT,GAMMA4,NC12,NSEN,P1 .AKC) 

PRINT*.*  C1T  *  GA y-H*  * 

CALL  PRNr  (  AKC  .NCI  2  .  PI  ) 

CALL  TFF ( CT  »C2  »  No  EN.NC22.1 *2 ) 

CALL  MMUL  CCTtGAMMA4,..C22,NSEN»Pl  .AKC) 

PRINT*. •  C2T  *  GAMMA4  • 

CALL  PRN~  <AKC.NC22.PD 

CALL  TFRC  CT ,C3«N2  EN.NC32 .1.2) 

CALL  MMUL  CCT.GAMMA4.NC32.NSEN. PI .AKC) 

PRINT*.*  C3T  *  GAMMA4  • 

CALL  PRNT  <AKC.NC32.PD 

PRINT*,*  CI23  SINGULAR  VALUES  * 

CALL  PRNT (SING. NR  V.l ) 

CALL  TFRC  TRT »GANM A4.NSEN.P1 .1 . 2 ) 

CALL  MMUL (TRT.GAMMA4.FI.NSEN. PI. RK) 

CALL  GMINVCP1.P1.RK.RK4  .J.TAPE) 

CALL  TFRC  CT .C4.NSEN.NP2  .1 .2  ) 

CALL  MMUL CTRT.C4.P1.NSEN.NP2, AKC) 

CALL  MMUL CCT,GAMMA4,\P2,NSEN,P1,K0B4) 

CALL  MMUL  CK0B4,PK4,NR2,P1 ,P1 »STOR> 

CALL  MMUL  CSTOR ,AKC ,NR2 ,P1 ,NR2,CTCC4 ) 

ENDIF 

THIS  CTCC4  WILL  BE  SUBSTITUTED  BACK  INTO  M?IC 
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n  n  o  oooono 


SYSTEM  4  TO  6ET  A  NEW  K4. 

NOW  THE  OBSERVER  GAIK  MATRIX  ,  K3 


CALL  TFR(CT,Ci,N:EN*f.C12*l,2  ) 

00  505  1=1, *C1 
00  505  JsltNSEN 

505  V(ItJ)  =  C*(I,J) 

CALL  TFR(CT,C2,N3EN,f;C22,l,2> 

DO  506  I - 1 1 NC 2 
00  506  J=1»N3EN 

506  V( I+NCl « J  )  =  CT  (I  ,  J) 

NRV  =  NCI  ♦  \C2 
PRINT*, •  V  <C1/C2>  IS  • 

CALL  PRNT  (VtNRVf  r,SEN> 

CALL  LSVOF* V* \COL *  NR V, NSEN »TEN,NCOL *-l *  SI NG» STOR , I E= ) 
P6INT*,»  • 

PRINT*,*  LSVOF  K3  IEP.  =  *,IER 
PRINT**//  |  • 

PRINT*,*  V  OUT  OF  LSVOF  IS  * 

CALL  PR  NT  <V,NSEf.,NS£M 
P2  =  NSEN  -  NTV 
IF  (P2.LT.1)  THEN 
DC  SC  7  ZsttNSEN 

507  G  AMMA3  (  1,1)  =  VCI,NSENI 
P2  =  1 
ELSE 

DO  506  I=1,NSEN 
DO  506  J=1,P2 

509  GAMMA  3(1*  J)  =  V(I,J*N?V> 

ENOIF 

PRINT*,*  TRANSFORMATION  MATRIX  GAMMA3  * 

CALL  PRNT  (GAMMA3 *  NSEN«P2 ) 

CHECK  TC  SEE  THAT  6AMMA3  IS  ORTHOGONAL  TO  Cl  AND  C2 

CALL  TFR ( CT *C1 , NS EN, NC12 , 1 ,2  ) 

CALL  MMUL (CT,GAMMA3, NCI2 ,N SEN, P2 ,AKC ) 

PPINT * , •  CIT  *  GAMMA3  • 

CALL  PRNT (AKC,NC12,P2) 

CALL  TFR«CT,C2,NSEN*<\C22,1,2> 

CALL  MMUL (CT , GAMMA3, <C22,NSEN,P2,AKC) 

PRINT*,*  C2T  *  GAHMA3  • 

CALL  FRNT (AKC,NC22 ,P2 ) 

C 

PRINT*,*  C12  SINGULAR  VALUES  • 

CALL  PRNT (SING«NRV«1) 

C 

CALL  TFR(TRT,GAMMA3,NSEN,P2,1,2) 

CALL  MMUL (TRT  «GAMMA3,P2 ,NSEN,P2,BK) 

CALL  GMINV(P2,P2,RK,FK3,J,TAPE> 

CALL  TFR(CT«C3,NSEN,NC32,1,2) 
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non  nnnnnonnn 


t- 


CALL  NMUL  CTRT  »C3 , P2 *f  SEN,  NC32  ,  AKC > 

CALL  NMUL(CT,GAMVA3,NC32,NSEN,P2,K3B3> 
CALL  MMUL  CK0B3,PK3,.NC32»P2,P2,  STOP) 
CALL  KMUL  CSTOR , AKC ,NC32 «P2,NC32,CTCC3 > 


CTCC3  WILL  BE  SUBSTITUTED  BACK  INTO 
NRIC-STSTEM  3  FOR  A  .\EW  K3. 


NOW  THE  OBSERVER  GAIN  MATRIX «  K2 


CALL  TFRC CT,C1 *NSEN,NC12,1,2 > 

DO  509  Is it NCI 
DO  509  J=  1  iNSEtN 

509  VCI.J)  s  CTCI t J) 

PRINT***  V  CC1>  IS  • 

CALL  PR  NT  CV, NCI, NSEN) 

CALL  LSVDFC  V*  NCOL  «  NCI  *  NSEN»TEN,NCOL  *“l»SING*STOR»IEp > 
PRINT *»•  * 

PRINT*,*  LSVOF  K2  IER  s  *.IEF 
PRINT •{// ) • 

PRINT*,*  V  OUT  OF  LSVOF  IS  * 

CALL  PR NT  CV,NSEN,NSEN) 

P3  =  NSEN  -  NCI 
IF  CP3.LT. I)  THEN 
DO  510  1=1, NSEN 

510  GAMMA2  C I, 1 1  =  VCI.NSEN) 

P3  =  1 

ELSE 

DO  511  1=1, NSEN 
DO  511  «J=  1  ,P3 

511  GAMMA2  C I, J)  =  VC  I » J*NC1 ) 

END  IF 

PRINT*,*  TRANSFORMATION  MATRIX  6AMHA2  • 

CALL  PR NT  C6AMMA2  »NSE.N»P3) 

CHECK  TO  SEE  THAT  GAMMA2  IS  ORTHOGONAL  TO  Cl 

CALL  MHUL CCT,GAMMA2,NC12,NSEN,P3,AKC» 

PRINT*,*  Cl T  *  GAMMA2  • 

CALL  PRNT CAKC,NC12,P3> 

PRINT*,*  Cl  SINGULAR  VALUES  • 

CALL  PRNT CSING, NCI, 1) 

CALL  TFR  C  TRT*GANMA2»  NS  EN,P3 ,1 *2 ) 

CALL  HMUL CTRT, GAMMA2.P3, NSEN, P3,RK> 

CALL  GMIMVCP3*P3,RK,RK2,J*7APE) 

CALL  TFRCCT.C2, NSEN, NC22, 1,2) 

CALL  MHUL  CTRT ,C 2»P3»N'EN ,NC22 , AKC I 
CALL  NMUL  CCT  ,GAMMA2  ,'.C22  ,N  SE  N*  P3  »K3B>  ) 

CALL  PMUL CKOB2,RK2,NC22,P3,P3,STOR) 

CALL  MMUL  CSTGfi , AKC ,NC22  ,P3,NC22,CTCC2> 
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CTCC2  UILL  BE  SUBSTITUTED  BACK  IN’O 
MRIC-SYSTEM  2  PC'*  A  NEU  K2 


NOW  THE  CONTROLLER  6A1N  MATRIX  *  G3 


IF  (DEC.EQ.4)  THEN 
DO  512  1=1, NR 
DC  512  J=1*NACT 

512  VCI,J>  =  B4 < I ♦NP *  J) 

PRINT  *  *  *  V  <80  I  £  • 

CALL  PRNT  <V*\R,NACT> 

CALL  LSVDF<V,NCOL,NR,f.ACT  »TEN»NC0L*-1*SINS*STCR*IER  ) 
PRINT***  • 

PRINT***  LSVDF  62  IER  =  **IER 
PPINT*<// >• 

PRINT *••  V  OUT  CF  LSVDF  IS  • 

CALL  P.RNT  <V*NACT,NACT> 

E2  =  f;ACT  -  NR 
IF  (E2.LT.1)  THEN 
DO  513  I=1.NACT 

513  T3 ( I « 1 )  =  V< I »NACr ) 

£2  =  1 
ELSE 

DO  51  *  1  = 1 ,NACT 
DO  514  J= I *E2 

514  T3( I* J)  =  V < I  * d* N -  ) 

ENOIF 

PRINT***  TRANSFCR  RATI  ON  MATRIX  T3  • 

CALL  PRNT  <T3*NACr*E2> 

CHECK  TO  SEE  THAT  T3  IS  0RTH060NAL  TO  B4 

NOTE!  IN  THIS  SEC T I CN *  BCG  IS  THE  WORK  AREA 
FOR  B  *  T.  IN  ALL  CASES  THESE  SHOULD 
BE  BLOCK  ZERO  MATRICES. 


CALL  MMUL  <B4  *  T  3*NP  2  »NACT *E2«BCG t 
PRINT*,*  B4  •  T3  • 

CALL  PRNT <BCG*NR2,E2> 

C 

PRINT***  B4  SINGULAR  VALUES  • 

CALL  PRNT (SING, NR, 1) 

C 

CALL  VMULFM<  T3  *T3  «NAC  r  «E2  *  E2 ,NCOL*NCOL*RK,NCOL  *IE3 > 

CALL  6MIN V<E2*E2,PK*nG3,d,TAPE) 

CALL  MMUL  <B3*T3,NC32,’.ACT  ,E2,K083> 

CALL  MMUL <K0B3 ,P53 ,NC32 ,E2 »E2 * SAT3 > 

CALL  VMULFP<SAT3, ’3,%C32*E2*NACT,\C0L,NC0L*K0B3*MC0L*:EP> 
CALL  VMUL FP( KC83*B 3 *NC32 ,NACT*NC 32 »NC2L*NC0L *S AT3*NC0L *IE3 ) 
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ENDIF 


THIS  SAT 3  MILL  BE  SUBSTITUTED  BACK  INTO  MRIC 
SYSTEM  3  FOR  A  NEW  G3. 


NOW  THE  CONTROLLER  GAIN  MATRIX*  G2 


DC  515  1=1, NC3 
OC  515  J=1.NACT 

515  V(I,J>  =  B3(I«-NC3,J> 

IF  (OEC.EQ.4)  THEN 
DC  516  I=1*NR 
00  516  J=1,NACT 

516  V(I*NC3.J)  =  B4(I*fJR,U) 

PRINT* * •  V  (83/B4)  IS  • 

NRV  =  NC3  ♦  NR 
ELSE 

PRINT*.*  V  CB3I  I'  • 

NFV  =  NC3 
ENDIF 

CALL  PRM  ( V.NRV. N ACT  ) 

CALL  LSVOF(V»NCQL,  NRV.fi  ACT  *TEN*NCCL.“1»SING*ST'?P.*IER  > 
PRINT*.*  • 

PRINT*.*  LSVDF  G2  IER  =  *»IE= 

PRINT  *  (  // I • 

PRINT*.*  V  OUT  OF  LSVDF  IS  • 

CALL  fRNT  <  V.  fiACT  «  NAC  T  ) 

E3  s  NACT  -  NCV 
IF  (E3.LT.1)  HEN 
DC  517  I  =1.NACT 

517  T2( I ,1 )  =  VCI.NACT) 

E3  =  1 
ELSE 

OC  516  I=1,NACT 
DC  516  J= 1 »E3 
8  T2CI.J)  =  V(  I ,  J*NRV) 

E)  DIF 

PRINT*,*  TRANSFORMATION  MATRIX  T2  • 

CALL  FRNT  (T2.NACT.E3) 

CHECK  TO  SEE  THAT  T2  IS  ORTHOGONAL  TO  B3  AND  B4 

CALL  MMUL  (B3.T2.NC32.NACT.E3 .BCG) 

PRINT*.*  B3  •  T  2  • 

CALL  PRNr (BCG.NC32  ,E3> 

IF  C0EC.EQ.4)  THE*,’ 

CALL  MMUL (B4.T2.NR2.NACT.E3.BCG) 

PRINT*,*  84  *  T2  • 

CALL  PRNT (8CG.NR2.E3) 

ENDIF 
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IF  CDEC.EQ.3)  THEN 

PRINT***  83  SINGULAR  VALUES  * 

ELSE 

PRINT***  B34  SINGULAR  VALUES  • 

Ei'.DIF 

CALL  PRNT  ( SING* NR  V  *11 

CALL  VMULFMCT2*T2»NAC- ,E3 ♦ E3 *NCOL»NCOL, RK* NCOL  ,IE9 1 
CALL  GMINVCE3.E3»rK,KG2  ,J,TAPE) 

CALL  MMUL  CB2*  T2*f.C22  *  .'.-ACT  ,E3  «K082) 

CALL  "MUL  CKOB2  *^G2*NC22»E3*E3»SAr2) 

CALL  VMULFPCS AT  2»  '  2  ,.\C22  *E3*  NACT  *  f.CCL  •  MCOL  *K~B2  ,NCOL  « I  E°  1 
CALL  VMULFPCKCB2»82,.*.C22»NAC7  *NC22  *NCOL*NCOL»SAT2*NCOL*IE°  ) 


SAT2  WILL  BE  SUBSTITUTED  BACK  INTO 
MRIC-SYSTEM  2  FOR  A  NEW  G2. 


NOW  THE  CONTROLLER  GAIN  MATRIX*  G1 


DG  515  I=1*NC2 
DC  519  J=1,NACT 

519  VCI,J)  =  B2CI*NC2,J) 

DO  52C  Is  1 *NC3 

DO  52G  Js 1  *  NACT 

520  V  <  I  +LC2  * J  )  =  B3<I*f.C3*J) 

IF  CDEC.EQ.4)  THEN 

DO  521  1=1 *NR 
DG  521  Js 1 «NACT 

521  VC  I ♦SC2*  'i  C3  »  J )  =  84(1  *NR  ,  J) 

PRINT***  VCB2/B3/B4)  IS  • 

NRV  =  NC2  ♦  NC3  ♦  NR 

ELSE 

PRINT*,*  V  CB2/83)  IS  * 

NRV  =  NC2  ♦  NC3 
E'.DIF 

CALL  PRNT  C V*NR V, NAC" 1 

CALL  LSV3FCV,NCOL*NRV,NACT*TEN,NCOL*-1*SING,STOR.IER) 
PRINT*,*  • 

PRINT*, •  LSVOF  G1  IER  =  »*IE* 

PRINT'C// >* 

PRINT** »  V  OUT  OF  LSVDF  IS  • 

CALL  PRNT  CV,NACT,NAC' > 

E4  =  NACT  -  NRV 
IF  CE4.LT, 11  THEN 
DC  522  Is  1, NACT 

522  T1(I«1)  s  VCItNACT) 

E4  =  1 

ELSE 

00  523  Isl,NACT 
00  523  Js 1 *C4 

523  TlCIvJ)  -  VCI * J*NR  V) 
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END  IF 

PRINT***  TRANSFORMATION  MATRIX  T1  * 

CALL  PRNT  CTl  t'JACT  ,E4  ) 

CHECK  TO  SEE  THAT  T1  IS  ORTHOGONAL  TO  B2»B3,B4 

CALL  MMUL  C92,T1,NC22,NACT,E4,BCG> 

PRINT***  B2  *  T1  • 

CALL  PRNT <BCG,NC22,E4> 

CALL  MMUL  <B3  ♦’  1  **!C32»  *  ACT  *  E4  ,BCG ) 

PRINT*,*  B3  *  » 

CALL  PRNT (BCG *NC3 2  ,£4  ) 

IF  C0EC.EQ.4)  T HE N 

CALL  MMUL  CB4,ti  f  N,»2  »."■  ACT  *E4  *  BCG  ) 

PRINT*,*  B4  *  T1  * 

CALL  PRNT  (BCG *NR2 »E4> 

ENDIF 

IF  (DCC.EQ.3)  THEN 

PRINT*,*  B23  SINGULAR  VALUES  • 

ELSE 

PRINT*,*  B234  SINGULAR  VALUES  • 

ENDIF 

CALL  PRNT  ( SING ,.\R  V » 1 ) 

CALL  VMULFM(Tl,n  ♦NAC*  ,E4,E4  ,NCCL,NC3L»RK,NCCL  ,IE'’  ) 

CALL  GMINV(E4,E4,-K,RG1 ,J,TAPE) 

CALL  MMUL  (Bl,Tl,f.C12,:.ACT,E4,KC8l> 

CALL  MMUL (K0B1,=GI,NCI2,C4*E4,SAT> 

CALL  VMULFP(3AT,ri,NC12,E4,NACT,NC0L,':C0L,K0Bl,NCCL,:ER) 
CALL  VMULFP(KCBl»Bl»\Ci2»NACT,fiC12»NCC'-»NCOL»SAT*r,COL,IEi:  > 


SAT1  MILL  BE  SUBSTITUTED  BACK  INTO  MRIC- 
SYSTEH  1  FOR  A  NEW  GI 


ZZ  =  1 
GOTO  115 
&00  CONTINUE 
ZZ  =  0 


THE  PROBLEM  IS  NOW  COMPLETE 


PPINT*<///)» 

PRINT*, •  THIS  RUN  HAS  BEEN  COMPLETED  * 
PRINT  *  C///1 • 


FROM  HERE  UE  CAfJ  START  OVER,  REARRANGE,  OR  STOP 
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PRINT*#*  ENTER  1  TO  CHANGE  THE  WEIGHTING  MATP I X  * 

Pc IN' •# •  * 

PRINT*. •  ENTER  2  TO  MAKE  A  FOUR  CONTROLLER  RUN  * 

PPIN**#*  * 

PRINT*#*  ENTER  2  ALSC#  TO  REARRANGE  MODES  FOR  *1  =  **N 

Pc INT  *  #  *  * 

PRINT*#  •  ENTER  3  TO  TERMINATE  THIS  JOB  • 

READtc  »  *)  Q 
PRINT*!//  )• 

PRINT*. Q 
IF  tQ.EO.l)  THEN 
GOTO  IOC 

ELSE  IF  ( Q.EQ.2 )  T HES 

GOTO  20 

ENDIF 

END 


Ill 


SUBROUTINE  FACTOR (N,A  ,S  .MR > 

C  A=5*S 

DIMENSION  A  <  1  )  ,$U  ) 

CQMNON/NAINB/NCOL ,NC0L1 
COHHON/IN  OUT/KOUT 
T0L=l.E-6 
MR  =  0 

NN=N*NC OL 
T  OL 1=0 • 

00  1  I=1,NN,NC0L1 
R=ABS(ACI  )) 

1  IF  CR . GT • T0L1  )  T0L1=R 

T0L1=T0L1 *1.E-12 
11  =  1 

DO  50  1=1  ,N 
IM1  =  I-1 

00  5  JU= I ,NN  ,NCOL 
5  S(JU>=0. 

10= II+IM1 

r  =a  <  r  o )  -o  ot  ( i  mi  ,  s  i  r  r  > .  s  u  i ) » 

IF  (A8SIR >.LT.(TOL»A(IO)*TOL1>>  GO  TO  50 
IF  <fl>  15  ,50  *  20 
15  MR  =  — 1 

WRITEIKlIUT  ilOOO) 

1000  FORMAT (37HQTRIED  TO  FACTOR  AN  INDEFINITE  MATRIX  > 
RETURN 

20  SUD)=SQRT(RI 
HR  =  MR  ♦  l 

IF  (I.EQ.N)  RETURN 
L  =  1 1 *NCOL 

00  25  JJ=L«NN,NCOL 
Ivl=  JJ+IMI 

25  S (IU)=(A( I J ) -DOT ( IM1 »SCII ) tSC JJ) ) )/StlO> 

50  II  =  II*NCOL 

RETURN 
END 


SUBROUTINE  FORMA ( A . 0 1 W, N »N2 , I C  > 

COMMON/MA 1 NR/NCOL 

REAL  A<NCCL,NCOL> ,W(17) ,0< 17) 

INTEGER  IC(N)»I»J»N»M 

DO  1  1=1, N2 

00  1  J=1.N2 

A ( I ,J)=0.0 

1  CONTINUE 
00  2  I = 1 »  N 
M  =  IC(I) 

A ( C I *N) , < I ♦N) ) =01 M  ) 

All  »  (  I  *  N  )  )  =  1.0 
A(d>N),I  )  =  -(U(M)**2) 

2  CONTINUE 
RETURN 
ENO 
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SCcKo-jr  fnt  F0ftMHt8#PHl  ,  N  .  N2  *  NACT  »  I  C  ) 

C0MM0N/MAIN8/NC0L 

REAL  8(NC0L»NC0L) .FHI (NCOLtNCOL) 

INTEGER  IC(N) ,NACT,N,M,I,J 

00  1  l=l»N2 

DO  1  J=1,NACT 

B(ItJ)  =  0.0 

CONTINUE 

00  2  1=1. N 

M  =  I C  <  I  > 

00  2  J=1»NACT 

B( (N*I ) ,J  )  =  PHI <M,J> 

CONTINUE 

RETURN 

END 


SUBROUT  INE  FORMC(C,PHIS  »N  *  N2  »NS  EN  » I  C  ) 
COM MON/HA  I NB/NCOL 

REAL  C(NCOLfNCOL) *  PHIS ( NCOL  » NCOL) 

INTEGER  IC(N) ,M,NSEN,N,N2, I »J 

DO  1  I =1 » NSEN 

DO  1  J=1«N2 

C(I  *J)  =  0.0 

CONTINUE 

00  2  1=1, NSEN 

00  2  J=1»N 

M  =  IC<J) 

CCI.JI  =  PHIS(Mfl) 

CONTINUE 

RETURN 

ENO 


SUBROUTINE  FORMQt Q , A ,N , IC ) 

COMMON/NAI N8/NC0L 

REAL  A(NCOL ) , 0 (NCOL »NCOL ) 

INTEGER  I  . J,K ,M,N,N2, IC (NCOL  ) 

N2  =  N  *  2 

DO  l  I  =1 ,  N2 

DO  I  J= 1 »  N2 

Q ( I  * J)  =  0.0 

CONTINUE 

DO  2  1=1, N 

J  =  C 

H  =  IC<  I) 

00  5  K=I-1,I 
0(1 ♦K , J>K  )  =  A ( M) 

CONTINUE 

RETURN 

ENO 


SUBROUTINE, GHINV(NR»NC»A.U»MR»MT) 
DIMENSION  A  (  1  )  tU (  l  ) 

COMMON  /fiA  (%N1  /  ND  i  M  «  NO  l  Ml  »S  (  1  ) 
COMMON/MA 1NB/NCOL .NCOLl 
COMMON/ IN COT/KOUT 
TOL=l.E-l 2 
MR =NC 

nrmi=nr-l 
TOL 1=1 «E"20 
JJ=1 

DO  100  J= 1 »NC 
F  AC  =00 T (NR  »  A  (  J  J ) .  A ( JJ  > ) 

JM1 =J-1 
JRM=J  J+NRM1 
JCM=JJ*JM1 
00  20  I=J J« JCM 
20  U(I>=0. 

U(JCM)=i.O 

(F  ( J.EQ. 1 )  SO  TO  54 
KK=1 

DO  30  K=l »  JM1 

IF  (S(K) . EQ.l  .0 )  GO  TO  30 

TEMP=-DCT (NR. A( JJ) , A (KK) ) 

CALL  VADD  (K  .  T EMP .  II  ( JJ >  .U(KK)  ) 

30  KK=KK>NCOL 
00  50  L=1 *2 
KK=  1 

DO  50  K  =  1  .  JML 
IF  (S(K).EQ.O.)  GO  TO  50 
TEHP=-OOT (NR,A( JJ) • A ( KK ) ) 

CALL  V  ADD (NR.TEMF  »  A( JJ ) *A ( KK ) ) 
CALL  V4CD(KfTEMP,U(JJ),U(KK>> 

50  KK=KX*NCOL 

TOL1=TOL*FAC 

FAC=DOT (NR.A(JJ).A(JJ)) 

54  IF  (FAC.GT.TOLl )  GO  TO  70 
DO  55  I  =  J  J.  JRM 

55  A ( I )  =0  • 

S (J)=0. 

KK  =  1 

DO  65  K=l .JP1 
IF  (S(K).EQ.O.)  GO  TO  65 
TEMP=-OOT (K.U(KK) . U ( J J ) ) 

CALL  VAOD( NR.TEMP »A(JJ> »A ( KK ) ) 

65  KK  =  KK  *NCO  L 

FAC=DOT ( J 'U(JJ) »U( JJ) ) 

MR  =  MR-1 
GO  70  75 
TO  S(J)=1.0 
KK=  1 

DO  72  K=l . JM1 
IF  (S(K).EQ.l.)  GO  TO  72 
TEHP  =  -OOT  <NR*iA(  JJ).A(KK)) 

CALL  VADO^IK  »T  EMP»  U  ( JJ  )  «U(KK)  > 
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72  KK  =  KK  *WCOL 

75  FAC=l , /SORT (FAC) 

DO  6C  I=JJtJRH 
80  A ( I )  =A ( I ) *F  AC 

DO  85  I=JJ»JCM 
85  L»(I  >=U<  1)  *FAC 
100  JJ=JJ+NCOL 

IF  (MR.EQ .NR.OR.MR.EB.HC>  GO  TO  120 
IF  (flT.NE.O)  WRITE  (KOUT  *110  )NR  »NC  »  MR 
110  FORMAT  <  13  ,1HX  »I2»8H  M  RANK, 12) 

120  NEND  =  NC*N  COL 
JJ=  1 

DO  135  J=1»NC 
DO  125  I  - 1 , NR 
t 1=1-0 
S (I  )=0. 

00  125  KK =JJ, MEND. NCOL 
125  S(I):$U)  ♦ACII>KK>*DCKKI 

1I  =  J 

DO  130  1=1, NR 
U(II)=S(I ) 

130  I  1=  1 1  ♦NCOL 
135  J J= JJ*NCO  LI 
RETURN 
END 


SUBROUTINE  INI  EC <N , A ,C  ,S, T  ) 
S=INTEGRAL  EA*C*EA  FROM  0  TO  T 
C  IS  DESTROYED 
DIMENSION  A< 1 )  ,C(l  )  ,S (l  ) 

COMMON/MA INI /NDIM,NOi Ml  , X ( 1  ) 
CQMHOM/MA  I  NB/NCOL  ,  NCOH 
COMMON/MA I N2/C0EF  C 100 ) 

NN=N*NCOL 
NM1 =N-l 
I  ND=0 

ANORM=XNOftH(N , A) 

DT  =  T 

5  IF  (ANORM*A8S(DT).LE«0.5)  GO  TO  10 
DT=DT/2. 

I N0=IND+1 
GO  TO  5 

10  DO  15  1=1 ,NN»NC0L 
J=I ♦NM1 
DO  15  JJ=I*J 

15  S< JJ)=OT*C(JJ) 

Tl=DT**2/2. 

DO  25  IT=3 ,15 

CALL  MMUL <A,C,N,N,N,X) 

DO  20  1  =  1  ,N 
I  I  =  ( I -1 > »  NCOL 
00  2Q  JJ= I ,NN«NCOL 
11=11*1 


C< JJ>  =( X( JJ)*X(II ) >*T1 
20  S  (  JJ>  -S  ( J  J)  +.C  ( JJ) 

25  T t =0T/FL0AT( IT) 

IF  (INO.CQ.O)  60  TO  100 
COEF (1 1 )- 1 aO 
00  30  1=1  .10 
11=11-1 

30  COEFI II >=DT*COEF( 11*1 ) /FLOAT (I > 
11  =  1 

00  40  1=1 .NN.NCOL 
J  =  I  ♦NMl 
00  35  JJ=I.J 

35  X <JJ)=A<JJI*C0EF<1> 

XU1>=X<I  IM-C0EF(2) 

40  I  I  =  1 1 ♦NCO  LI 
00  55  L=3  .11 
CALL  MHUL (A.X  .N.N.N.C) 

11  =  1 

T l  =  C0EF <L  ) 

DO  55  1=1. NN.NCOL 
J=I ♦NMl 
00  SO  JJ=  I  ,J 
50  X(JJ)=C(JJ) 

XUI)=X<I1)*T1 
55  I  I  =1 1  *NC0L1 
C  X=EXP<A*OT) 

L  =0 

60  L=L*1 

CALL  MMUL (X .S .N.N ,W»C ) 

11=1 

DO  90  1=1  ,N 
J=II 

IF  (I.EQ.l)  60  TO  75 
DO  70  JJ=  I  .II  .NCOL 
S  < JJ)=S( J  ) 

70  J=J>1 

75  00  85  JJ= I  .N 

KK=  JJ 

DO  80  K=l. NN.NCOL 
S(J)=S( J)*C(K)*X(KK) 

80  KK=KK*NCOL 
85  J=J>NCOL 

00  87  JJ= It NN.NCOL 
87  C(JJ>=X(JJ) 

90  1 1  =  1 1 >NCO  L 

IF  (L.EQ. INO)  60  TO  100 
CALL  MHUL <C.C.N»N.N»X> 

GO  TO  60 
100  CONTINUE 
RETURN 
ENO 
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SUBROUT  INE  MLI NEQ  (  N  » A  ,  C  »X  ,  TOL  » I ER  ) 
SOLVES  /“  •  X*X  *  *C=C 
A  AMD  X  CAM  BE  IN  SAKE  LOCATION 
ANSWER  RETURNED  IN  C  AND  X 
DIMENSION  A<l».C<l>,)m> 

COMMON/MA I N8/NC0L  tNCOLl 
COMHON/MA IN3/F<1> 

AOV=TOL»l *E-6 
DTs.5 
DT1=0. 

NN=N*N COL 

DO  5  11  =  1  *  NN • NCOL1 
5  OT 1 =OT 1-A ( 1 1 ) 

oti=oti/n 

IF  (OT1.GT.4.0)  OT=OT»4.0/DT1 
11=1 

DO  20  1=1 »N 
DO  15  J J= I «NN  »NCOL 
15  X<JJ)=DT*A(JJ) 

X(II)=X(II)-.5 
20  II=II*NC0L1 

CALL  GM  IN  V  (  N  *  N  f  X  «  F  »MR  ,  0  ) 

I  £R=4 

IF  (MR.NE.N)  RETURN 
CALL  MMUL  <C  t  F  »N  »N  *N ,X  ) 

C  INITIALIZATION  OF  X»F 

1=1 

DO  40  1 1= l  .NN.NCOL 
J=I1 

IF  (  l  «  EG*  1  )  GO  TO  30 
DO  25  JJ= I fll .NCOL 
CCU)  =  C( JJ  » 

25  J=J*1 

30  IO=J 

DO  35  JJ= II « NNtNCOL 
C( J)  =  DT*OOT (N.F<I I  liX( JJ) 1 
35  J=J*l 

F  <  J  D>  =F  <1  D)  >1  *0 
40  1=1*1 

DO  90  IT=1»20 
NEZ  =0 

CALL  HMUL (C»F«NtN»N»X } 

1=1 

II  =  1 
J  =  1 

60  TO  70 
60  J=I  I 

DO  65  JU= I » 1 1  iNCOL 
C(J)=C(JJ) 

65  J=J*1 

70  IO=J 

DTI =C( J) 

00  75  JJ= II iNN*NCOL 
C(J)=C( J) *00T (N»F(I1) tX(JJ) > 
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80  x(JO)=F(JJ> 

IF  ( ABS  CC  (ID)  ).6T  .  1  .E150 )  60  TO  95 

IF  (AeS<C(IO>-OTl>.LT-(ADv^TOL*ABs<C(ID)))>  NEZ=NEZ-1 
1  =  1-1 

1 1-II-NCO  L 
IF  (l.LE.N)  60  TO  60 
IF  <MEZ.£Q-N)  GO  TO  150 
CALL  MHUL (X  »X  *N,N  ,N  »F  ) 

90  CONTINUE 
95  I£R  =  1 

RETURN 

150  CONTINUE 
NH1=N-1 

DO  155  I - 1 • NN  *NCOL 
1 1 ~  X ♦NMl 
00  155  JJ  =  I,II 
155  X(Jd»  sCisJJi 
1ER-0 
RETURN 
END 


SUBROUTINE  NMUL (X  1 7  ,Nl  *N2  ,  N3  ,Z  > 
COHMON/NA1NO/NCOL 

DIMENSION  X<NC0L«1  )*Y(NCOL«U  tZ(NCOLtl) 
DO  3  J=1,N3 
DO  2  1=1, N1 
5  =0  • 

DO  1  K=1,N2 

1  S=S-X(I »K  >*Y(K* J» 

2  ZU*J>=5 

3  CONTINUE 
END 


SUBROUTINE  NRlC(N»A'S'QfX'Z'TOL'lER> 
DIMENSION  A(1 »  ,S( 1  )  ,Q(1  )  ,X( 1  )  ,ZC1 ) 
COMMON/MAI  N1/N0IM.ND1M1  ,FU  ) 
COMMON/MA I Nfl/NCOL ,NC0L1 
COMMON /MA I N2/TR (1  ) 

COMMON/ IN  OUT /KOUT 

A0V=T0L*1 .E-6 

NN=N* WOOL 

NM1=M-1 

IN0=1 

COUNT =0. 

IF  (lER.CQ.L)  COUNT =99. 

IF  (ICR.EQ.l )  MR=N 
IF  < ICR.EO* i >  60  TO  100 
7  la-1, 

300  CONTINUE 


non 


F 


1  Eft  =0 

COUNT=COUNT>1 . 

OC  15  I =1 tN 
00  15  J  =  I  «NH  tNCOL 

15  X(J)=-S(J> 

CALL  1NTE6(N,A.X,2,T1 ) 

CALL  FACTOR  (N*2»X* MR) 

I£R=1 

IF  (HR.LT.O)  60  TO  200 
IER  =  0 

CALL  GMINV(N*NtX»2*MRtO) 

CALL  TFR( TR.Z.N.N,1 .2 > 

CALL  HMUL (2«TRtN»N*N.X) 

00  18  11=1 .NN.NCOLl 

I=II 

00  17  J=I I «NN  *NCOL 
X(J)=CX(J)>X(I)>/2. 

XC1 )=X(J) 

17  I  =I>1 
10  CONTINUE 

100  CONTINUE 

DO  16  1=1 *N 

16  TR<I)=-1.0 
A+SX  IS  STABLE 

POSSIBLE  UNCONTROLLABILITY  IF  MR.NE.N 
Jin  OIL  LOW  IS  A  NUTTY  MATH  PROF 
T0L1=T0L/10. 

MAX I  T=AO 

00  AO  I  T=  1 1  HA  XIT 
IF  (IER.EQ.l)  GO  TO  101 
CALL  MHUL (S«X tN«H <N»F) 

CALL  MHUL (X.F.N.NtN.Z) 

DO  20  I =1 (NNtNCOL 
l I=I*NH1 
DO  20  J=I til 
X(J)=A(J)-FCJ) 

20  2C J)=Z< J> ♦QtJ> 

101  CONTINUE 
IER=0 

CALL  MLINEQ(N*X*Z *X«T0L1 tIER) 

IF  (IER.NE.O)  60  TO  200 
L  =0 

C1=0.0 

11=1 

DO  25  1=1 »N 

IF  <ABS(X  (IJ)-*TRCI  ))#LT«(AOV»TOL*X(I1)  ))  L=L^1 
TRC I )=X (I  I ) 

11=11 ♦NCOLl 
25  C1  =  C1>TR<D 

IF  (ABS(Cl) .GT.1.E20)  GO  TO  50 
IF  CL .NE.N)  GO  TO  AO 
CALL  GMINVCN»N»2»F*MR»0) 

CALL  HMUL (S ♦ X » N ,N * Ni 2 > 

DO  30  1=1  »NN » MCOL 
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IXs J ♦WMl 
DO  30  J  =  I  ,11 
30  Z(J)=A< J»-Z< J) 

IF  (MR.NE.N)  WRITECKOUT »35)MR 

35  FORMAT<27HORICCATI  SOLN-IS  PSD—  RANK  ,33) 

SO  TO  65 
40  CONTINUE 

WRl TE ( K0U7  » 45  )  MAXlT 

45  FORMAT (27H0RICCATI  NON- CONVERGENT  IN  ITERATIONS) 

GO  TO  60 

50  WRIT£(KOOT«55)IT»T1 

55  FORMAT (30H0RICCATI  BLOW-UP  AT  ITERATION  ,I2»12H  INITIAL  T=  ,F10.5) 
60  I ER  =  L 

65  RETURN 

20 0  IF  (IND.EQ.2)  GO  TO  250 
IF  (COUNT .GE .10 • )  RETURN 
Tl=Tl/<2«* ‘COUNT ) 

1MD=2 
GO  TO  300 

250  Tl=Tl*<2«*  *C0UN7) 

IND=1 
GO  TO  300 
ENO 


SUBROUTINE  PRNT (MAT  ,N ,  N) 
C0MM0N/MA1NB/NC  OL 
REAL  MAT( NCOLtNCOL) 

INTEGER  N ,t  * 0 « K,M 
PRINT*,*  * 

IF  (M.GT.12)  GOTO  2 
DO  1  1=1, N 

PRINT*(lX,l2F10.4>*,(MAT(I*J)*J=l,h) 

1  CONTINUE 
GOTO  10 

2  CONTINUE 

IF  (M.GT.24)  THEN 
CALL  PRNT  XL( MAT,N,M) 

RETURN 

ENOIF 

00  3  1  =  1,  N 

PRINT* (IX ,12F10.4)*,(MAT(I»J) ,3=1,12) 

3  CONTINUE 
PRINT •(//  )» 

00  4  1=1, N 

PRINT  *  <  IX  .12F10.4)  *,{MAT(I  »3),J*13,*4) 

4  CONTINUE 

10  PRINT* (///)* 

RETURN 
ENO 


SUBROUTINE  PRNT  XL  (  MAT  »N  »M ) 

COMMON/MA  TN^'KiCCL 
COMMCN/MA l NA/NDA 
REAL  MAT(NOA.NOA) 

INTE6ER  1 » J.K.L.M.N 
PRINT*. •  * 

DO  1  L  =  1 .  M.12 
K  :  L  Ml 

IF  <M-L.LT.ll >  K  =  M 
00  2  1=1* N 

PRINT  MIX  .12F10.4)  • ,  <  MAT  C I  .  J  )  ,  J=L  .K  ) 

CONTINUE 

PRINT  M //  )  • 

CONTINUE 
PRINT  *(///) • 

RETURN 

END 


SUBROUTINE  TFR<X,A.N,M,K,I > 

1  GIVES  X  =  A.  2  GIVES  X  =  A  TRANSPOSE 

3  GIVES  X  =  A  AS  A  VECTOR 

4  GIVES  A  =  X  WHERE  X  WAS  A  VECTOR 

01 MENS I  ON  X(l).A(l) 

C OMMO M / MA I NB / NC CL 

JS=<K-1  )*NCOL*M 

JENO=M*NCOL 

GO  TO  <10  ,30.50.70  ).l 

00  20  11=1* N 

DO  20  JJ= II , JENO.NCOL 

X  (JJ)=A(JJ*-J5  ) 

RETURN 

00  40  11-1 tH 
KK=<II-1>  »  NCOL 
00  *0  JJ=1.M 
LL=< JJ-1> *NCOL+II 
X(KK*JJ>= A(LL*JS> 

RETURN 

KK=0 

00  60  11=1 .JENO.NCOL 
LL=S1*N-1 
DO  60  J J= I I «LL 
KK=KKM 

X(KK)=A<J J+JS> 

RETURN 
KK=M*NM 
00  80  11=1. M 
LLMM-IIJ *NC0L»1 
00  80  X J=  1  *N 
KK=KK-1 
JJ=LL*N-I J 
A  <  JJ>  JS)  =  X<  KK  ) 

RETURN 

END 


«• 


F  UHC  H  CN  001  (NR  ,A  ,8  ) 

DIMENSION  A ( 1  )  ,8(1  ) 
0OT=0. 

DO  1  1=1, NR 
1  DOT -DOT  ♦A(I)*0(I) 

RETURN 
END 


SUBROUTINE  VAOO(N,C1  ,A*B> 
DIMENSION  A(1 )  ,8(1  ) 

DO  1  1  =  1, N 

1  A(I)=A(I> +C1*B(I) 

RETURN 

END 


FUNCTION  XNORM(N.A) 

C  COMPUTES  AN  APPROXIMATION  TO  NORM  OF  A  --  NOT 
DIMENSION  A(l) 

COMMON/MAI NB/NCOL , NC0L1 

NN=N*NCOL 

C1  =  0. 

TR=A(1» 

IF  (N.E6.1)  GO  TO  20 
1-2 

DO  10  11=  NC0L1  ,NN , NCOL 
J  =  II 

DO  5  JJ=1 ,11, NCOL 
Cl=Ct  »A8S (A ( J ) »  A  (  d  J ) ) 

5  d  =  J*l 

TR=TR*A(J> 

10  1=1*1 

TR=TR/FLOAT(N) 

DO  IS  11=1 ,NN ,NC0L1 
15  Cl=Cl*(A( II »-TR»* *2 
20  XN0RM=ABS(TR)*SQRT(C1) 

RETURN 

END 


BOUND 


123 


Appendix  E 

Line-of-Sight  and  Defocus  Algorithm 
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The  equations  relating  optical  surface  motion  are  given  below  in 
Eqs  E-2,  E-3,  and  E-4.  These  form  the  line-of-sight  and  defocus  cri¬ 
teria  as 

LOSX  -  Y/F,  LOSY  =  X/F,  DEFOCUS  =  Z  (E-l) 

where 


X 

*  A. 

-X 

+  X_  -  R  *0Y  +  A. 

*  0Y 

-  2T  '  0Y  |  +  X  -  X- 

(E-2) 

1 

L  p 

t  p  p  2 

s 

tj  t  f 

Y 

■  A. 

-Y 

+  Y_  R  • 0X  -  A_ 

*  0X 

+  2T  *  0Xj  +  Y  -  Y, 

(E-3) 

1 

P 

t  -  p  p  2 

s 

tj  t  f 

Z 

»  A_ 

Z 

-  2  Z  +  Z  "1  +  Z  - 

Z  , 

(E-4) 

3 

P 

s  tj  t 

f 

F  *  8.051  ”  focal  length 

where 

A1  -  0.2987 

A2  -  93.90 

A3  -  0.0892 

R  *  53.9 
P 

T  -  66.95 

The  terms  X^,  Y^,  Z^,  QX^,  0Y^,  for  I  ■  p,  s,  t,  f  refer  to  the  trans¬ 
lations  and  rotations  in  the  global  X,  Y,  and  Z  directions  of  the  primary 
(p),  secondary  (s),  tertiary  (t)  and  focal  plane(f)  elements.  The  coeffi¬ 
cients  Aj,  T,  and  F  are  functions  of  the  radius  of  curvature  of  the  mirrors. 
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These  are  given  by 


Rp  -  \  +  2“l  +  '2>] 


R 


F  = 


R  R 

tp 


4T  -  2R^ 


T  = 


+  ^  +  t2 


where 


( E— 5 ) 

(E-6) 

(E-7) 

(E-8) 

(E-9) 


Rp  -  radius  of  curvature  of  the  primary  mirror 


Rt  =  radius  of  curvature  of  the  tertiary  mirror 


t^  *  axial  distance  from  primary  to  secondary  mirrors 
t£  »  axial  distance  from  secondary  to  tertiary  mirrors 
The  expressions  for  the  translation  and  rotation  of  each  mirror  in  global 


coordinates  may  be  formed  in  terms  of  the  displacements  at  the  nodes  which 
support  the  mirrors.  These  are  given  by  the  following  equations  in  which  the 


numerical  subscripts  indicate  specific  support  nodes. 
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Primary  Mirror: 


XP- 

X34  +  1.25  ft34  -  Y 

35^ 

(E-10) 

Y  - 
P 

0.50  (Y34  +  Y35) 

(E-ll) 

Z  - 

P 

-0.2143  (Z34  +  Z35) 

+  0.7143 

(Z28  +  Z30) 

(E-12) 

0X 

P 

*  0.0714  (Z34  +  Z35) 

+  0.0714 

(Z28  +  Z30) 

(E-13) 

0Y 

P 

-  0.125  (Z34  -  Z35) 

(E-14) 

ez 

p 

*  0.125  (Y35  -  Yj4) 

(E-15) 

Secondary  Mirror: 


Z 

s 


(E-16) 

(E-17) 

(E-18) 


ex  - 

s 

0X4O 

(E-19) 

0Y  - 
3 

0Y4O 

(E-20) 

ez  =■ 

s 

0Z4O 

(E-21) 

Teritiary  Mirror: 


Xt  -  X27  +  0.3750  (Y29  -  Y27)  (E-22) 
Yt  -  0.50  (Y27  +  Y29)  (E-23) 
Zt  -  0.7143  (Z2?  +  Z29)  -  0.2143  (Z32  +  Z^)  (E-24) 
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